In this study we explored the mode of genetic inheritance of Varroa destructor, a parasitic mite of honeybees, so far known as a haplodiploid species with a arrhenotokous parthenogenesis reproductive mode. The mite shows contrasting phenomena of highly inbreeding life style (sib-mating), yet maintaining relatively high genetic variation.
The experimental setup consisted by a three-generational pedigree
with sib-mating.
Sample size: 30 families, total of 223 individuals 3 different colonies,
from OIST apiary
For details of mite collection and pedigree construction, please see the
original manuscript.
library("tidyverse")
library("plyr")
library("dplyr")
library("ggplot2")
library("scales")
library("ggpubr")
library("gridExtra")
library("grid")
library("GGally")
library("vcfR") # for extracting genotype data from a vcf file
library("data.table")
library("stringr")
library("janitor")
library("gmodels")
library("rstatix")
library("freqtables")
library("broom")
library("DescTools") # for the Goodness-of-Fit test, 3 variables; and for Breslow-Day Test for Homogeneity of the Odds Ratios
library("vcd") # for the woolf test testign log odds for each pair
library("patchwork") # for gathering the plots
#library("fuzzyjoin") # to join tables based on a string in a column
#library("aspi") # Repeated G–tests of Goodness-of-Fit, work only for 2 variables..
#library("RVAideMemoire") # Repeated G–tests of Goodness-of-Fit, work only for 2 variables...
#library("InfoTrad")
#library("ggthemes") # for more colors in the ggplot
knitr::opts_chunk$set(echo = TRUE, warning = FALSE, message = FALSE,
fig.width = 8,
fig.asp = 0.6,
out.width = "100%")
#fig.width = 6,fig.asp = 0.8,out.width = "100%"
Extract genotypes for each site and individual. The metadata for all samples can be found in here.
vcf <- read.vcfR("/Users/nuriteliash/Documents/GitHub/varroa-linkage-map/data/vcf_filter/Q40BIALLDP16HDP40mis.5Chr7/Q40BIALLDP16HDP40mis.5Chr7.recode.vcf", verbose = FALSE )
vcf
## ***** Object of Class vcfR *****
## 223 samples
## 7 CHROMs
## 35,169 variants
## Object size: 293.5 Mb
## 0 percent missing data
## ***** ***** *****
# extract the genotype for each site in each individual
gt <- extract.gt(vcf, element = "GT")
gt <- as.data.frame(t(gt)) %>%
rownames_to_column("ID")
#clean the ID names
gt$ID <- sub("_[^_]+$", "", gt$ID)
The individual ID name nomenclature is composed of 3 parts, separated
by an underscore (_).
The first part is the family ID, the second is the unique name of the
individual, and the the third part indicates its generation, sex, and
its position in relation to the foundress mite (generation F0).
For example:
Individual ID 240_240a_son, belong to family 240, has
a unique name of 240a and is the son of the foundress
mite of family 240, that is, its a male of F1 generation.
Individual ID 240_241b_grndat, belong to family 240,
has a unique name of 241b and is the granddaughter of
the foundress mite of family 240, that is, its a female of F2 generation
(and the daughter of 240_240a_son).
From the 223 individuals of 30 families, we excluded non adults mites
(for which its hard to determine the sex).
Eventually, we kept 144 individuals (adult males and females in F0, F1
and F2 generations) of 26 families, and all 35,169 biallelic sites.
For example, these are the first 6 sites (13565-25513) of chromosome (NW_019211454.1), in one family (family number 240).
table <- gt %>%
t() %>%
as.data.frame() %>%
row_to_names(row_number = 1) %>%
dplyr::select(contains(c("son", "dat", "fnd"))) # keep only adults of F0, F1 and F2
table %>% select(contains(c("240_"))) %>% head()
## 240_240a_son 240_241c_grnson 240_241b_grndat
## NW_019211454.1_13565 0/1 <NA> <NA>
## NW_019211454.1_18037 0/0 <NA> <NA>
## NW_019211454.1_20998 <NA> <NA> <NA>
## NW_019211454.1_24935 <NA> 0/0 <NA>
## NW_019211454.1_25470 <NA> <NA> <NA>
## NW_019211454.1_25513 <NA> <NA> <NA>
## 240_241a_grndat 240_241_dat 240_240_fnd
## NW_019211454.1_13565 0/1 0/1 0/0
## NW_019211454.1_18037 0/0 0/0 0/0
## NW_019211454.1_20998 0/0 0/0 0/0
## NW_019211454.1_24935 0/0 0/0 0/0
## NW_019211454.1_25470 0/0 0/0 0/0
## NW_019211454.1_25513 0/0 <NA> 0/0
The family members include:
Each individual, was genotyped for each site with one of the three genotypes:
For more information about the mapping, variant calling and variant filtration, please see the Snakemake pipeline in here.
In the following code we viewed the F2 offspring genotype of all possible nine crosses between F1 male and females.
Then we cross heterozygotic sites, to explore the mode of genetic inheritance in varroa mite:
For each crossing we have expected genotype-ratio of the F2, based on the following assumptions under arrhenotokous haplodiploid system:
To detect false sites
we expect all F2 offspring to be homozygotic (0/0) like their parents (F1)
x <- dat %>%
dplyr::select(c("sample","gt","sex", "exp","prop", "total")) %>%
mutate(type = "exp") %>%
dplyr::rename(count = exp)
y <- dat %>%
dplyr::select(c("sample","gt","sex", "obs","total")) %>%
mutate(type = "obs") %>%
dplyr::rename(count = obs) %>%
mutate(prop = (count/total))
df_male <- rbind(x, y) %>%
dplyr::filter(sex == "male")
df_fem <- rbind(x, y) %>%
dplyr::filter(sex == "female")
p_male_pooled <- ggplot(df_male, aes(fill=gt, y=prop, x=type)) +
geom_bar(position="fill", stat="identity", ) +
xlab("F2 males") +
ylab("Proportion of F2 genotype") +
labs(title = "F2 Males of F1 male (0/0) x female (0/0)",
subtitle = paste0("Pooled sites = ", df_male %>% dplyr::filter(gt == "0/1") %>% dplyr::filter(type == "obs") %>% summarise(sum(total)) %>% as.numeric())) +
scale_x_discrete("type", labels = c("Expected","Observed")) +
labs(fill = "Genotype") +
theme_classic() +
theme(axis.text.x = element_text(size=15),
axis.title.x = element_blank())
p_female_pooled <- ggplot(df_fem, aes(fill=gt, y=prop, x=type)) +
geom_bar(position="fill", stat="identity", ) +
xlab("F2 females") +
ylab("Proportion of F2 genotype") +
labs(title = "F2 females of F1 male (0/0) x female (0/0)",
subtitle = paste0("Pooled sites = ", df_fem %>% dplyr::filter(gt == "0/1") %>% dplyr::filter(type == "obs") %>% summarise(sum(total)) %>% as.numeric())) +
scale_x_discrete("type", labels = c("Expected","Observed")) +
labs(fill = "Genotype") +
theme_classic() +
theme(axis.text.x = element_text(size=15),
axis.title.x = element_blank())
grid.arrange(p_female_pooled, p_male_pooled, ncol = 2)
# plot ratio for each individual:
# males:
ggplot(df_male, aes(fill=gt, y=prop, x=type)) +
geom_bar(position="fill", stat="identity", ) +
xlab("F2 siblings") +
ylab("Proportion of F2 genotype") +
labs(title = "F2 Males of F1 male (0/0) x female (0/0)") +
# geom_text(aes(label = n), position = position_stack(vjust = 0.5)) +
geom_text(data = filter(df_male, gt == "1/1"), aes(y=prop, x=type, label=total), vjust = 0) +
facet_wrap(~ sample) +
labs(fill = "Genotype") +
theme_classic() +
theme(axis.text.x = element_text(size=15),
axis.title.x = element_blank())
# and females:
ggplot(df_fem, aes(fill=gt, y=prop, x=type)) +
geom_bar(position="fill", stat="identity", ) +
xlab("F2 siblings") +
ylab("Proportion of F2 genotype") +
labs(title = "F2 Females of F1 male (0/0) x female (0/0)") +
# geom_text(aes(label = n), position = position_stack(vjust = 0.5)) +
geom_text(data = filter(df_fem, gt == "1/1"), aes(y=prop, x=type, label=total), vjust = 0) +
facet_wrap(~ sample) +
labs(fill = "Genotype") +
theme_classic() +
theme(axis.text.x = element_text(size=15),
axis.title.x = element_blank())
we expect the F0 foundress to be homozygotic (0/0), like her son
# plot the pooled data for all F0 foundress with male (0/0) x female (0/0) offspring
ggplot(observed, aes(fill=gt, y=prop, x = type)) +
geom_bar(position="fill", stat="identity") +
ylab("Proportion of F0 genotype") +
labs(title = "F0 foundress of F1 male (0/0) x female (0/0)",
subtitle = paste0("Pooled sites = ", observed %>% dplyr::filter(gt == "0/1") %>% summarise(sum(total)) %>% as.numeric())) +
theme_classic() +
theme(axis.text.x = element_blank(), axis.title.x = element_blank()) +
labs(fill = "Genotype")
# plot for each individual:
ggplot(observed, aes(fill=gt, y=prop, x=type)) +
geom_bar(position="fill", stat="identity", ) +
xlab("Observed") +
ylab("Proportion of F0 genotype") +
labs(title = "F0 foundress of F1 male (0/0) x female (0/0)") +
geom_text(data = filter(observed, gt == "1/1"), aes(y=prop, x=type, label=total), vjust = 0) +
facet_wrap(~ sample) +
labs(fill = "Genotype") +
theme_classic() +
theme(axis.text.x = element_blank(), axis.title.x = element_blank())
we expect all F2 offspring to be homozygotic (1/1) like their parents (F1)
# set the families (include only families with at least one adult F2)
family = grep("grndat|grnson",gt$ID, value=TRUE) %>%
str_extract("[^_]+") %>%
unique()
# define a list to put all the data frames in
obs <- list()
for (fam in family) {
obs[[fam]] <- table %>%
dplyr::select(starts_with(fam)) %>%
dplyr::filter_at(vars(matches("_son")), all_vars(. == "1/1")) %>%
dplyr::filter_at(vars(matches("_dat")), all_vars(. == "1/1")) %>%
dplyr::select(contains("grn")) %>%
tidyr::pivot_longer(everything()) %>%
dplyr::rename(sample = name, gt = value) %>%
dplyr::count(sample, gt, .drop = FALSE) %>%
dplyr::filter(gt %in% c("0/0", "1/1", "0/1")) %>%
mutate(n = as.numeric(n)) %>%
group_by(sample) %>%
mutate(total = as.numeric(sum(n))) %>%
dplyr::rename(obs = n) %>%
mutate(sex = case_when(
grepl("son", sample) ~ "male",
grepl("dat", sample) ~ "female"))
}
# bind all families together, to a final data frame containing all observed counts
observed <- do.call("rbind", obs)
# state the expected site proportion of each genotype and sex under a perfect mendelian segregation of F2
fem_exp <- data.frame(sex = "female", gt = c("0/0", "0/1", "1/1") , prop = c(0,0,1))
male_exp <- data.frame(sex = "male",gt = c("0/0", "0/1", "1/1"),prop = c(0,0,1))
mendel <- bind_rows(fem_exp,male_exp) %>% mutate(sex = as.character(sex))%>% mutate(gt = as.character(gt))
# mutate(type = "exp")
# using on the total genotypes count of sites per sample, calculate the expected counts, and add them to the observed counts into a final dat table
dat <- observed %>%
#mutate(sample = as.character(sample)) %>%
mutate(gt = as.character(gt)) %>%
left_join(mendel, by = c("gt", "sex")) %>%
mutate(exp = total*prop)
x <- dat %>%
dplyr::select(c("sample","gt","sex", "exp","prop", "total")) %>%
mutate(type = "exp") %>%
dplyr::rename(count = exp)
y <- dat %>%
dplyr::select(c("sample","gt","sex", "obs","total")) %>%
mutate(type = "obs") %>%
dplyr::rename(count = obs) %>%
mutate(prop = (count/total))
df_male <- rbind(x, y) %>%
dplyr::filter(sex == "male")
df_fem <- rbind(x, y) %>%
dplyr::filter(sex == "female")
p_male_pooled <- ggplot(df_male, aes(fill=gt, y=prop, x=type)) +
geom_bar(position="fill", stat="identity", ) +
xlab("F2 males") +
ylab("Proportion of F2 genotype") +
labs(title = "F2 Males of F1 male (1/1) x female (1/1)",
subtitle = paste0("Pooled sites = ", df_male %>% dplyr::filter(gt == "0/0") %>% dplyr::filter(type == "obs") %>% summarise(sum(total)) %>% as.numeric())) +
scale_x_discrete("type", labels = c("Expected","Observed")) +
labs(fill = "Genotype") +
theme_classic() +
theme(axis.text.x = element_text(size=15),
axis.title.x = element_blank())
p_female_pooled <- ggplot(df_fem, aes(fill=gt, y=prop, x=type)) +
geom_bar(position="fill", stat="identity", ) +
xlab("F2 females") +
ylab("Proportion of F2 genotype") +
labs(title = "F2 females of F1 male (1/1) x female (1/1)",
subtitle = paste0("Pooled sites = ", df_fem %>% dplyr::filter(gt == "0/0") %>% dplyr::filter(type == "obs") %>% summarise(sum(total)) %>% as.numeric())) +
scale_x_discrete("type", labels = c("Expected","Observed")) +
labs(fill = "Genotype") +
theme_classic() +
theme(axis.text.x = element_text(size=15),
axis.title.x = element_blank())
grid.arrange(p_female_pooled, p_male_pooled, ncol = 2)
ggplot(df_male, aes(fill=gt, y=prop, x=type)) +
geom_bar(position="fill", stat="identity", ) +
xlab("F2 siblings") +
ylab("Proportion of F2 genotype") +
labs(title = "F2 Males of F1 male (1/1) x female (1/1)") +
# geom_text(aes(label = n), position = position_stack(vjust = 0.5)) +
geom_text(data = filter(df_male, gt == "1/1"), aes(y=prop, x=type, label=total), vjust = 0) +
facet_wrap(~ sample) +
labs(fill = "Genotype") +
theme_classic() +
theme(axis.text.x = element_text(size=15),
axis.title.x = element_blank())
ggplot(df_fem, aes(fill=gt, y=prop, x=type)) +
geom_bar(position="fill", stat="identity", ) +
xlab("F2 siblings") +
ylab("Proportion of F2 genotype") +
labs(title = "F2 Females of F1 male (1/1) x female (1/1)") +
# geom_text(aes(label = n), position = position_stack(vjust = 0.5)) +
geom_text(data = filter(df_fem, gt == "1/1"), aes(y=prop, x=type, label=total), vjust = 0) +
facet_wrap(~ sample) +
labs(fill = "Genotype") +
theme_classic() +
theme(axis.text.x = element_text(size=15),
axis.title.x = element_blank())
we expect the F0 foundress to be homozygotic (1/1), like her son
# plot the pooled data for all F0 foundress with male (1/1) x female (1/1) offspring
ggplot(observed, aes(fill=gt, y=prop, x = type)) +
geom_bar(position="fill", stat="identity") +
ylab("Proportion of F0 genotype") +
labs(title = "F0 foundress of F1 male (1/1) x female (1/1)",
subtitle = paste0("Pooled sites = ", observed %>% dplyr::filter(gt == "0/1") %>% summarise(sum(total)) %>% as.numeric())) +
theme_classic() +
theme(axis.text.x = element_blank(), axis.title.x = element_blank()) +
labs(fill = "Genotype")
# plot for each individual:
ggplot(observed, aes(fill=gt, y=prop, x=type)) +
geom_bar(position="fill", stat="identity", ) +
xlab("Observed") +
ylab("Proportion of F0 genotype") +
labs(title = "F0 foundress of F1 male (1/1) x female (1/1)") +
geom_text(data = filter(observed, gt == "1/1"), aes(y=prop, x=type, label=total), vjust = 0) +
facet_wrap(~ sample) +
labs(fill = "Genotype") +
theme_classic() +
theme(axis.text.x = element_blank(), axis.title.x = element_blank())
We expect zero sites for this cross, since there is no paternal
inheritance to the males.
Indeed, there are only 103 sites in the F2 pooled females, and 71 for
the F2 pooled males:
# set the families (include only families with at least one adult F2)
family = grep("grndat|grnson",gt$ID, value=TRUE) %>%
str_extract("[^_]+") %>%
unique()
# define a list to put all the data frames in
obs <- list()
for (fam in family) {
obs[[fam]] <- table %>%
dplyr::select(starts_with(fam)) %>%
dplyr::filter_at(vars(matches("_son")), all_vars(. == "0/0")) %>%
dplyr::filter_at(vars(matches("_dat")), all_vars(. == "1/1")) %>%
dplyr::select(contains("grn")) %>%
tidyr::pivot_longer(everything()) %>%
dplyr::rename(sample = name, gt = value) %>%
dplyr::count(sample, gt, .drop = FALSE) %>%
dplyr::filter(gt %in% c("0/0", "1/1", "0/1")) %>%
mutate(n = as.numeric(n)) %>%
group_by(sample) %>%
mutate(total = as.numeric(sum(n))) %>%
dplyr::rename(obs = n) %>%
mutate(sex = case_when(
grepl("son", sample) ~ "male",
grepl("dat", sample) ~ "female"))
}
# bind all families together, to a final data frame containing all observed counts
observed <- do.call("rbind", obs)
# state the expected site proportion of each genotype and sex under a perfect mendelian segregation of F2
fem_exp <- data.frame(sex = "female", gt = c("0/0", "0/1", "1/1") , prop = c(0,0,0))
male_exp <- data.frame(sex = "male",gt = c("0/0", "0/1", "1/1"),prop = c(0,0,0))
mendel <- bind_rows(fem_exp,male_exp) %>% mutate(sex = as.character(sex))%>% mutate(gt = as.character(gt))
# mutate(type = "exp")
# using on the total genotypes count of sites per sample, calculate the expected counts, and add them to the observed counts into a final dat table
dat <- observed %>%
#mutate(sample = as.character(sample)) %>%
mutate(gt = as.character(gt)) %>%
left_join(mendel, by = c("gt", "sex")) %>%
mutate(exp = total*prop)
x <- dat %>%
dplyr::select(c("sample","gt","sex", "exp","prop", "total")) %>%
mutate(type = "exp") %>%
dplyr::rename(count = exp)
y <- dat %>%
dplyr::select(c("sample","gt","sex", "obs","total")) %>%
mutate(type = "obs") %>%
dplyr::rename(count = obs) %>%
mutate(prop = (count/total))
df_male <- rbind(x, y) %>%
dplyr::filter(sex == "male")
df_fem <- rbind(x, y) %>%
dplyr::filter(sex == "female")
p_male_pooled <- ggplot(df_male, aes(fill=gt, y=prop, x=type)) +
geom_bar(position="fill", stat="identity", ) +
xlab("F2 males") +
ylab("Proportion of F2 genotype") +
labs(title = "F2 Males of F1 male (0/0) x female (1/1)",
subtitle = paste0("Pooled sites = ", df_male %>% dplyr::filter(gt == "0/0") %>% dplyr::filter(type == "obs") %>% summarise(sum(total)) %>% as.numeric())) +
scale_x_discrete("type", labels = c("Expected","Observed")) +
labs(fill = "Genotype") +
theme_classic() +
theme(axis.text.x = element_text(size=15),
axis.title.x = element_blank())
p_female_pooled <- ggplot(df_fem, aes(fill=gt, y=prop, x=type)) +
geom_bar(position="fill", stat="identity", ) +
xlab("F2 females") +
ylab("Proportion of F2 genotype") +
labs(title = "F2 females of F1 male (0/0) x female (1/1)",
subtitle = paste0("Pooled sites = ", df_fem %>% dplyr::filter(gt == "0/0") %>% dplyr::filter(type == "obs") %>% summarise(sum(total)) %>% as.numeric())) +
scale_x_discrete("type", labels = c("Expected","Observed")) +
labs(fill = "Genotype") +
theme_classic() +
theme(axis.text.x = element_text(size=15),
axis.title.x = element_blank())
grid.arrange(p_female_pooled, p_male_pooled, ncol = 2)
ggplot(df_male, aes(fill=gt, y=prop, x=type)) +
geom_bar(position="fill", stat="identity", ) +
xlab("F2 siblings") +
ylab("Proportion of F2 genotype") +
labs(title = "F2 Males of F1 male (0/0) x female (1/1)") +
# geom_text(aes(label = n), position = position_stack(vjust = 0.5)) +
geom_text(data = filter(df_male, gt == "1/1"), aes(y=prop, x=type, label=total), vjust = 0) +
facet_wrap(~ sample) +
labs(fill = "Genotype") +
theme_classic() +
theme(axis.text.x = element_text(size=15),
axis.title.x = element_blank())
ggplot(df_fem, aes(fill=gt, y=prop, x=type)) +
geom_bar(position="fill", stat="identity", ) +
xlab("F2 siblings") +
ylab("Proportion of F2 genotype") +
labs(title = "F2 Females of F1 male (0/0) x female (1/1)") +
# geom_text(aes(label = n), position = position_stack(vjust = 0.5)) +
geom_text(data = filter(df_fem, gt == "1/1"), aes(y=prop, x=type, label=total), vjust = 0) +
facet_wrap(~ sample) +
labs(fill = "Genotype") +
theme_classic() +
theme(axis.text.x = element_text(size=15),
axis.title.x = element_blank())
we expect ?????…
# plot the pooled data for all F0 foundress with male (0/0) x female (1/1) offspring
ggplot(observed, aes(fill=gt, y=prop, x = type)) +
geom_bar(position="fill", stat="identity") +
ylab("Proportion of F0 genotype") +
labs(title = "F0 foundress of F1 male (0/0) x female (1/1)",
subtitle = paste0("Pooled sites = ", observed %>% dplyr::filter(gt == "0/1") %>% summarise(sum(total)) %>% as.numeric())) +
theme_classic() +
theme(axis.text.x = element_blank(), axis.title.x = element_blank()) +
labs(fill = "Genotype")
# plot for each individual:
ggplot(observed, aes(fill=gt, y=prop, x=type)) +
geom_bar(position="fill", stat="identity", ) +
xlab("Observed") +
ylab("Proportion of F0 genotype") +
labs(title = "F0 foundress of F1 male (0/0) x female (1/1)") +
geom_text(data = filter(observed, gt == "1/1"), aes(y=prop, x=type, label=total), vjust = 0) +
facet_wrap(~ sample) +
labs(fill = "Genotype") +
theme_classic() +
theme(axis.text.x = element_blank(), axis.title.x = element_blank())
We expect zero sites for this cross, since there is no paternal
inheritance to the males.
Indeed, there are only 98 sites in the F2 pooled females, and 75 for the
F2 pooled males:
# set the families (include only families with at least one adult F2)
family = grep("grndat|grnson",gt$ID, value=TRUE) %>%
str_extract("[^_]+") %>%
unique()
# define a list to put all the data frames in
obs <- list()
for (fam in family) {
obs[[fam]] <- table %>%
dplyr::select(starts_with(fam)) %>%
dplyr::filter_at(vars(matches("_son")), all_vars(. == "1/1")) %>%
dplyr::filter_at(vars(matches("_dat")), all_vars(. == "0/0")) %>%
dplyr::select(contains("grn")) %>%
tidyr::pivot_longer(everything()) %>%
dplyr::rename(sample = name, gt = value) %>%
dplyr::count(sample, gt, .drop = FALSE) %>%
dplyr::filter(gt %in% c("0/0", "1/1", "0/1")) %>%
mutate(n = as.numeric(n)) %>%
group_by(sample) %>%
mutate(total = as.numeric(sum(n))) %>%
dplyr::rename(obs = n) %>%
mutate(sex = case_when(
grepl("son", sample) ~ "male",
grepl("dat", sample) ~ "female"))
}
# bind all families together, to a final data frame containing all observed counts
observed <- do.call("rbind", obs)
# state the expected site proportion of each genotype and sex under a perfect mendelian segregation of F2
fem_exp <- data.frame(sex = "female", gt = c("0/0", "0/1", "1/1") , prop = c(0,0,0))
male_exp <- data.frame(sex = "male",gt = c("0/0", "0/1", "1/1"),prop = c(0,0,0))
mendel <- bind_rows(fem_exp,male_exp) %>% mutate(sex = as.character(sex))%>% mutate(gt = as.character(gt))
# mutate(type = "exp")
# using on the total genotypes count of sites per sample, calculate the expected counts, and add them to the observed counts into a final dat table
dat <- observed %>%
#mutate(sample = as.character(sample)) %>%
mutate(gt = as.character(gt)) %>%
left_join(mendel, by = c("gt", "sex")) %>%
mutate(exp = total*prop)
x <- dat %>%
dplyr::select(c("sample","gt","sex", "exp","prop", "total")) %>%
mutate(type = "exp") %>%
dplyr::rename(count = exp)
y <- dat %>%
dplyr::select(c("sample","gt","sex", "obs","total")) %>%
mutate(type = "obs") %>%
dplyr::rename(count = obs) %>%
mutate(prop = (count/total))
df_male <- rbind(x, y) %>%
dplyr::filter(sex == "male")
df_fem <- rbind(x, y) %>%
dplyr::filter(sex == "female")
p_male_pooled <- ggplot(df_male, aes(fill=gt, y=prop, x=type)) +
geom_bar(position="fill", stat="identity", ) +
xlab("F2 males") +
ylab("Proportion of F2 genotype") +
labs(title = "F2 Males of F1 male (1/1) x female (0/0)",
subtitle = paste0("Pooled sites = ", df_male %>% dplyr::filter(gt == "0/0") %>% dplyr::filter(type == "obs") %>% summarise(sum(total)) %>% as.numeric())) +
scale_x_discrete("type", labels = c("Expected","Observed")) +
labs(fill = "Genotype") +
theme_classic() +
theme(axis.text.x = element_text(size=15),
axis.title.x = element_blank())
p_female_pooled <- ggplot(df_fem, aes(fill=gt, y=prop, x=type)) +
geom_bar(position="fill", stat="identity", ) +
xlab("F2 females") +
ylab("Proportion of F2 genotype") +
labs(title = "F2 females of F1 male (1/1) x female (0/0)",
subtitle = paste0("Pooled sites = ", df_fem %>% dplyr::filter(gt == "0/0") %>% dplyr::filter(type == "obs") %>% summarise(sum(total)) %>% as.numeric()))+
scale_x_discrete("type", labels = c("Expected","Observed")) +
labs(fill = "Genotype") +
theme_classic() +
theme(axis.text.x = element_text(size=15),
axis.title.x = element_blank())
grid.arrange(p_female_pooled, p_male_pooled, ncol = 2)
ggplot(df_male, aes(fill=gt, y=prop, x=type)) +
geom_bar(position="fill", stat="identity", ) +
xlab("F2 siblings") +
ylab("Proportion of F2 genotype") +
labs(title = "F2 Males of F1 male (1/1) x female (0/0)") +
# geom_text(aes(label = n), position = position_stack(vjust = 0.5)) +
geom_text(data = filter(df_male, gt == "1/1"), aes(y=prop, x=type, label=total), vjust = 0) +
facet_wrap(~ sample) +
labs(fill = "Genotype") +
theme_classic() +
theme(axis.text.x = element_text(size=15),
axis.title.x = element_blank())
ggplot(df_fem, aes(fill=gt, y=prop, x=type)) +
geom_bar(position="fill", stat="identity", ) +
xlab("F2 siblings") +
ylab("Proportion of F2 genotype") +
labs(title = "F2 Females of F1 male (1/1) x female (0/0)") +
# geom_text(aes(label = n), position = position_stack(vjust = 0.5)) +
geom_text(data = filter(df_fem, gt == "1/1"), aes(y=prop, x=type, label=total), vjust = 0) +
facet_wrap(~ sample) +
labs(fill = "Genotype") +
theme_classic() +
theme(axis.text.x = element_text(size=15),
axis.title.x = element_blank())
we expect ?????…
# plot the pooled data for all F0 foundress with male (1/1) x female (0/0) offspring
ggplot(observed, aes(fill=gt, y=prop, x = type)) +
geom_bar(position="fill", stat="identity") +
ylab("Proportion of F0 genotype") +
labs(title = "F0 foundress of F1 male (1/1) x female (0/0)",
subtitle = paste0("Pooled sites = ", observed %>% dplyr::filter(gt == "0/1") %>% summarise(sum(total)) %>% as.numeric())) +
theme_classic() +
theme(axis.text.x = element_blank(), axis.title.x = element_blank()) +
labs(fill = "Genotype")
# plot for each individual:
ggplot(observed, aes(fill=gt, y=prop, x=type)) +
geom_bar(position="fill", stat="identity", ) +
xlab("Observed") +
ylab("Proportion of F0 genotype") +
labs(title = "F0 foundress of F1 male (1/1) x female (0/0)") +
geom_text(data = filter(observed, gt == "1/1"), aes(y=prop, x=type, label=total), vjust = 0) +
facet_wrap(~ sample) +
labs(fill = "Genotype") +
theme_classic() +
theme(axis.text.x = element_blank(), axis.title.x = element_blank())
Based the ‘Control crossing’ results, I made a vector with all ‘False sites’. These sites will be excluded from the following analyses (5-9).
# based on the former control crossing - make a vector of the false sites:
# deduct these sites form the 'table' file
We expect F2 females’ sites to be evenly distributed between 0/1 and 0/0 genotypes, and F2 males to be evenly distributed between 0/0 and 1/1 (as they are hemizygote).
# set the families (include only families with at least one adult F2)
family = grep("grndat|grnson",gt$ID, value=TRUE) %>%
str_extract("[^_]+") %>%
unique()
# define a list to put all the data frames in
obs <- list()
for (fam in family) {
obs[[fam]] <- table %>%
dplyr::select(starts_with(fam)) %>%
dplyr::filter_at(vars(matches("_son")), all_vars(. == "0/0")) %>%
dplyr::filter_at(vars(matches("_dat")), all_vars(. == "0/1")) %>%
dplyr::select(contains("grn")) %>%
tidyr::pivot_longer(everything()) %>%
dplyr::rename(sample = name, gt = value) %>%
dplyr::count(sample, gt, .drop = FALSE) %>%
dplyr::filter(gt %in% c("0/0", "1/1", "0/1")) %>%
mutate(n = as.numeric(n)) %>%
group_by(sample) %>%
mutate(total = as.numeric(sum(n))) %>%
dplyr::rename(obs = n) %>%
mutate(sex = case_when(
grepl("son", sample) ~ "male",
grepl("dat", sample) ~ "female"))
}
# bind all families together, to a final data frame containing all observed counts
observed <- do.call("rbind", obs)
# state the expected site proportion of each genotype and sex under a perfect mendelian segregation of F2
fem_exp <- data.frame(sex = "female", gt = c("0/0", "0/1", "1/1") , prop = c(0.5, 0.5,0))
male_exp <- data.frame(sex = "male",gt = c("0/0", "0/1", "1/1"),prop = c(0.5, 0,0.5))
mendel <- bind_rows(fem_exp,male_exp) %>% mutate(sex = as.character(sex))%>% mutate(gt = as.character(gt))
# mutate(type = "exp")
# using on the total genotypes count of sites per sample, calculate the expected counts, and add them to the observed counts into a final dat table
dat <- observed %>%
#mutate(sample = as.character(sample)) %>%
mutate(gt = as.character(gt)) %>%
left_join(mendel, by = c("gt", "sex")) %>%
mutate(exp = total*prop)
x <- dat %>%
dplyr::select(c("sample","gt","sex", "exp","prop", "total")) %>%
mutate(type = "exp") %>%
dplyr::rename(count = exp)
y <- dat %>%
dplyr::select(c("sample","gt","sex", "obs","total")) %>%
mutate(type = "obs") %>%
dplyr::rename(count = obs) %>%
mutate(prop = (count/total))
df_male <- rbind(x, y) %>%
dplyr::filter(sex == "male")
df_fem <- rbind(x, y) %>%
dplyr::filter(sex == "female")
p_male_pooled <- ggplot(df_male, aes(fill=gt, y=prop, x=type)) +
geom_bar(position="fill", stat="identity", ) +
xlab("F2 males") +
ylab("Proportion of F2 genotype") +
labs(title = "F2 Males of F1 male (0/0) x female (0/1)",
subtitle = paste0("Pooled sites = ", df_male %>% dplyr::filter(gt == "0/0") %>% dplyr::filter(type == "obs") %>% summarise(sum(total)) %>% as.numeric())) +
scale_x_discrete("type", labels = c("Expected","Observed")) +
labs(fill = "Genotype") +
theme_classic() +
theme(axis.text.x = element_text(size=15),
axis.title.x = element_blank())
p_female_pooled <- ggplot(df_fem, aes(fill=gt, y=prop, x=type)) +
geom_bar(position="fill", stat="identity", ) +
xlab("F2 females") +
ylab("Proportion of F2 genotype") +
labs(title = "F2 females of F1 male (0/0) x female (0/1)",
subtitle = paste0("Pooled sites = ", df_fem %>% dplyr::filter(gt == "0/0") %>% dplyr::filter(type == "obs") %>% summarise(sum(total)) %>% as.numeric())) +
scale_x_discrete("type", labels = c("Expected","Observed")) +
labs(fill = "Genotype") +
theme_classic() +
theme(axis.text.x = element_text(size=15),
axis.title.x = element_blank())
grid.arrange(p_female_pooled, p_male_pooled, ncol = 2)
ggplot(df_male, aes(fill=gt, y=prop, x=type)) +
geom_bar(position="fill", stat="identity", ) +
xlab("F2 siblings") +
ylab("Proportion of F2 genotype") +
labs(title = "F2 Males of F1 male (0/0) x female (0/1)") +
# geom_text(aes(label = n), position = position_stack(vjust = 0.5)) +
geom_text(data = filter(df_male, gt == "1/1"), aes(y=prop, x=type, label=total), vjust = 0) +
facet_wrap(~ sample) +
labs(fill = "Genotype") +
theme_classic()
ggplot(df_fem, aes(fill=gt, y=prop, x=type)) +
geom_bar(position="fill", stat="identity", ) +
xlab("F2 siblings") +
ylab("Proportion of F2 genotype") +
labs(title = "F2 Females of F1 male (0/0) x female (0/1)") +
# geom_text(aes(label = n), position = position_stack(vjust = 0.5)) +
geom_text(data = filter(df_fem, gt == "1/1"), aes(y=prop, x=type, label=total), vjust = 0) +
facet_wrap(~ sample) +
labs(fill = "Genotype") +
theme_classic()
we expect ?????…
# plot the pooled data for all F0 foundress with male (0/0) x female (0/1) offspring
ggplot(observed, aes(fill=gt, y=prop, x = type)) +
geom_bar(position="fill", stat="identity") +
ylab("Proportion of F0 genotype") +
labs(title = "F0 foundress of F1 male (0/0) x female (0/1)",
subtitle = paste0("Pooled sites = ", observed %>% dplyr::filter(gt == "0/1") %>% summarise(sum(total)) %>% as.numeric())) +
theme_classic() +
theme(axis.text.x = element_blank(), axis.title.x = element_blank()) +
labs(fill = "Genotype")
# plot for each individual:
ggplot(observed, aes(fill=gt, y=prop, x=type)) +
geom_bar(position="fill", stat="identity", ) +
xlab("Observed") +
ylab("Proportion of F0 genotype") +
labs(title = "F0 foundress of F1 male (0/0) x female (0/1)") +
geom_text(data = filter(observed, gt == "1/1"), aes(y=prop, x=type, label=total), vjust = 0) +
facet_wrap(~ sample) +
labs(fill = "Genotype") +
theme_classic() +
theme(axis.text.x = element_blank(), axis.title.x = element_blank())
Based the ‘Control crossing’ results, I made a vector with all ‘False sites’. These sites will be excluded from the following analyses (5-9).
We expect F2 females’ sites to be evenly distributed between 0/1 and 1/1 genotypes, and F2 males to be evenly distributed between 0/0 and 1/1 (as they are hemizygote).
# set the families (include only families with at least one adult F2)
family = grep("grndat|grnson",gt$ID, value=TRUE) %>%
str_extract("[^_]+") %>%
unique()
# define a list to put all the data frames in
obs <- list()
for (fam in family) {
obs[[fam]] <- table %>%
dplyr::select(starts_with(fam)) %>%
dplyr::filter_at(vars(matches("_son")), all_vars(. == "1/1")) %>%
dplyr::filter_at(vars(matches("_dat")), all_vars(. == "0/1")) %>%
dplyr::select(contains("grn")) %>%
tidyr::pivot_longer(everything()) %>%
dplyr::rename(sample = name, gt = value) %>%
dplyr::count(sample, gt, .drop = FALSE) %>%
dplyr::filter(gt %in% c("0/0", "1/1", "0/1")) %>%
mutate(n = as.numeric(n)) %>%
group_by(sample) %>%
mutate(total = as.numeric(sum(n))) %>%
dplyr::rename(obs = n) %>%
mutate(sex = case_when(
grepl("son", sample) ~ "male",
grepl("dat", sample) ~ "female"))
}
# bind all families together, to a final data frame containing all observed counts
observed <- do.call("rbind", obs)
# state the expected site proportion of each genotype and sex under a perfect mendelian segregation of F2
fem_exp <- data.frame(sex = "female", gt = c("0/0", "0/1", "1/1") , prop = c(0,0.5,1))
male_exp <- data.frame(sex = "male",gt = c("0/0", "0/1", "1/1"),prop = c(0.5,0,0.5))
mendel <- bind_rows(fem_exp,male_exp) %>% mutate(sex = as.character(sex))%>% mutate(gt = as.character(gt))
# mutate(type = "exp")
# using on the total genotypes count of sites per sample, calculate the expected counts, and add them to the observed counts into a final dat table
dat <- observed %>%
#mutate(sample = as.character(sample)) %>%
mutate(gt = as.character(gt)) %>%
left_join(mendel, by = c("gt", "sex")) %>%
mutate(exp = total*prop)
x <- dat %>%
dplyr::select(c("sample","gt","sex", "exp","prop", "total")) %>%
mutate(type = "exp") %>%
dplyr::rename(count = exp)
y <- dat %>%
dplyr::select(c("sample","gt","sex", "obs","total")) %>%
mutate(type = "obs") %>%
dplyr::rename(count = obs) %>%
mutate(prop = (count/total))
df_male <- rbind(x, y) %>%
dplyr::filter(sex == "male")
df_fem <- rbind(x, y) %>%
dplyr::filter(sex == "female")
p_male_pooled <- ggplot(df_male, aes(fill=gt, y=prop, x=type)) +
geom_bar(position="fill", stat="identity", ) +
xlab("F2 males") +
ylab("Proportion of F2 genotype") +
labs(title = "F2 Males of F1 male (1/1) x female (0/1)",
subtitle = paste0("Pooled sites = ", df_male %>% dplyr::filter(gt == "0/0") %>% dplyr::filter(type == "obs") %>% summarise(sum(total)) %>% as.numeric())) +
scale_x_discrete("type", labels = c("Expected","Observed")) +
labs(fill = "Genotype") +
theme_classic() +
theme(axis.text.x = element_text(size=15),
axis.title.x = element_blank())
p_female_pooled <- ggplot(df_fem, aes(fill=gt, y=prop, x=type)) +
geom_bar(position="fill", stat="identity", ) +
xlab("F2 females") +
ylab("Proportion of F2 genotype") +
labs(title = "F2 females of F1 male (1/1) x female (0/1)",
subtitle = paste0("Pooled sites = ", df_fem %>% dplyr::filter(gt == "0/0") %>% dplyr::filter(type == "obs") %>% summarise(sum(total)) %>% as.numeric())) +
scale_x_discrete("type", labels = c("Expected","Observed")) +
labs(fill = "Genotype") +
theme_classic() +
theme(axis.text.x = element_text(size=15),
axis.title.x = element_blank())
grid.arrange(p_female_pooled, p_male_pooled, ncol = 2)
ggplot(df_male, aes(fill=gt, y=prop, x=type)) +
geom_bar(position="fill", stat="identity", ) +
xlab("F2 siblings") +
ylab("Proportion of F2 genotype") +
labs(title = "F2 Males of F1 male (1/1) x female (0/1)") +
# geom_text(aes(label = n), position = position_stack(vjust = 0.5)) +
geom_text(data = filter(df_male, gt == "1/1"), aes(y=prop, x=type, label=total), vjust = 0) +
facet_wrap(~ sample) +
labs(fill = "Genotype") +
theme_classic()
ggplot(df_fem, aes(fill=gt, y=prop, x=type)) +
geom_bar(position="fill", stat="identity", ) +
xlab("F2 siblings") +
ylab("Proportion of F2 genotype") +
labs(title = "F2 Females of F1 male (1/1) x female (0/1)") +
# geom_text(aes(label = n), position = position_stack(vjust = 0.5)) +
geom_text(data = filter(df_fem, gt == "1/1"), aes(y=prop, x=type, label=total), vjust = 0) +
facet_wrap(~ sample) +
labs(fill = "Genotype") +
theme_classic()
we expect ?????…
# plot the pooled data for all F0 foundress with male (1/1) x female (0/1) offspring
ggplot(observed, aes(fill=gt, y=prop, x = type)) +
geom_bar(position="fill", stat="identity") +
ylab("Proportion of F0 genotype") +
labs(title = "F0 foundress of F1 male (1/1) x female (0/1)",
subtitle = paste0("Pooled sites = ", observed %>% dplyr::filter(gt == "0/1") %>% summarise(sum(total)) %>% as.numeric())) +
theme_classic() +
theme(axis.text.x = element_blank(), axis.title.x = element_blank()) +
labs(fill = "Genotype")
# plot for each individual:
ggplot(observed, aes(fill=gt, y=prop, x=type)) +
geom_bar(position="fill", stat="identity", ) +
xlab("Observed") +
ylab("Proportion of F0 genotype") +
labs(title = "F0 foundress of F1 male (1/1) x female (0/1)") +
geom_text(data = filter(observed, gt == "1/1"), aes(y=prop, x=type, label=total), vjust = 0) +
facet_wrap(~ sample) +
labs(fill = "Genotype") +
theme_classic() +
theme(axis.text.x = element_blank(), axis.title.x = element_blank())
Based the ‘Control crossing’ results, I made a vector with all ‘False sites’. These sites will be excluded from the following analyses (5-9).
The former crosses of heterozygotic females (5 and 6) show that F2
males can be heterozygotic and carry two alleles, like their
mother.
But are these sites real? and if they are, can they be
transmitted to their daughters?
Looking at the F2 offspring of crossing heterozygotic F1 males with
females, we tested 3 possible predictions:
# set the families (include only families with at least one adult F2)
family = grep("grndat|grnson",gt$ID, value=TRUE) %>%
str_extract("[^_]+") %>%
unique()
# define a list to put all the data frames in
obs <- list()
for (fam in family) {
obs[[fam]] <- table %>%
dplyr::select(starts_with(fam)) %>%
dplyr::filter_at(vars(matches("_son")), all_vars(. == "0/1")) %>%
dplyr::filter_at(vars(matches("_dat")), all_vars(. == "0/0")) %>%
dplyr::select(contains("grn")) %>%
tidyr::pivot_longer(everything()) %>%
dplyr::rename(sample = name, gt = value) %>%
dplyr::count(sample, gt, .drop = FALSE) %>%
dplyr::filter(gt %in% c("0/0", "1/1", "0/1")) %>%
mutate(n = as.numeric(n)) %>%
group_by(sample) %>%
mutate(total = as.numeric(sum(n))) %>%
dplyr::rename(obs = n) %>%
mutate(sex = case_when(
grepl("son", sample) ~ "male",
grepl("dat", sample) ~ "female"))
}
# bind all families together, to a final data frame containing all observed counts
observed <- do.call("rbind", obs)
# state the expected site proportion of each genotype and sex under the different hypothesis:
# Heterozygotic sites are real:Males can carry and inherit two alleles
fem_sitesReal <- data.frame(type = "sitesReal", sex = "female", gt = c("0/0", "0/1", "1/1") , prop = c(0.5, 0.5, 0))
# Heterozygotic sites are false:Males can carry and inherit only one allele
# the heterozygotic sites are actually (0)
fem_sitesFalse_0 <- data.frame(type = "sitesFalse_0",sex = "female", gt = c("0/0", "0/1", "1/1") , prop = c(1,0,0))
# the heterozygotic sites are actually (1)
fem_sitesFalse_1 <- data.frame(type = "sitesFalse_1", sex = "female", gt = c("0/0", "0/1", "1/1") , prop = c(0,1,0))
# males F2 genotype is expected to stay the same , as it is produced from unfertilized egg, hence have no paternal inheritance
male_sitesReal <- data.frame(type = "sitesReal",sex = "male",gt = c("0/0", "0/1", "1/1"),prop = c(1,0,0))
male_sitesFalse_0 <- data.frame(type = "sitesFalse_0",sex = "male",gt = c("0/0", "0/1", "1/1"),prop = c(1,0,0))
male_sitesFalse_1 <- data.frame(type = "sitesFalse_1", sex = "male",gt = c("0/0", "0/1", "1/1"),prop = c(1,0,0))
expected_prop <- bind_rows(fem_sitesReal,fem_sitesFalse_0,fem_sitesFalse_1,male_sitesReal,male_sitesFalse_0,male_sitesFalse_1) %>% mutate(sex = as.character(sex))%>% mutate(gt = as.character(gt))
# mutate(type = "exp")
# using on the total genotypes count of sites per sample, calculate the expected counts, and add them to the observed counts into a final dat table
expected <- observed %>%
mutate(gt = as.character(gt)) %>%
full_join(expected_prop, by = c("gt", "sex")) %>%
mutate(count = total*prop)
x <- observed %>%
dplyr::select(c("sample","sex","gt","total", "obs")) %>%
mutate(type = "Observed") %>%
dplyr::rename(count = obs) %>%
mutate(prop = (count/total))
dat <- bind_rows(x, expected) %>%
dplyr::select(-obs)
# subset for males and females, and reorder the "types" levels:
dat_fem <- dat %>% dplyr::filter(sex == "female")
level_order_fem <- factor(dat_fem$type, level = c("sitesFalse_0", "sitesFalse_1", "sitesReal", "Observed"))
dat_male <- dat %>% dplyr::filter(sex == "male")
level_order_male <- factor(dat_male$type, level = c("sitesFalse_0", "sitesFalse_1", "sitesReal", "Observed"))
p_male_pooled <- ggplot(dat_male, aes(fill=gt, y=prop, x=level_order_male)) +
geom_bar(position="fill", stat="identity", width = 0.8) +
ylab("Proportion of F2 genotype") +
labs(title = paste0("F2 Male, pooled sites = ", dat_male %>% dplyr::filter(type == "Observed") %>% summarise(sum(count)) %>% as.numeric())) +
labs(fill = "Genotype:") +
theme_classic() +
theme(axis.text.x = element_text(size=12, angle = 45,hjust = 1),
axis.title.x = element_blank())
p_female_pooled <- ggplot(dat_fem, aes(fill=gt, y=prop, x=level_order_fem)) +
geom_bar(position="fill", stat="identity", width = 0.8) +
ylab("Proportion of F2 genotype") +
labs(title = paste0("F2 Female, pooled sites = ", dat_fem %>% dplyr::filter(type == "Observed") %>% summarise(sum(count)) %>% as.numeric())) +
labs(fill = "Genotype:") +
theme_classic() +
theme(axis.text.x = element_text(size=12, angle = 45,hjust = 1),
axis.title.x = element_blank())
combined <- p_female_pooled + p_male_pooled & theme(legend.position = "bottom")
combined + plot_layout(guides = "collect") + plot_annotation(
title = 'F2 genotype. F1 Cross: male (0/1) x female (0/0)')
# plot each individual:
ggplot(dat_male, aes(fill=gt, y=prop, x=level_order_male)) +
geom_bar(position="fill", stat="identity", ) +
ylab("Proportion of F2 genotype") +
labs(title = "F2 male genotype. F1 Cross: male (0/1) x female (0/0)", fill = "Genotype:") +
geom_text(data = dat_male %>% filter(gt == "1/1") %>% filter(type == "Observed"), aes(x = 4, y = 0.8, label = total), inherit.aes = FALSE) +
facet_wrap(~ sample) +
theme_classic()+
theme(legend.position = "bottom", axis.title.x = element_blank(), axis.text.x = element_text(size=8, angle = 45,hjust = 1))
ggplot(dat_fem, aes(fill=gt, y=prop, x=level_order_fem)) +
geom_bar(position="fill", stat="identity", ) +
ylab("Proportion of F2 genotype") +
labs(title = "F2 female genotype. F1 Cross: male (0/1) x female (0/0)", fill = "Genotype:") +
geom_text(data = dat_fem %>% filter(gt == "1/1") %>% filter(type == "Observed"), aes(x = 4, y = 0.8, label = total), inherit.aes = FALSE) +
facet_wrap(~ sample) +
theme_classic()+
theme(legend.position = "bottom", axis.title.x = element_blank(), axis.text.x = element_text(size=8, angle = 45,hjust = 1))
we expect ?????…
# plot the pooled data for all F0 foundress with male (0/1) x female (0/0) offspring
ggplot(observed, aes(fill=gt, y=prop, x = type)) +
geom_bar(position="fill", stat="identity") +
ylab("Proportion of F0 genotype") +
labs(title = "F0 foundress of F1 male (0/1) x female (0/0)",
subtitle = paste0("Pooled sites = ", observed %>% dplyr::filter(gt == "0/1") %>% summarise(sum(total)) %>% as.numeric())) +
theme_classic() +
theme(axis.text.x = element_blank(), axis.title.x = element_blank()) +
labs(fill = "Genotype")
# plot for each individual:
ggplot(observed, aes(fill=gt, y=prop, x=type)) +
geom_bar(position="fill", stat="identity", ) +
xlab("Observed") +
ylab("Proportion of F0 genotype") +
labs(title = "F0 foundress of F1 male (0/1) x female (0/0)") +
geom_text(data = filter(observed, gt == "1/1"), aes(y=prop, x=type, label=total), vjust = 0) +
facet_wrap(~ sample) +
labs(fill = "Genotype") +
theme_classic() +
theme(axis.text.x = element_blank(), axis.title.x = element_blank())
Based the ‘Control crossing’ results, I made a vector with all ‘False sites’. These sites will be excluded from the following analyses (5-9).
# set the families (include only families with at least one adult F2)
family = grep("grndat|grnson",gt$ID, value=TRUE) %>%
str_extract("[^_]+") %>%
unique()
# define a list to put all the data frames in
obs <- list()
for (fam in family) {
obs[[fam]] <- table %>%
dplyr::select(starts_with(fam)) %>%
dplyr::filter_at(vars(matches("_son")), all_vars(. == "0/1")) %>%
dplyr::filter_at(vars(matches("_dat")), all_vars(. == "1/1")) %>%
dplyr::select(contains("grn")) %>%
tidyr::pivot_longer(everything()) %>%
dplyr::rename(sample = name, gt = value) %>%
dplyr::count(sample, gt, .drop = FALSE) %>%
dplyr::filter(gt %in% c("0/0", "1/1", "0/1")) %>%
mutate(n = as.numeric(n)) %>%
group_by(sample) %>%
mutate(total = as.numeric(sum(n))) %>%
dplyr::rename(obs = n) %>%
mutate(sex = case_when(
grepl("son", sample) ~ "male",
grepl("dat", sample) ~ "female"))
}
# bind all families together, to a final data frame containing all observed counts
observed <- do.call("rbind", obs)
# state the expected site proportion of each genotype and sex under the different hypothesis:
# Heterozygotic sites are real:Males can carry and inherit two alleles
fem_sitesReal <- data.frame(type = "sitesReal", sex = "female", gt = c("0/0", "0/1", "1/1") , prop = c(0,0.5,0.5))
# Heterozygotic sites are false:Males can carry and inherit only one allele
# the heterozygotic sites are actually (0)
fem_sitesFalse_0 <- data.frame(type = "sitesFalse_0",sex = "female", gt = c("0/0", "0/1", "1/1") , prop = c(0,1,0))
# the heterozygotic sites are actually (1)
fem_sitesFalse_1 <- data.frame(type = "sitesFalse_1", sex = "female", gt = c("0/0", "0/1", "1/1") , prop = c(0,0,1))
# males F2 genotype is expected to stay the same , as it is produced from unfertilized egg, hence have no paternal inheritance
male_sitesReal <- data.frame(type = "sitesReal",sex = "male",gt = c("0/0", "0/1", "1/1"),prop = c(0,0,1))
male_sitesFalse_0 <- data.frame(type = "sitesFalse_0",sex = "male",gt = c("0/0", "0/1", "1/1"),prop = c(0,0,1))
male_sitesFalse_1 <- data.frame(type = "sitesFalse_1", sex = "male",gt = c("0/0", "0/1", "1/1"),prop = c(0,0,1))
expected_prop <- bind_rows(fem_sitesReal,fem_sitesFalse_0,fem_sitesFalse_1,male_sitesReal,male_sitesFalse_0,male_sitesFalse_1) %>% mutate(sex = as.character(sex))%>% mutate(gt = as.character(gt))
# mutate(type = "exp")
# using on the total genotypes count of sites per sample, calculate the expected counts, and add them to the observed counts into a final dat table
expected <- observed %>%
mutate(gt = as.character(gt)) %>%
full_join(expected_prop, by = c("gt", "sex")) %>%
mutate(count = total*prop)
x <- observed %>%
dplyr::select(c("sample","sex","gt","total", "obs")) %>%
mutate(type = "Observed") %>%
dplyr::rename(count = obs) %>%
mutate(prop = (count/total))
dat <- bind_rows(x, expected) %>%
dplyr::select(-obs)
# subset for males and females, and reorder the "types" levels:
dat_fem <- dat %>% dplyr::filter(sex == "female")
level_order_fem <- factor(dat_fem$type, level = c("sitesFalse_0", "sitesFalse_1", "sitesReal", "Observed"))
dat_male <- dat %>% dplyr::filter(sex == "male")
level_order_male <- factor(dat_male$type, level = c("sitesFalse_0", "sitesFalse_1", "sitesReal", "Observed"))
p_male_pooled <- ggplot(dat_male, aes(fill=gt, y=prop, x=level_order_male)) +
geom_bar(position="fill", stat="identity", width = 0.8) +
ylab("Proportion of F2 genotype") +
labs(title = paste0("F2 Male, pooled sites = ", dat_male %>% dplyr::filter(type == "Observed") %>% summarise(sum(count)) %>% as.numeric())) +
labs(fill = "Genotype:") +
theme_classic() +
theme(axis.text.x = element_text(size=12, angle = 45,hjust = 1),
axis.title.x = element_blank())
p_female_pooled <- ggplot(dat_fem, aes(fill=gt, y=prop, x=level_order_fem)) +
geom_bar(position="fill", stat="identity", width = 0.8) +
ylab("Proportion of F2 genotype") +
labs(title = paste0("F2 Female, pooled sites = ", dat_fem %>% dplyr::filter(type == "Observed") %>% summarise(sum(count)) %>% as.numeric())) +
labs(fill = "Genotype:") +
theme_classic() +
theme(axis.text.x = element_text(size=12, angle = 45,hjust = 1),
axis.title.x = element_blank())
combined <- p_female_pooled + p_male_pooled & theme(legend.position = "bottom")
combined + plot_layout(guides = "collect") + plot_annotation(
title = 'F2 genotype. F1 Cross: male (0/1) x female (1/1)')
# plot each individual:
ggplot(dat_male, aes(fill=gt, y=prop, x=level_order_male)) +
geom_bar(position="fill", stat="identity", ) +
ylab("Proportion of F2 genotype") +
labs(title = "F2 male genotype. F1 Cross: male (0/1) x female (1/1)", fill = "Genotype:") +
geom_text(data = dat_male %>% filter(gt == "1/1") %>% filter(type == "Observed"), aes(x = 4, y = 0.8, label = total), inherit.aes = FALSE) +
facet_wrap(~ sample) +
theme_classic()+
theme(legend.position = "bottom", axis.title.x = element_blank(), axis.text.x = element_text(size=8, angle = 45,hjust = 1))
ggplot(dat_fem, aes(fill=gt, y=prop, x=level_order_fem)) +
geom_bar(position="fill", stat="identity", ) +
ylab("Proportion of F2 genotype") +
labs(title = "F2 female genotype. F1 Cross: male (0/1) x female (1/1)", fill = "Genotype:") +
geom_text(data = dat_fem %>% filter(gt == "1/1") %>% filter(type == "Observed"), aes(x = 4, y = 0.8, label = total), inherit.aes = FALSE) +
facet_wrap(~ sample) +
theme_classic()+
theme(legend.position = "bottom", axis.title.x = element_blank(), axis.text.x = element_text(size=8, angle = 45,hjust = 1))
we expect ?????…
# plot the pooled data for all F0 foundress with male (0/1) x female (1/1) offspring
ggplot(observed, aes(fill=gt, y=prop, x = type)) +
geom_bar(position="fill", stat="identity") +
ylab("Proportion of F0 genotype") +
labs(title = "F0 foundress of F1 male (0/1) x female (1/1)",
subtitle = paste0("Pooled sites = ", observed %>% dplyr::filter(gt == "0/1") %>% summarise(sum(total)) %>% as.numeric())) +
theme_classic() +
theme(axis.text.x = element_blank(), axis.title.x = element_blank()) +
labs(fill = "Genotype")
# plot for each individual:
ggplot(observed, aes(fill=gt, y=prop, x=type)) +
geom_bar(position="fill", stat="identity", ) +
xlab("Observed") +
ylab("Proportion of F0 genotype") +
labs(title = "F0 foundress of F1 male (0/1) x female (1/1)") +
geom_text(data = filter(observed, gt == "1/1"), aes(y=prop, x=type, label=total), vjust = 0) +
facet_wrap(~ sample) +
labs(fill = "Genotype") +
theme_classic() +
theme(axis.text.x = element_blank(), axis.title.x = element_blank())
Based the ‘Control crossing’ results, I made a vector with all ‘False sites’. These sites will be excluded from the following analyses (5-9).
# set the families (include only families with at least one adult F2)
family = grep("grndat|grnson",gt$ID, value=TRUE) %>%
str_extract("[^_]+") %>%
unique()
# define a list to put all the data frames in
obs <- list()
for (fam in family) {
obs[[fam]] <- table %>%
dplyr::select(starts_with(fam)) %>%
dplyr::filter_at(vars(matches("_son")), all_vars(. == "0/1")) %>%
dplyr::filter_at(vars(matches("_dat")), all_vars(. == "0/1")) %>%
dplyr::select(contains("grn")) %>%
tidyr::pivot_longer(everything()) %>%
dplyr::rename(sample = name, gt = value) %>%
dplyr::count(sample, gt, .drop = FALSE) %>%
dplyr::filter(gt %in% c("0/0", "1/1", "0/1")) %>%
mutate(n = as.numeric(n)) %>%
group_by(sample) %>%
mutate(total = as.numeric(sum(n))) %>%
dplyr::rename(obs = n) %>%
mutate(sex = case_when(
grepl("son", sample) ~ "male",
grepl("dat", sample) ~ "female"))
}
# bind all families together, to a final data frame containing all observed counts
observed <- do.call("rbind", obs)
# state the expected site proportion of each genotype and sex under the different hypothesis:
# Heterozygotic sites are real:Males can carry and inherit two alleles
fem_sitesReal <- data.frame(type = "sitesReal", sex = "female", gt = c("0/0", "0/1", "1/1") , prop = c(0.25, 0.5, 0.25))
# Heterozygotic sites are false:Males can carry and inherit only one allele
# the heterozygotic sites are actually (0)
fem_sitesFalse_0 <- data.frame(type = "sitesFalse_0",sex = "female", gt = c("0/0", "0/1", "1/1") , prop = c(0.5, 0.5, 0))
# the heterozygotic sites are actually (1)
fem_sitesFalse_1 <- data.frame(type = "sitesFalse_1", sex = "female", gt = c("0/0", "0/1", "1/1") , prop = c(0, 0.5, 0.5))
# males F2 genotype is expected to stay the same , as it is produced from unfertilized egg, hence have no paternal inheritance
male_sitesReal <- data.frame(type = "sitesReal",sex = "male",gt = c("0/0", "0/1", "1/1"),prop = c(0.5, 0, 0.5))
male_sitesFalse_0 <- data.frame(type = "sitesFalse_0",sex = "male",gt = c("0/0", "0/1", "1/1"),prop = c(0.5, 0, 0.5))
male_sitesFalse_1 <- data.frame(type = "sitesFalse_1", sex = "male",gt = c("0/0", "0/1", "1/1"),prop = c(0.5, 0, 0.5))
expected_prop <- bind_rows(fem_sitesReal,fem_sitesFalse_0,fem_sitesFalse_1,male_sitesReal,male_sitesFalse_0,male_sitesFalse_1) %>% mutate(sex = as.character(sex))%>% mutate(gt = as.character(gt))
# mutate(type = "exp")
# using on the total genotypes count of sites per sample, calculate the expected counts, and add them to the observed counts into a final dat table
expected <- observed %>%
mutate(gt = as.character(gt)) %>%
full_join(expected_prop, by = c("gt", "sex")) %>%
mutate(count = total*prop)
x <- observed %>%
dplyr::select(c("sample","sex","gt","total", "obs")) %>%
mutate(type = "Observed") %>%
dplyr::rename(count = obs) %>%
mutate(prop = (count/total))
dat <- bind_rows(x, expected) %>%
dplyr::select(-obs)
# subset for males and females, and reorder the "types" levels:
dat_fem <- dat %>% dplyr::filter(sex == "female")
level_order_fem <- factor(dat_fem$type, level = c("sitesFalse_0", "sitesFalse_1", "sitesReal", "Observed"))
dat_male <- dat %>% dplyr::filter(sex == "male")
level_order_male <- factor(dat_male$type, level = c("sitesFalse_0", "sitesFalse_1", "sitesReal", "Observed"))
p_male_pooled <- ggplot(dat_male, aes(fill=gt, y=prop, x=level_order_male)) +
geom_bar(position="fill", stat="identity", width = 0.8) +
ylab("Proportion of F2 genotype") +
labs(title = paste0("F2 Male, pooled sites = ", dat_male %>% dplyr::filter(type == "Observed") %>% summarise(sum(count)) %>% as.numeric())) +
labs(fill = "Genotype:") +
theme_classic() +
theme(axis.text.x = element_text(size=12, angle = 45,hjust = 1),
axis.title.x = element_blank())
p_female_pooled <- ggplot(dat_fem, aes(fill=gt, y=prop, x=level_order_fem)) +
geom_bar(position="fill", stat="identity", width = 0.8) +
ylab("Proportion of F2 genotype") +
labs(title = paste0("F2 Female, pooled sites = ", dat_fem %>% dplyr::filter(type == "Observed") %>% summarise(sum(count)) %>% as.numeric())) +
labs(fill = "Genotype:") +
theme_classic() +
theme(axis.text.x = element_text(size=12, angle = 45,hjust = 1),
axis.title.x = element_blank())
combined <- p_female_pooled + p_male_pooled & theme(legend.position = "bottom")
combined + plot_layout(guides = "collect") + plot_annotation(
title = 'F2 genotype. F1 Cross: male (0/1) x female (0/1)')
# plot each individual:
ggplot(dat_male, aes(fill=gt, y=prop, x=level_order_male)) +
geom_bar(position="fill", stat="identity", ) +
ylab("Proportion of F2 genotype") +
labs(title = "F2 male genotype. F1 Cross: male (0/1) x female (0/1)", fill = "Genotype:") +
geom_text(data = dat_male %>% filter(gt == "1/1") %>% filter(type == "Observed"), aes(x = 4, y = 0.8, label = total), inherit.aes = FALSE) +
facet_wrap(~ sample) +
theme_classic()+
theme(legend.position = "bottom", axis.title.x = element_blank(), axis.text.x = element_text(size=8, angle = 45,hjust = 1))
ggplot(dat_fem, aes(fill=gt, y=prop, x=level_order_fem)) +
geom_bar(position="fill", stat="identity", ) +
ylab("Proportion of F2 genotype") +
labs(title = "F2 female genotype. F1 Cross: male (0/1) x female (0/1)", fill = "Genotype:") +
geom_text(data = dat_fem %>% filter(gt == "1/1") %>% filter(type == "Observed"), aes(x = 4, y = 0.8, label = total), inherit.aes = FALSE) +
facet_wrap(~ sample) +
theme_classic()+
theme(legend.position = "bottom", axis.title.x = element_blank(), axis.text.x = element_text(size=8, angle = 45,hjust = 1))
we expect ?????…
# plot the pooled data for all F0 foundress with male (0/1) x female (1/1) offspring
ggplot(observed, aes(fill=gt, y=prop, x = type)) +
geom_bar(position="fill", stat="identity") +
ylab("Proportion of F0 genotype") +
labs(title = "F0 foundress of F1 male (0/1) x female (0/1)",
subtitle = paste0("Pooled sites = ", observed %>% dplyr::filter(gt == "1/1") %>% summarise(sum(total)) %>% as.numeric())) +
theme_classic() +
theme(axis.text.x = element_blank(), axis.title.x = element_blank()) +
labs(fill = "Genotype")
# plot for each individual:
ggplot(observed, aes(fill=gt, y=prop, x=type)) +
geom_bar(position="fill", stat="identity", ) +
xlab("Observed") +
ylab("Proportion of F0 genotype") +
labs(title = "F0 foundress of F1 male (0/1) x female (0/1)") +
geom_text(data = filter(observed, gt == "1/1"), aes(y=prop, x=type, label=total), vjust = 0) +
facet_wrap(~ sample) +
labs(fill = "Genotype") +
theme_classic() +
theme(axis.text.x = element_blank(), axis.title.x = element_blank())
We did 2 types of a-parametric tests for independence, testing the predicted counts against the observed counts:
for one family, 284
for one family, 338
for one family, 110
all genotypes for the 3 families
old code for false sites
# set the families (include only families with at least one adult F2)
family = grep("grndat|grnson",gt$ID, value=TRUE) %>%
str_extract("[^_]+") %>%
unique()
# define a list to put all the data frames in
list_grndat <- list()
for (fam in family) {
list_grndat[[fam]] <- table %>%
dplyr::select(starts_with(fam)) %>%
dplyr::filter_at(vars(matches("_son")), all_vars(. == "0/0")) %>%
dplyr::filter_at(vars(matches("_dat")), all_vars(. == "0/0")) %>%
dplyr::select(contains("grndat")) %>%
rownames_to_column("site")
}
#%>% dplyr::filter(across(everything(), ~ !str_detect(., "0/0")))%>% na.omit()
list_grndat$`240` %>% view()
# for male F2
list_grnson <- list()
for (fam in family) {
list_grnson[[fam]] <- table %>%
dplyr::select(starts_with(fam)) %>%
dplyr::filter_at(vars(matches("_son")), all_vars(. == "0/0")) %>%
dplyr::filter_at(vars(matches("_dat")), all_vars(. == "0/0")) %>%
dplyr::select(contains("grnson")) %>%
rownames_to_column("site")
}
#%>% dplyr::filter(across(everything(), ~ !str_detect(., "0/0"))) %>% na.omit()
list_grnson$`600`
## site 600_601a_grnson
## 1 NW_019211454.1_4912263 0/0
## 2 NW_019211454.1_4920532 0/0
## 3 NW_019211454.1_4921004 0/0
## 4 NW_019211454.1_4952477 0/0
## 5 NW_019211454.1_4966874 <NA>
## 6 NW_019211454.1_4975120 <NA>
## 7 NW_019211454.1_4978792 0/0
## 8 NW_019211454.1_4980011 <NA>
## 9 NW_019211454.1_5013515 0/0
## 10 NW_019211454.1_5021396 0/0
## 11 NW_019211454.1_5095835 0/0
## 12 NW_019211454.1_5122236 <NA>
## 13 NW_019211454.1_5166671 <NA>
## 14 NW_019211454.1_5175616 <NA>
## 15 NW_019211454.1_5184381 0/0
## 16 NW_019211454.1_5214324 0/0
## 17 NW_019211454.1_5229776 <NA>
## 18 NW_019211454.1_5257386 0/0
## 19 NW_019211454.1_5261860 <NA>
## 20 NW_019211454.1_5266356 0/0
## 21 NW_019211454.1_5268357 0/0
## 22 NW_019211454.1_5279110 0/0
## 23 NW_019211454.1_5918916 0/0
## 24 NW_019211454.1_6033069 0/0
## 25 NW_019211454.1_6041909 0/0
## 26 NW_019211454.1_6167320 <NA>
## 27 NW_019211454.1_6226168 0/0
## 28 NW_019211454.1_6251092 0/0
## 29 NW_019211454.1_6259775 <NA>
## 30 NW_019211454.1_6366438 <NA>
## 31 NW_019211454.1_6453874 <NA>
## 32 NW_019211454.1_6479391 0/0
## 33 NW_019211454.1_6503569 0/0
## 34 NW_019211454.1_6509668 <NA>
## 35 NW_019211454.1_6512550 0/0
## 36 NW_019211454.1_6512723 <NA>
## 37 NW_019211454.1_6516631 <NA>
## 38 NW_019211454.1_6538442 <NA>
## 39 NW_019211454.1_8203978 <NA>
## 40 NW_019211454.1_9513244 <NA>
## 41 NW_019211454.1_9526241 0/0
## 42 NW_019211454.1_9526415 0/0
## 43 NW_019211454.1_10787954 0/0
## 44 NW_019211454.1_10789989 0/0
## 45 NW_019211454.1_12645630 0/0
## 46 NW_019211454.1_12703153 0/0
## 47 NW_019211454.1_12852621 <NA>
## 48 NW_019211454.1_12860342 <NA>
## 49 NW_019211454.1_12869979 <NA>
## 50 NW_019211454.1_12874184 <NA>
## 51 NW_019211454.1_12889655 <NA>
## 52 NW_019211454.1_12943805 <NA>
## 53 NW_019211454.1_13005638 0/0
## 54 NW_019211454.1_13010879 0/0
## 55 NW_019211454.1_13013741 <NA>
## 56 NW_019211454.1_13044739 <NA>
## 57 NW_019211454.1_13056888 <NA>
## 58 NW_019211454.1_13063554 0/0
## 59 NW_019211454.1_13073686 <NA>
## 60 NW_019211454.1_13100815 <NA>
## 61 NW_019211454.1_13107115 0/0
## 62 NW_019211454.1_13146372 <NA>
## 63 NW_019211454.1_13160204 <NA>
## 64 NW_019211454.1_13177229 <NA>
## 65 NW_019211454.1_13186306 <NA>
## 66 NW_019211454.1_13204249 0/0
## 67 NW_019211454.1_14865144 0/0
## 68 NW_019211454.1_14955656 0/0
## 69 NW_019211454.1_14999663 <NA>
## 70 NW_019211454.1_14999711 <NA>
## 71 NW_019211454.1_15001520 0/0
## 72 NW_019211454.1_15003776 <NA>
## 73 NW_019211454.1_15011010 <NA>
## 74 NW_019211454.1_15037118 <NA>
## 75 NW_019211454.1_15153923 <NA>
## 76 NW_019211454.1_15158713 0/0
## 77 NW_019211454.1_15161035 0/0
## 78 NW_019211454.1_15168967 0/0
## 79 NW_019211454.1_15181015 0/0
## 80 NW_019211454.1_15610546 <NA>
## 81 NW_019211454.1_15648222 0/0
## 82 NW_019211454.1_15648701 0/0
## 83 NW_019211454.1_15701177 <NA>
## 84 NW_019211454.1_15712292 0/0
## 85 NW_019211454.1_16244152 0/0
## 86 NW_019211454.1_16246384 0/0
## 87 NW_019211454.1_16252347 <NA>
## 88 NW_019211454.1_16252834 0/0
## 89 NW_019211454.1_16266943 0/0
## 90 NW_019211454.1_16297722 <NA>
## 91 NW_019211454.1_16312031 0/0
## 92 NW_019211454.1_16312209 0/0
## 93 NW_019211454.1_16340860 <NA>
## 94 NW_019211454.1_16348147 <NA>
## 95 NW_019211454.1_16350959 0/0
## 96 NW_019211454.1_16353460 <NA>
## 97 NW_019211454.1_16357130 0/0
## 98 NW_019211454.1_16368943 <NA>
## 99 NW_019211454.1_16380802 <NA>
## 100 NW_019211454.1_16387023 0/0
## 101 NW_019211454.1_16389486 0/0
## 102 NW_019211454.1_16410682 <NA>
## 103 NW_019211454.1_16419671 0/0
## 104 NW_019211454.1_16461135 <NA>
## 105 NW_019211454.1_16469261 0/0
## 106 NW_019211454.1_16472031 0/0
## 107 NW_019211454.1_16483028 0/0
## 108 NW_019211454.1_16502287 0/0
## 109 NW_019211454.1_16519175 <NA>
## 110 NW_019211454.1_16548057 0/0
## 111 NW_019211454.1_16548075 0/0
## 112 NW_019211454.1_16548390 0/0
## 113 NW_019211454.1_16548759 0/0
## 114 NW_019211454.1_16549009 0/0
## 115 NW_019211454.1_16550594 <NA>
## 116 NW_019211454.1_16555740 <NA>
## 117 NW_019211454.1_16557116 0/0
## 118 NW_019211454.1_16557423 0/0
## 119 NW_019211454.1_16557649 0/0
## 120 NW_019211454.1_16561544 <NA>
## 121 NW_019211454.1_16574668 <NA>
## 122 NW_019211454.1_16580293 <NA>
## 123 NW_019211454.1_16645606 0/0
## 124 NW_019211454.1_16647761 <NA>
## 125 NW_019211454.1_16647867 0/0
## 126 NW_019211454.1_16655076 0/0
## 127 NW_019211454.1_16656654 <NA>
## 128 NW_019211454.1_16668875 0/0
## 129 NW_019211454.1_16679234 <NA>
## 130 NW_019211454.1_16689923 0/0
## 131 NW_019211454.1_16693222 0/0
## 132 NW_019211454.1_16698054 0/0
## 133 NW_019211454.1_16969415 0/0
## 134 NW_019211454.1_16978613 0/0
## 135 NW_019211454.1_16995141 0/0
## 136 NW_019211454.1_17114649 <NA>
## 137 NW_019211454.1_17129653 <NA>
## 138 NW_019211454.1_17146268 <NA>
## 139 NW_019211454.1_17148834 <NA>
## 140 NW_019211454.1_17182202 <NA>
## 141 NW_019211454.1_17193528 <NA>
## 142 NW_019211454.1_17198346 0/0
## 143 NW_019211454.1_17198496 <NA>
## 144 NW_019211454.1_17211037 0/0
## 145 NW_019211454.1_17217149 <NA>
## 146 NW_019211454.1_17219103 0/0
## 147 NW_019211454.1_17223506 0/0
## 148 NW_019211454.1_17243529 0/0
## 149 NW_019211454.1_17246645 <NA>
## 150 NW_019211454.1_17248571 0/0
## 151 NW_019211454.1_17257077 0/0
## 152 NW_019211454.1_17260462 0/0
## 153 NW_019211454.1_17271714 0/0
## 154 NW_019211454.1_17290015 <NA>
## 155 NW_019211454.1_17291980 <NA>
## 156 NW_019211454.1_17301347 0/0
## 157 NW_019211454.1_17314319 0/0
## 158 NW_019211454.1_17320038 <NA>
## 159 NW_019211454.1_17322790 0/0
## 160 NW_019211454.1_17331569 0/0
## 161 NW_019211454.1_17344167 <NA>
## 162 NW_019211454.1_17350818 <NA>
## 163 NW_019211454.1_17352533 <NA>
## 164 NW_019211454.1_17374556 0/0
## 165 NW_019211454.1_17374579 0/0
## 166 NW_019211454.1_17375864 0/0
## 167 NW_019211454.1_17376941 <NA>
## 168 NW_019211454.1_17541597 0/0
## 169 NW_019211454.1_17551873 0/0
## 170 NW_019211454.1_17558887 0/0
## 171 NW_019211454.1_17576766 0/0
## 172 NW_019211454.1_17578196 0/0
## 173 NW_019211454.1_17583341 0/0
## 174 NW_019211454.1_17606030 0/0
## 175 NW_019211454.1_17610400 0/0
## 176 NW_019211454.1_17622107 0/0
## 177 NW_019211454.1_17772729 <NA>
## 178 NW_019211454.1_17803877 <NA>
## 179 NW_019211454.1_17837365 0/0
## 180 NW_019211454.1_17837385 0/0
## 181 NW_019211454.1_17837976 0/0
## 182 NW_019211454.1_17844163 0/0
## 183 NW_019211454.1_17954895 <NA>
## 184 NW_019211454.1_17955731 0/0
## 185 NW_019211454.1_17955924 <NA>
## 186 NW_019211454.1_17965364 <NA>
## 187 NW_019211454.1_17982176 0/0
## 188 NW_019211454.1_17986906 0/0
## 189 NW_019211454.1_17999036 <NA>
## 190 NW_019211454.1_18109983 <NA>
## 191 NW_019211454.1_18113968 <NA>
## 192 NW_019211454.1_18906639 <NA>
## 193 NW_019211454.1_21189838 <NA>
## 194 NW_019211454.1_21205822 <NA>
## 195 NW_019211454.1_23149732 <NA>
## 196 NW_019211454.1_23162575 <NA>
## 197 NW_019211454.1_23188311 0/0
## 198 NW_019211454.1_23192907 <NA>
## 199 NW_019211454.1_23199148 <NA>
## 200 NW_019211454.1_23212140 <NA>
## 201 NW_019211454.1_23218946 0/0
## 202 NW_019211454.1_23251760 0/0
## 203 NW_019211454.1_23254977 0/0
## 204 NW_019211454.1_23264224 0/0
## 205 NW_019211454.1_23268905 <NA>
## 206 NW_019211454.1_23272086 <NA>
## 207 NW_019211454.1_23276013 0/0
## 208 NW_019211454.1_23278009 0/0
## 209 NW_019211454.1_23298719 <NA>
## 210 NW_019211454.1_23299008 <NA>
## 211 NW_019211454.1_23300769 <NA>
## 212 NW_019211454.1_23300881 <NA>
## 213 NW_019211454.1_23301117 0/0
## 214 NW_019211454.1_23333824 <NA>
## 215 NW_019211454.1_23366149 0/0
## 216 NW_019211454.1_23382553 0/0
## 217 NW_019211454.1_23391385 0/0
## 218 NW_019211454.1_23410544 <NA>
## 219 NW_019211454.1_23425501 0/0
## 220 NW_019211454.1_23442161 0/0
## 221 NW_019211454.1_23577877 <NA>
## 222 NW_019211454.1_23595719 <NA>
## 223 NW_019211454.1_23624888 <NA>
## 224 NW_019211454.1_23634167 <NA>
## 225 NW_019211454.1_23637875 0/0
## 226 NW_019211454.1_23644002 <NA>
## 227 NW_019211454.1_24570146 <NA>
## 228 NW_019211454.1_24570473 <NA>
## 229 NW_019211454.1_24578003 <NA>
## 230 NW_019211454.1_24683480 0/0
## 231 NW_019211454.1_24727073 0/0
## 232 NW_019211454.1_24775786 0/0
## 233 NW_019211454.1_24819583 0/0
## 234 NW_019211454.1_27474697 <NA>
## 235 NW_019211454.1_27476484 0/0
## 236 NW_019211454.1_27491198 0/0
## 237 NW_019211454.1_27498234 0/0
## 238 NW_019211454.1_27498752 0/0
## 239 NW_019211454.1_27520650 0/0
## 240 NW_019211454.1_27531194 <NA>
## 241 NW_019211454.1_27540325 <NA>
## 242 NW_019211454.1_27580634 0/0
## 243 NW_019211454.1_27581718 0/0
## 244 NW_019211454.1_27582068 0/0
## 245 NW_019211454.1_27593680 0/0
## 246 NW_019211454.1_27594931 <NA>
## 247 NW_019211454.1_27596087 0/0
## 248 NW_019211454.1_27596126 0/0
## 249 NW_019211454.1_27612215 0/0
## 250 NW_019211454.1_27617579 0/0
## 251 NW_019211454.1_27699731 <NA>
## 252 NW_019211454.1_27727090 <NA>
## 253 NW_019211454.1_27735379 <NA>
## 254 NW_019211454.1_27738749 0/0
## 255 NW_019211454.1_27751954 0/0
## 256 NW_019211454.1_27753097 0/0
## 257 NW_019211454.1_27758233 <NA>
## 258 NW_019211454.1_27758304 0/0
## 259 NW_019211454.1_27758399 0/0
## 260 NW_019211454.1_27790477 <NA>
## 261 NW_019211454.1_27795450 <NA>
## 262 NW_019211454.1_27816715 <NA>
## 263 NW_019211454.1_27892109 0/0
## 264 NW_019211454.1_27892155 0/0
## 265 NW_019211454.1_27959399 0/0
## 266 NW_019211454.1_27961020 <NA>
## 267 NW_019211454.1_27996384 0/0
## 268 NW_019211454.1_28024213 0/0
## 269 NW_019211454.1_28039755 <NA>
## 270 NW_019211454.1_28065845 0/0
## 271 NW_019211454.1_28066125 <NA>
## 272 NW_019211454.1_29231990 0/0
## 273 NW_019211454.1_29232287 0/0
## 274 NW_019211454.1_29261601 0/0
## 275 NW_019211454.1_29919924 0/0
## 276 NW_019211454.1_30101274 <NA>
## 277 NW_019211454.1_30402966 0/0
## 278 NW_019211454.1_30419096 <NA>
## 279 NW_019211454.1_30426006 0/0
## 280 NW_019211454.1_30431348 0/0
## 281 NW_019211454.1_30435231 <NA>
## 282 NW_019211454.1_30447033 <NA>
## 283 NW_019211454.1_30454380 <NA>
## 284 NW_019211454.1_30454704 <NA>
## 285 NW_019211454.1_30457049 0/0
## 286 NW_019211454.1_30470808 <NA>
## 287 NW_019211454.1_30480822 0/0
## 288 NW_019211454.1_30498655 <NA>
## 289 NW_019211454.1_30502656 <NA>
## 290 NW_019211454.1_30505451 0/0
## 291 NW_019211454.1_30519580 0/0
## 292 NW_019211454.1_30520762 <NA>
## 293 NW_019211454.1_30527028 0/0
## 294 NW_019211454.1_30527215 0/0
## 295 NW_019211454.1_30528862 <NA>
## 296 NW_019211454.1_30756276 <NA>
## 297 NW_019211454.1_30781127 0/0
## 298 NW_019211454.1_30802390 <NA>
## 299 NW_019211454.1_30815022 0/0
## 300 NW_019211454.1_30850228 0/0
## 301 NW_019211454.1_31019741 0/0
## 302 NW_019211454.1_31117387 0/0
## 303 NW_019211454.1_32039254 0/0
## 304 NW_019211454.1_32054483 0/0
## 305 NW_019211454.1_32099381 <NA>
## 306 NW_019211454.1_32136801 <NA>
## 307 NW_019211454.1_32140314 <NA>
## 308 NW_019211454.1_32142383 0/0
## 309 NW_019211454.1_32143449 0/0
## 310 NW_019211454.1_32145009 0/0
## 311 NW_019211454.1_32156007 <NA>
## 312 NW_019211454.1_32166326 <NA>
## 313 NW_019211454.1_32171625 0/0
## 314 NW_019211454.1_32175316 0/0
## 315 NW_019211454.1_32186517 <NA>
## 316 NW_019211454.1_32209139 <NA>
## 317 NW_019211454.1_32278692 <NA>
## 318 NW_019211454.1_32293756 0/0
## 319 NW_019211454.1_32297314 <NA>
## 320 NW_019211454.1_32317836 0/0
## 321 NW_019211454.1_32964012 0/0
## 322 NW_019211454.1_33338832 0/0
## 323 NW_019211454.1_33352032 0/0
## 324 NW_019211454.1_33356322 0/0
## 325 NW_019211454.1_33357679 <NA>
## 326 NW_019211454.1_33372762 <NA>
## 327 NW_019211454.1_33421667 <NA>
## 328 NW_019211454.1_33448354 <NA>
## 329 NW_019211454.1_33452617 0/0
## 330 NW_019211454.1_33457401 0/0
## 331 NW_019211454.1_33470836 0/0
## 332 NW_019211454.1_33481705 0/0
## 333 NW_019211454.1_34661141 0/0
## 334 NW_019211454.1_34696147 <NA>
## 335 NW_019211454.1_34710662 0/0
## 336 NW_019211454.1_34735114 0/0
## 337 NW_019211454.1_34750224 <NA>
## 338 NW_019211454.1_34755204 <NA>
## 339 NW_019211454.1_34767703 0/0
## 340 NW_019211454.1_34769592 0/0
## 341 NW_019211454.1_34780333 0/0
## 342 NW_019211454.1_34780346 0/0
## 343 NW_019211454.1_34809110 0/0
## 344 NW_019211454.1_34843865 <NA>
## 345 NW_019211454.1_34871348 <NA>
## 346 NW_019211454.1_34947904 0/0
## 347 NW_019211454.1_35028182 <NA>
## 348 NW_019211454.1_35029455 <NA>
## 349 NW_019211454.1_35038046 <NA>
## 350 NW_019211454.1_35081471 <NA>
## 351 NW_019211454.1_35300566 0/0
## 352 NW_019211454.1_35531494 <NA>
## 353 NW_019211454.1_35537358 0/0
## 354 NW_019211454.1_35635742 <NA>
## 355 NW_019211454.1_35636953 0/0
## 356 NW_019211454.1_35657554 0/0
## 357 NW_019211454.1_36225107 0/0
## 358 NW_019211454.1_36231049 0/0
## 359 NW_019211454.1_36249867 0/0
## 360 NW_019211454.1_36303013 <NA>
## 361 NW_019211454.1_36303334 <NA>
## 362 NW_019211454.1_36305588 <NA>
## 363 NW_019211454.1_36305722 0/0
## 364 NW_019211454.1_36306986 <NA>
## 365 NW_019211454.1_36309014 0/0
## 366 NW_019211454.1_36312828 <NA>
## 367 NW_019211454.1_36322719 0/0
## 368 NW_019211454.1_36325864 0/0
## 369 NW_019211454.1_36327293 0/0
## 370 NW_019211454.1_36400301 <NA>
## 371 NW_019211454.1_36429883 0/0
## 372 NW_019211454.1_36431797 0/0
## 373 NW_019211454.1_36431811 0/0
## 374 NW_019211454.1_36434432 0/0
## 375 NW_019211454.1_36437879 <NA>
## 376 NW_019211454.1_36439773 0/0
## 377 NW_019211454.1_36456925 0/0
## 378 NW_019211454.1_36467442 0/0
## 379 NW_019211454.1_36476993 <NA>
## 380 NW_019211454.1_36485502 <NA>
## 381 NW_019211454.1_36487316 <NA>
## 382 NW_019211454.1_36488320 0/0
## 383 NW_019211454.1_36489299 <NA>
## 384 NW_019211454.1_36490354 0/0
## 385 NW_019211454.1_36492502 <NA>
## 386 NW_019211454.1_36517090 0/0
## 387 NW_019211454.1_36559306 <NA>
## 388 NW_019211454.1_36579567 <NA>
## 389 NW_019211454.1_36581439 0/0
## 390 NW_019211454.1_36592385 <NA>
## 391 NW_019211454.1_36594808 <NA>
## 392 NW_019211454.1_36625536 0/0
## 393 NW_019211454.1_36627365 0/0
## 394 NW_019211454.1_36651633 0/0
## 395 NW_019211454.1_36707318 0/0
## 396 NW_019211454.1_36757145 <NA>
## 397 NW_019211454.1_36761381 0/0
## 398 NW_019211454.1_36779314 <NA>
## 399 NW_019211454.1_36785002 0/0
## 400 NW_019211454.1_36786796 0/0
## 401 NW_019211454.1_36791692 0/1
## 402 NW_019211454.1_36799139 0/0
## 403 NW_019211454.1_36827402 0/0
## 404 NW_019211454.1_36827809 0/0
## 405 NW_019211454.1_36848731 <NA>
## 406 NW_019211454.1_36852527 <NA>
## 407 NW_019211454.1_36854992 <NA>
## 408 NW_019211454.1_36859160 0/0
## 409 NW_019211454.1_36860378 <NA>
## 410 NW_019211454.1_36866653 <NA>
## 411 NW_019211454.1_36870521 0/0
## 412 NW_019211454.1_36873027 <NA>
## 413 NW_019211454.1_36874395 0/0
## 414 NW_019211454.1_36881402 <NA>
## 415 NW_019211454.1_36885217 0/0
## 416 NW_019211454.1_36922811 <NA>
## 417 NW_019211454.1_36935960 0/0
## 418 NW_019211454.1_36941774 <NA>
## 419 NW_019211454.1_36961149 <NA>
## 420 NW_019211454.1_36963866 <NA>
## 421 NW_019211454.1_36970514 0/0
## 422 NW_019211454.1_36974827 0/0
## 423 NW_019211454.1_37003339 <NA>
## 424 NW_019211454.1_37003387 <NA>
## 425 NW_019211454.1_37010589 <NA>
## 426 NW_019211454.1_37011724 <NA>
## 427 NW_019211454.1_37024601 0/0
## 428 NW_019211454.1_37025330 <NA>
## 429 NW_019211454.1_37032799 <NA>
## 430 NW_019211454.1_37039254 <NA>
## 431 NW_019211454.1_37047454 0/0
## 432 NW_019211454.1_37049914 0/0
## 433 NW_019211454.1_37056997 <NA>
## 434 NW_019211454.1_37070881 0/0
## 435 NW_019211454.1_37100466 <NA>
## 436 NW_019211454.1_37128833 0/0
## 437 NW_019211454.1_37130569 <NA>
## 438 NW_019211454.1_37133936 0/0
## 439 NW_019211454.1_37134076 0/0
## 440 NW_019211454.1_37135404 0/0
## 441 NW_019211454.1_37135625 0/0
## 442 NW_019211454.1_37138864 <NA>
## 443 NW_019211454.1_37138887 <NA>
## 444 NW_019211454.1_37143549 0/0
## 445 NW_019211454.1_37753628 <NA>
## 446 NW_019211454.1_37910377 <NA>
## 447 NW_019211454.1_37938433 <NA>
## 448 NW_019211454.1_37940118 <NA>
## 449 NW_019211454.1_37940348 0/0
## 450 NW_019211454.1_37947399 0/0
## 451 NW_019211454.1_37948672 <NA>
## 452 NW_019211454.1_37987970 <NA>
## 453 NW_019211454.1_37987998 <NA>
## 454 NW_019211454.1_37988822 <NA>
## 455 NW_019211454.1_37989288 <NA>
## 456 NW_019211454.1_37992108 0/0
## 457 NW_019211454.1_37993898 <NA>
## 458 NW_019211454.1_37994047 0/0
## 459 NW_019211454.1_37994181 <NA>
## 460 NW_019211454.1_37998762 <NA>
## 461 NW_019211454.1_38000338 <NA>
## 462 NW_019211454.1_38001316 <NA>
## 463 NW_019211454.1_38001809 0/0
## 464 NW_019211454.1_38004005 <NA>
## 465 NW_019211454.1_38008590 <NA>
## 466 NW_019211454.1_38046727 <NA>
## 467 NW_019211454.1_38624087 0/0
## 468 NW_019211454.1_38639387 0/0
## 469 NW_019211454.1_38642178 0/0
## 470 NW_019211454.1_38646340 0/0
## 471 NW_019211454.1_38653750 0/0
## 472 NW_019211454.1_38693611 <NA>
## 473 NW_019211454.1_38694566 <NA>
## 474 NW_019211454.1_38710500 0/0
## 475 NW_019211454.1_38719649 0/0
## 476 NW_019211454.1_38719683 0/0
## 477 NW_019211454.1_38719882 <NA>
## 478 NW_019211454.1_38723980 0/0
## 479 NW_019211454.1_38728405 <NA>
## 480 NW_019211454.1_38731284 <NA>
## 481 NW_019211454.1_38732279 0/0
## 482 NW_019211454.1_38734590 0/0
## 483 NW_019211454.1_38760743 0/0
## 484 NW_019211454.1_38760785 0/0
## 485 NW_019211454.1_38760800 0/0
## 486 NW_019211454.1_38772689 0/0
## 487 NW_019211454.1_39463236 0/0
## 488 NW_019211454.1_39511091 <NA>
## 489 NW_019211454.1_39719447 0/0
## 490 NW_019211454.1_39762977 0/0
## 491 NW_019211454.1_39788462 0/0
## 492 NW_019211454.1_39798196 <NA>
## 493 NW_019211454.1_39808695 <NA>
## 494 NW_019211454.1_39810900 <NA>
## 495 NW_019211454.1_40034273 <NA>
## 496 NW_019211454.1_40070100 0/0
## 497 NW_019211454.1_40070316 0/0
## 498 NW_019211454.1_40089652 0/0
## 499 NW_019211454.1_40090061 0/0
## 500 NW_019211454.1_40114751 <NA>
## 501 NW_019211454.1_40126972 0/0
## 502 NW_019211454.1_40129662 0/0
## 503 NW_019211454.1_40132976 0/0
## 504 NW_019211454.1_40144234 <NA>
## 505 NW_019211454.1_40147921 <NA>
## 506 NW_019211454.1_40147937 <NA>
## 507 NW_019211454.1_40157895 0/0
## 508 NW_019211454.1_40158668 0/0
## 509 NW_019211454.1_40158739 0/0
## 510 NW_019211454.1_40162517 0/0
## 511 NW_019211454.1_40167978 0/0
## 512 NW_019211454.1_40169514 0/0
## 513 NW_019211454.1_40170666 0/0
## 514 NW_019211454.1_40187999 0/0
## 515 NW_019211454.1_40206130 0/0
## 516 NW_019211454.1_40208386 0/0
## 517 NW_019211454.1_40211570 0/0
## 518 NW_019211454.1_40227396 <NA>
## 519 NW_019211454.1_40241170 0/0
## 520 NW_019211454.1_40243531 <NA>
## 521 NW_019211454.1_40250427 0/0
## 522 NW_019211454.1_40261072 0/0
## 523 NW_019211454.1_40266318 0/0
## 524 NW_019211454.1_40276818 <NA>
## 525 NW_019211454.1_40280330 0/0
## 526 NW_019211454.1_40626334 0/0
## 527 NW_019211454.1_40630404 <NA>
## 528 NW_019211454.1_40631887 0/0
## 529 NW_019211454.1_40638403 0/0
## 530 NW_019211454.1_40666283 0/0
## 531 NW_019211454.1_40680609 0/0
## 532 NW_019211454.1_40683288 <NA>
## 533 NW_019211454.1_40693798 <NA>
## 534 NW_019211454.1_40698206 0/0
## 535 NW_019211454.1_40717536 0/0
## 536 NW_019211454.1_40721776 0/0
## 537 NW_019211454.1_40722145 <NA>
## 538 NW_019211454.1_40725411 0/0
## 539 NW_019211454.1_40725787 0/0
## 540 NW_019211454.1_41101184 <NA>
## 541 NW_019211454.1_41155118 0/0
## 542 NW_019211454.1_41157815 0/0
## 543 NW_019211454.1_41158461 <NA>
## 544 NW_019211454.1_41176200 0/0
## 545 NW_019211454.1_41178977 0/0
## 546 NW_019211454.1_41187641 <NA>
## 547 NW_019211454.1_41189701 0/0
## 548 NW_019211454.1_41195384 0/0
## 549 NW_019211454.1_41198484 <NA>
## 550 NW_019211454.1_41200614 0/0
## 551 NW_019211454.1_41200809 0/1
## 552 NW_019211454.1_41201413 0/0
## 553 NW_019211454.1_41213762 <NA>
## 554 NW_019211454.1_41220506 0/0
## 555 NW_019211454.1_41230368 0/0
## 556 NW_019211454.1_41240538 0/0
## 557 NW_019211454.1_41244157 0/0
## 558 NW_019211454.1_41250910 <NA>
## 559 NW_019211454.1_41252294 0/0
## 560 NW_019211454.1_41271694 0/0
## 561 NW_019211454.1_41275567 <NA>
## 562 NW_019211454.1_41276956 <NA>
## 563 NW_019211454.1_41276969 0/0
## 564 NW_019211454.1_41347283 <NA>
## 565 NW_019211454.1_41379837 <NA>
## 566 NW_019211454.1_41389149 0/0
## 567 NW_019211454.1_44066323 0/0
## 568 NW_019211454.1_44075125 0/0
## 569 NW_019211454.1_44103469 0/0
## 570 NW_019211454.1_44105552 0/0
## 571 NW_019211454.1_44112775 <NA>
## 572 NW_019211454.1_44140420 0/0
## 573 NW_019211454.1_46214928 <NA>
## 574 NW_019211454.1_48486198 0/0
## 575 NW_019211454.1_48491589 0/0
## 576 NW_019211454.1_48525461 <NA>
## 577 NW_019211454.1_49033072 0/0
## 578 NW_019211454.1_49033419 0/0
## 579 NW_019211454.1_49035336 <NA>
## 580 NW_019211454.1_49043289 0/0
## 581 NW_019211454.1_49057108 <NA>
## 582 NW_019211454.1_49378114 0/0
## 583 NW_019211454.1_50191124 0/0
## 584 NW_019211454.1_52386002 <NA>
## 585 NW_019211454.1_53431088 <NA>
## 586 NW_019211454.1_56129168 <NA>
## 587 NW_019211454.1_56754553 <NA>
## 588 NW_019211454.1_56790992 0/0
## 589 NW_019211454.1_58461103 0/0
## 590 NW_019211454.1_60040889 0/0
## 591 NW_019211454.1_60045747 0/0
## 592 NW_019211454.1_60052833 0/0
## 593 NW_019211454.1_60070105 0/0
## 594 NW_019211454.1_60084987 <NA>
## 595 NW_019211454.1_60086106 0/0
## 596 NW_019211454.1_60090982 <NA>
## 597 NW_019211454.1_60091330 <NA>
## 598 NW_019211454.1_60094151 0/0
## 599 NW_019211454.1_60094842 0/0
## 600 NW_019211454.1_60096936 <NA>
## 601 NW_019211454.1_60116457 0/0
## 602 NW_019211454.1_60116528 <NA>
## 603 NW_019211454.1_60119713 0/0
## 604 NW_019211454.1_60120460 <NA>
## 605 NW_019211454.1_60126262 <NA>
## 606 NW_019211454.1_60143176 <NA>
## 607 NW_019211454.1_60150897 <NA>
## 608 NW_019211454.1_60178344 0/0
## 609 NW_019211454.1_60199503 0/0
## 610 NW_019211454.1_60214441 <NA>
## 611 NW_019211454.1_60214456 <NA>
## 612 NW_019211454.1_60227203 0/0
## 613 NW_019211454.1_60233060 0/0
## 614 NW_019211454.1_60294325 0/0
## 615 NW_019211454.1_60295124 <NA>
## 616 NW_019211454.1_60300810 <NA>
## 617 NW_019211454.1_60300845 <NA>
## 618 NW_019211454.1_60302269 0/0
## 619 NW_019211454.1_60309499 0/0
## 620 NW_019211454.1_60314137 <NA>
## 621 NW_019211454.1_60870511 <NA>
## 622 NW_019211454.1_60871712 <NA>
## 623 NW_019211454.1_60886131 0/0
## 624 NW_019211454.1_60892756 <NA>
## 625 NW_019211454.1_60898037 <NA>
## 626 NW_019211454.1_60903426 0/0
## 627 NW_019211454.1_60910004 0/0
## 628 NW_019211454.1_60918838 0/0
## 629 NW_019211454.1_60936792 0/0
## 630 NW_019211454.1_60999522 0/0
## 631 NW_019211454.1_61005443 <NA>
## 632 NW_019211454.1_61010017 0/0
## 633 NW_019211454.1_61036882 0/0
## 634 NW_019211454.1_61059853 0/0
## 635 NW_019211454.1_61060838 0/0
## 636 NW_019211454.1_61060872 0/0
## 637 NW_019211454.1_61062560 0/0
## 638 NW_019211454.1_61064207 0/0
## 639 NW_019211454.1_61066588 0/0
## 640 NW_019211454.1_61071808 <NA>
## 641 NW_019211454.1_61100520 <NA>
## 642 NW_019211454.1_62221894 0/0
## 643 NW_019211454.1_62223198 0/0
## 644 NW_019211454.1_62228815 <NA>
## 645 NW_019211454.1_62253745 <NA>
## 646 NW_019211454.1_62297656 0/0
## 647 NW_019211454.1_62362013 0/0
## 648 NW_019211454.1_62379659 <NA>
## 649 NW_019211454.1_62384719 <NA>
## 650 NW_019211454.1_62742862 <NA>
## 651 NW_019211454.1_62744523 0/0
## 652 NW_019211454.1_62755203 0/0
## 653 NW_019211454.1_62760791 <NA>
## 654 NW_019211454.1_62766087 0/0
## 655 NW_019211454.1_62845696 <NA>
## 656 NW_019211454.1_66357056 <NA>
## 657 NW_019211454.1_67642464 <NA>
## 658 NW_019211454.1_68070552 <NA>
## 659 NW_019211454.1_68081673 <NA>
## 660 NW_019211454.1_68208481 0/0
## 661 NW_019211454.1_68210796 0/0
## 662 NW_019211454.1_68249627 <NA>
## 663 NW_019211454.1_68249776 0/0
## 664 NW_019211454.1_68250522 <NA>
## 665 NW_019211454.1_68252004 0/0
## 666 NW_019211454.1_68252460 <NA>
## 667 NW_019211454.1_68268811 <NA>
## 668 NW_019211454.1_68303585 <NA>
## 669 NW_019211454.1_68310183 0/0
## 670 NW_019211454.1_68334629 <NA>
## 671 NW_019211454.1_68385302 <NA>
## 672 NW_019211454.1_68390801 <NA>
## 673 NW_019211454.1_68392311 <NA>
## 674 NW_019211454.1_68410266 0/0
## 675 NW_019211454.1_68426227 0/0
## 676 NW_019211454.1_68436675 0/0
## 677 NW_019211454.1_68524498 0/0
## 678 NW_019211454.1_68539594 0/0
## 679 NW_019211454.1_68539621 0/0
## 680 NW_019211454.1_68565197 0/0
## 681 NW_019211454.1_69692222 0/0
## 682 NW_019211454.1_69694549 <NA>
## 683 NW_019211454.1_69694715 0/0
## 684 NW_019211454.1_69702475 0/0
## 685 NW_019211454.1_69765615 0/0
## 686 NW_019211454.1_69772933 0/0
## 687 NW_019211454.1_69777852 0/0
## 688 NW_019211454.1_69795954 <NA>
## 689 NW_019211454.1_70186883 0/0
## 690 NW_019211454.1_72789477 0/0
## 691 NW_019211454.1_72789486 0/0
## 692 NW_019211454.1_72789502 0/0
## 693 NW_019211454.1_72789511 0/0
## 694 NW_019211454.1_74126558 <NA>
## 695 NW_019211454.1_74128294 <NA>
## 696 NW_019211454.1_74290774 0/0
## 697 NW_019211454.1_74292728 0/0
## 698 NW_019211454.1_74341144 <NA>
## 699 NW_019211454.1_74352633 0/0
## 700 NW_019211454.1_74375274 <NA>
## 701 NW_019211454.1_74383872 <NA>
## 702 NW_019211454.1_74386880 0/0
## 703 NW_019211454.1_74397070 <NA>
## 704 NW_019211454.1_74668449 0/0
## 705 NW_019211454.1_75034547 <NA>
## 706 NW_019211454.1_75046501 <NA>
## 707 NW_019211454.1_75059444 <NA>
## 708 NW_019211454.1_75090669 <NA>
## 709 NW_019211454.1_75101747 0/0
## 710 NW_019211454.1_75101827 0/0
## 711 NW_019211454.1_75101887 0/0
## 712 NW_019211454.1_75112840 0/0
## 713 NW_019211454.1_75133554 <NA>
## 714 NW_019211454.1_75136801 0/0
## 715 NW_019211454.1_75147647 0/0
## 716 NW_019211454.1_75149807 <NA>
## 717 NW_019211454.1_75158568 0/0
## 718 NW_019211454.1_75163728 0/0
## 719 NW_019211454.1_75194641 <NA>
## 720 NW_019211454.1_75195645 <NA>
## 721 NW_019211454.1_76223546 <NA>
## 722 NW_019211454.1_76223565 0/0
## 723 NW_019211454.1_76234397 <NA>
## 724 NW_019211454.1_76591412 <NA>
## 725 NW_019211455.1_297151 <NA>
## 726 NW_019211455.1_1723309 0/0
## 727 NW_019211455.1_1735588 0/0
## 728 NW_019211455.1_1755178 0/0
## 729 NW_019211455.1_1757742 <NA>
## 730 NW_019211455.1_1759287 <NA>
## 731 NW_019211455.1_1786467 <NA>
## 732 NW_019211455.1_1790136 <NA>
## 733 NW_019211455.1_1796024 0/0
## 734 NW_019211455.1_1799665 0/0
## 735 NW_019211455.1_1817006 <NA>
## 736 NW_019211455.1_1817248 <NA>
## 737 NW_019211455.1_1818976 <NA>
## 738 NW_019211455.1_1827439 0/0
## 739 NW_019211455.1_1844582 <NA>
## 740 NW_019211455.1_1998846 0/0
## 741 NW_019211455.1_2002301 0/0
## 742 NW_019211455.1_2003141 <NA>
## 743 NW_019211455.1_2018035 0/0
## 744 NW_019211455.1_2024149 0/0
## 745 NW_019211455.1_2035966 <NA>
## 746 NW_019211455.1_2037845 <NA>
## 747 NW_019211455.1_2049751 <NA>
## 748 NW_019211455.1_2055201 0/0
## 749 NW_019211455.1_2086517 0/0
## 750 NW_019211455.1_2090967 0/0
## 751 NW_019211455.1_2105994 0/0
## 752 NW_019211455.1_2156648 0/0
## 753 NW_019211455.1_2269993 0/0
## 754 NW_019211455.1_2270186 0/0
## 755 NW_019211455.1_2273866 0/0
## 756 NW_019211455.1_6296845 0/0
## 757 NW_019211455.1_6297935 <NA>
## 758 NW_019211455.1_6310894 <NA>
## 759 NW_019211455.1_6311013 <NA>
## 760 NW_019211455.1_6313858 <NA>
## 761 NW_019211455.1_6328519 <NA>
## 762 NW_019211455.1_6341511 <NA>
## 763 NW_019211455.1_6342371 <NA>
## 764 NW_019211455.1_6388395 <NA>
## 765 NW_019211455.1_6862184 <NA>
## 766 NW_019211455.1_6892496 0/0
## 767 NW_019211455.1_6901533 0/0
## 768 NW_019211455.1_7375407 0/0
## 769 NW_019211455.1_7621250 <NA>
## 770 NW_019211455.1_7760562 0/0
## 771 NW_019211455.1_9623242 0/0
## 772 NW_019211455.1_9624267 0/0
## 773 NW_019211455.1_9986181 <NA>
## 774 NW_019211455.1_10308072 0/0
## 775 NW_019211455.1_10351749 0/0
## 776 NW_019211455.1_10712687 0/0
## 777 NW_019211455.1_10724719 <NA>
## 778 NW_019211455.1_10739248 0/0
## 779 NW_019211455.1_10787412 0/0
## 780 NW_019211455.1_10797398 0/0
## 781 NW_019211455.1_10802538 <NA>
## 782 NW_019211455.1_10815261 <NA>
## 783 NW_019211455.1_10834864 <NA>
## 784 NW_019211455.1_10837041 0/0
## 785 NW_019211455.1_10842021 0/0
## 786 NW_019211455.1_12991126 <NA>
## 787 NW_019211455.1_17149864 0/0
## 788 NW_019211455.1_17185737 0/0
## 789 NW_019211455.1_17252612 0/0
## 790 NW_019211455.1_17563362 0/0
## 791 NW_019211455.1_17585370 0/0
## 792 NW_019211455.1_17595372 0/0
## 793 NW_019211455.1_17653912 0/0
## 794 NW_019211455.1_17673820 0/0
## 795 NW_019211455.1_17678811 0/0
## 796 NW_019211455.1_17740334 <NA>
## 797 NW_019211455.1_17797673 <NA>
## 798 NW_019211455.1_17817512 <NA>
## 799 NW_019211455.1_17879102 0/0
## 800 NW_019211455.1_17882985 0/0
## 801 NW_019211455.1_17896957 0/0
## 802 NW_019211455.1_17906927 <NA>
## 803 NW_019211455.1_17982258 0/0
## 804 NW_019211455.1_18353510 0/0
## 805 NW_019211455.1_22222784 <NA>
## 806 NW_019211455.1_23246093 0/0
## 807 NW_019211455.1_23495473 0/0
## 808 NW_019211455.1_24612052 0/0
## 809 NW_019211455.1_25486809 0/0
## 810 NW_019211455.1_25625408 <NA>
## 811 NW_019211455.1_25882252 0/0
## 812 NW_019211455.1_25905800 <NA>
## 813 NW_019211455.1_25972186 <NA>
## 814 NW_019211455.1_25980384 0/0
## 815 NW_019211455.1_25996836 <NA>
## 816 NW_019211455.1_26000220 0/0
## 817 NW_019211455.1_26031305 <NA>
## 818 NW_019211455.1_26130475 0/0
## 819 NW_019211455.1_26465563 <NA>
## 820 NW_019211455.1_27706793 <NA>
## 821 NW_019211455.1_27775552 <NA>
## 822 NW_019211455.1_27797560 <NA>
## 823 NW_019211455.1_27819122 0/0
## 824 NW_019211455.1_27826549 0/0
## 825 NW_019211455.1_27826921 <NA>
## 826 NW_019211455.1_27835138 <NA>
## 827 NW_019211455.1_27835266 <NA>
## 828 NW_019211455.1_27835677 0/0
## 829 NW_019211455.1_27838161 0/0
## 830 NW_019211455.1_27840881 0/0
## 831 NW_019211455.1_27957912 0/0
## 832 NW_019211455.1_28957975 0/0
## 833 NW_019211455.1_28986332 <NA>
## 834 NW_019211455.1_28986700 0/0
## 835 NW_019211455.1_29111946 0/0
## 836 NW_019211455.1_29120354 <NA>
## 837 NW_019211455.1_29127976 <NA>
## 838 NW_019211455.1_29135396 <NA>
## 839 NW_019211455.1_29140383 0/0
## 840 NW_019211455.1_29142218 <NA>
## 841 NW_019211455.1_29157976 0/0
## 842 NW_019211455.1_29158671 0/0
## 843 NW_019211455.1_29177263 <NA>
## 844 NW_019211455.1_29180050 <NA>
## 845 NW_019211455.1_29185410 0/0
## 846 NW_019211455.1_29186191 0/0
## 847 NW_019211455.1_29200301 0/0
## 848 NW_019211455.1_29202343 0/0
## 849 NW_019211455.1_29260943 <NA>
## 850 NW_019211455.1_29261900 0/0
## 851 NW_019211455.1_29353407 <NA>
## 852 NW_019211455.1_29362379 <NA>
## 853 NW_019211455.1_29387221 0/0
## 854 NW_019211455.1_30201474 0/0
## 855 NW_019211455.1_30227266 0/0
## 856 NW_019211455.1_30231935 0/0
## 857 NW_019211455.1_30738434 0/0
## 858 NW_019211455.1_31207747 <NA>
## 859 NW_019211455.1_32296843 0/0
## 860 NW_019211455.1_32337677 0/0
## 861 NW_019211455.1_32339825 0/0
## 862 NW_019211455.1_32339846 0/0
## 863 NW_019211455.1_32343219 <NA>
## 864 NW_019211455.1_32346622 <NA>
## 865 NW_019211455.1_32356611 <NA>
## 866 NW_019211455.1_32361819 <NA>
## 867 NW_019211455.1_32363967 <NA>
## 868 NW_019211455.1_32369078 <NA>
## 869 NW_019211455.1_32376932 <NA>
## 870 NW_019211455.1_32388092 <NA>
## 871 NW_019211455.1_32390675 <NA>
## 872 NW_019211455.1_32418640 0/0
## 873 NW_019211455.1_32419186 <NA>
## 874 NW_019211455.1_32426852 <NA>
## 875 NW_019211455.1_32430044 0/0
## 876 NW_019211455.1_32434705 0/0
## 877 NW_019211455.1_32441005 0/0
## 878 NW_019211455.1_32445017 0/0
## 879 NW_019211455.1_32468374 0/0
## 880 NW_019211455.1_32469897 0/0
## 881 NW_019211455.1_32470012 0/0
## 882 NW_019211455.1_32474203 <NA>
## 883 NW_019211455.1_32493142 <NA>
## 884 NW_019211455.1_32496172 <NA>
## 885 NW_019211455.1_32497410 0/0
## 886 NW_019211455.1_32507379 0/0
## 887 NW_019211455.1_32512867 0/0
## 888 NW_019211455.1_32515507 0/0
## 889 NW_019211455.1_32522896 0/0
## 890 NW_019211455.1_32523059 <NA>
## 891 NW_019211455.1_32553392 0/0
## 892 NW_019211455.1_32881736 0/0
## 893 NW_019211455.1_32894397 <NA>
## 894 NW_019211455.1_32894669 0/0
## 895 NW_019211455.1_32938799 0/0
## 896 NW_019211455.1_32952007 0/0
## 897 NW_019211455.1_32963757 0/0
## 898 NW_019211455.1_32970140 <NA>
## 899 NW_019211455.1_32971481 <NA>
## 900 NW_019211455.1_32971605 0/0
## 901 NW_019211455.1_32981224 <NA>
## 902 NW_019211455.1_32984952 0/0
## 903 NW_019211455.1_32985410 <NA>
## 904 NW_019211455.1_32994239 <NA>
## 905 NW_019211455.1_33027539 0/0
## 906 NW_019211455.1_33054097 0/0
## 907 NW_019211455.1_33057384 0/0
## 908 NW_019211455.1_33069852 <NA>
## 909 NW_019211455.1_33077797 0/0
## 910 NW_019211455.1_33081788 <NA>
## 911 NW_019211455.1_33081937 <NA>
## 912 NW_019211455.1_33084526 <NA>
## 913 NW_019211455.1_33085106 0/0
## 914 NW_019211455.1_33087714 <NA>
## 915 NW_019211455.1_33092936 <NA>
## 916 NW_019211455.1_33102611 <NA>
## 917 NW_019211455.1_33114831 0/0
## 918 NW_019211455.1_33123530 0/0
## 919 NW_019211455.1_33126806 0/0
## 920 NW_019211455.1_33130026 <NA>
## 921 NW_019211455.1_33200929 0/0
## 922 NW_019211455.1_33226353 0/0
## 923 NW_019211455.1_33262612 0/0
## 924 NW_019211455.1_33292229 0/0
## 925 NW_019211455.1_33293808 0/0
## 926 NW_019211455.1_33299254 <NA>
## 927 NW_019211455.1_33300015 0/0
## 928 NW_019211455.1_33302561 0/0
## 929 NW_019211455.1_33313586 0/0
## 930 NW_019211455.1_33319912 <NA>
## 931 NW_019211455.1_33320213 0/0
## 932 NW_019211455.1_33321615 <NA>
## 933 NW_019211455.1_33350023 0/0
## 934 NW_019211455.1_33352594 0/0
## 935 NW_019211455.1_33359991 <NA>
## 936 NW_019211455.1_33362669 <NA>
## 937 NW_019211455.1_33374791 0/0
## 938 NW_019211455.1_33377643 <NA>
## 939 NW_019211455.1_33403086 <NA>
## 940 NW_019211455.1_33423204 0/0
## 941 NW_019211455.1_33438484 <NA>
## 942 NW_019211455.1_33439183 0/0
## 943 NW_019211455.1_33466919 0/0
## 944 NW_019211455.1_33476515 0/0
## 945 NW_019211455.1_33486164 <NA>
## 946 NW_019211455.1_33532935 0/0
## 947 NW_019211455.1_33556976 0/0
## 948 NW_019211455.1_33603645 0/0
## 949 NW_019211455.1_33630706 0/0
## 950 NW_019211455.1_33635799 0/0
## 951 NW_019211455.1_33745725 <NA>
## 952 NW_019211455.1_33751908 0/0
## 953 NW_019211455.1_33752085 0/0
## 954 NW_019211455.1_33843457 0/0
## 955 NW_019211455.1_33885647 <NA>
## 956 NW_019211455.1_33954535 0/0
## 957 NW_019211455.1_33967120 0/0
## 958 NW_019211455.1_33970776 <NA>
## 959 NW_019211455.1_43992022 0/0
## 960 NW_019211455.1_44010288 0/0
## 961 NW_019211455.1_44018236 0/0
## 962 NW_019211455.1_44080482 <NA>
## 963 NW_019211455.1_44094732 <NA>
## 964 NW_019211455.1_44115058 0/0
## 965 NW_019211455.1_44126474 <NA>
## 966 NW_019211455.1_44150586 <NA>
## 967 NW_019211455.1_44150866 <NA>
## 968 NW_019211455.1_44168728 0/0
## 969 NW_019211455.1_44169811 0/0
## 970 NW_019211455.1_44273021 0/0
## 971 NW_019211455.1_44345826 0/0
## 972 NW_019211455.1_44379285 0/0
## 973 NW_019211455.1_44425219 <NA>
## 974 NW_019211455.1_45942338 0/0
## 975 NW_019211455.1_45947753 0/0
## 976 NW_019211455.1_45952601 <NA>
## 977 NW_019211455.1_45999743 0/0
## 978 NW_019211455.1_46023005 <NA>
## 979 NW_019211455.1_46056725 0/0
## 980 NW_019211455.1_46143785 0/0
## 981 NW_019211455.1_46173048 0/0
## 982 NW_019211455.1_46212761 0/0
## 983 NW_019211455.1_46230879 <NA>
## 984 NW_019211455.1_46232773 <NA>
## 985 NW_019211455.1_46280279 0/0
## 986 NW_019211455.1_46352577 0/0
## 987 NW_019211455.1_46358307 0/0
## 988 NW_019211455.1_46398917 0/0
## 989 NW_019211455.1_46402800 <NA>
## 990 NW_019211455.1_46403999 <NA>
## 991 NW_019211455.1_46422811 0/1
## 992 NW_019211455.1_46422983 0/0
## 993 NW_019211455.1_46438010 0/0
## 994 NW_019211455.1_46463185 <NA>
## 995 NW_019211455.1_46482603 0/0
## 996 NW_019211455.1_46490910 <NA>
## 997 NW_019211455.1_46540183 <NA>
## 998 NW_019211455.1_46569030 <NA>
## 999 NW_019211455.1_46628280 <NA>
## 1000 NW_019211455.1_46635397 0/0
## 1001 NW_019211455.1_46645131 0/0
## 1002 NW_019211455.1_46648663 0/0
## 1003 NW_019211455.1_46657450 0/0
## 1004 NW_019211455.1_46715628 0/0
## 1005 NW_019211455.1_46758223 0/0
## 1006 NW_019211455.1_46758853 <NA>
## 1007 NW_019211455.1_46835426 0/0
## 1008 NW_019211455.1_46840769 <NA>
## 1009 NW_019211455.1_46856747 0/0
## 1010 NW_019211455.1_46856806 0/0
## 1011 NW_019211455.1_46884511 0/0
## 1012 NW_019211455.1_46898007 0/0
## 1013 NW_019211455.1_46914680 <NA>
## 1014 NW_019211455.1_46918514 <NA>
## 1015 NW_019211455.1_46953584 0/0
## 1016 NW_019211455.1_46984811 0/0
## 1017 NW_019211455.1_47022470 0/0
## 1018 NW_019211455.1_47042771 0/0
## 1019 NW_019211455.1_47077891 <NA>
## 1020 NW_019211455.1_47144157 <NA>
## 1021 NW_019211455.1_47179282 <NA>
## 1022 NW_019211455.1_47179293 <NA>
## 1023 NW_019211455.1_47179374 <NA>
## 1024 NW_019211455.1_47193481 0/0
## 1025 NW_019211455.1_47202010 0/0
## 1026 NW_019211455.1_47202462 0/0
## 1027 NW_019211455.1_47212866 <NA>
## 1028 NW_019211455.1_47237571 0/0
## 1029 NW_019211455.1_47248373 <NA>
## 1030 NW_019211455.1_47401099 0/0
## 1031 NW_019211455.1_47425962 <NA>
## 1032 NW_019211455.1_47452310 0/0
## 1033 NW_019211455.1_47701640 <NA>
## 1034 NW_019211455.1_48198008 0/0
## 1035 NW_019211455.1_48208542 0/0
## 1036 NW_019211455.1_48245617 <NA>
## 1037 NW_019211455.1_48258947 <NA>
## 1038 NW_019211455.1_48258992 <NA>
## 1039 NW_019211455.1_51129038 <NA>
## 1040 NW_019211455.1_51307592 0/0
## 1041 NW_019211455.1_52825989 0/0
## 1042 NW_019211455.1_52846975 <NA>
## 1043 NW_019211455.1_52879066 0/0
## 1044 NW_019211455.1_52882124 <NA>
## 1045 NW_019211455.1_52886915 0/0
## 1046 NW_019211455.1_52890367 0/0
## 1047 NW_019211455.1_52893316 <NA>
## 1048 NW_019211455.1_52911389 0/0
## 1049 NW_019211455.1_52937516 0/0
## 1050 NW_019211455.1_52943625 <NA>
## 1051 NW_019211455.1_52943738 0/0
## 1052 NW_019211455.1_52947377 <NA>
## 1053 NW_019211455.1_52995750 0/0
## 1054 NW_019211455.1_53540635 0/0
## 1055 NW_019211455.1_54380106 <NA>
## 1056 NW_019211455.1_54381574 <NA>
## 1057 NW_019211455.1_54386280 0/0
## 1058 NW_019211455.1_54389650 <NA>
## 1059 NW_019211455.1_54395668 <NA>
## 1060 NW_019211455.1_54402567 0/0
## 1061 NW_019211455.1_54411340 <NA>
## 1062 NW_019211455.1_54421016 <NA>
## 1063 NW_019211455.1_54491786 0/0
## 1064 NW_019211455.1_54491868 <NA>
## 1065 NW_019211455.1_54493357 <NA>
## 1066 NW_019211455.1_54498426 <NA>
## 1067 NW_019211455.1_54527577 0/0
## 1068 NW_019211455.1_54553702 0/0
## 1069 NW_019211455.1_54576375 <NA>
## 1070 NW_019211455.1_54607362 <NA>
## 1071 NW_019211455.1_54624652 0/0
## 1072 NW_019211455.1_54644131 0/0
## 1073 NW_019211455.1_54651895 0/0
## 1074 NW_019211455.1_54657681 0/0
## 1075 NW_019211455.1_54669062 <NA>
## 1076 NW_019211455.1_54669579 0/0
## 1077 NW_019211455.1_54687546 0/0
## 1078 NW_019211455.1_54717899 0/0
## 1079 NW_019211455.1_54724953 0/0
## 1080 NW_019211455.1_54732453 0/0
## 1081 NW_019211455.1_54737739 0/0
## 1082 NW_019211455.1_54737924 0/0
## 1083 NW_019211455.1_54758026 <NA>
## 1084 NW_019211455.1_54759583 <NA>
## 1085 NW_019211455.1_54762916 <NA>
## 1086 NW_019211455.1_54782102 0/0
## 1087 NW_019211455.1_54783615 0/0
## 1088 NW_019211455.1_54863497 0/0
## 1089 NW_019211455.1_54879180 0/0
## 1090 NW_019211455.1_54883036 0/0
## 1091 NW_019211455.1_54965397 0/0
## 1092 NW_019211455.1_54979264 0/0
## 1093 NW_019211455.1_54983705 0/0
## 1094 NW_019211455.1_55648597 0/0
## 1095 NW_019211455.1_55689498 0/0
## 1096 NW_019211455.1_55705842 0/0
## 1097 NW_019211455.1_55776197 <NA>
## 1098 NW_019211455.1_55802770 <NA>
## 1099 NW_019211455.1_55855535 0/0
## 1100 NW_019211455.1_55857271 0/0
## 1101 NW_019211455.1_55862495 <NA>
## 1102 NW_019211455.1_55873586 0/0
## 1103 NW_019211455.1_55895887 <NA>
## 1104 NW_019211455.1_55899517 <NA>
## 1105 NW_019211455.1_55904098 <NA>
## 1106 NW_019211455.1_55941737 <NA>
## 1107 NW_019211455.1_55958958 0/0
## 1108 NW_019211455.1_55998586 0/0
## 1109 NW_019211455.1_56308420 <NA>
## 1110 NW_019211455.1_56962613 <NA>
## 1111 NW_019211456.1_852213 <NA>
## 1112 NW_019211456.1_862934 0/0
## 1113 NW_019211456.1_1514091 0/0
## 1114 NW_019211456.1_1557845 0/0
## 1115 NW_019211456.1_1557913 0/0
## 1116 NW_019211456.1_1575938 <NA>
## 1117 NW_019211456.1_1576213 0/0
## 1118 NW_019211456.1_1581618 <NA>
## 1119 NW_019211456.1_1583683 0/0
## 1120 NW_019211456.1_1634695 0/0
## 1121 NW_019211456.1_1657692 <NA>
## 1122 NW_019211456.1_1662405 0/0
## 1123 NW_019211456.1_1671159 <NA>
## 1124 NW_019211456.1_1682880 0/0
## 1125 NW_019211456.1_3857748 0/0
## 1126 NW_019211456.1_3886580 0/0
## 1127 NW_019211456.1_3906070 <NA>
## 1128 NW_019211456.1_3932173 0/0
## 1129 NW_019211456.1_3938341 0/0
## 1130 NW_019211456.1_3941809 0/0
## 1131 NW_019211456.1_3960878 0/0
## 1132 NW_019211456.1_3983269 0/0
## 1133 NW_019211456.1_3996604 <NA>
## 1134 NW_019211456.1_4007274 0/0
## 1135 NW_019211456.1_4010857 0/0
## 1136 NW_019211456.1_4035995 0/0
## 1137 NW_019211456.1_4049358 <NA>
## 1138 NW_019211456.1_4061699 0/0
## 1139 NW_019211456.1_4067032 0/0
## 1140 NW_019211456.1_4092318 0/0
## 1141 NW_019211456.1_4097747 0/0
## 1142 NW_019211456.1_4098080 0/0
## 1143 NW_019211456.1_4123215 0/0
## 1144 NW_019211456.1_4138181 <NA>
## 1145 NW_019211456.1_4303393 0/0
## 1146 NW_019211456.1_5959685 <NA>
## 1147 NW_019211456.1_8208869 0/0
## 1148 NW_019211456.1_9135241 <NA>
## 1149 NW_019211456.1_9568584 <NA>
## 1150 NW_019211456.1_9614543 <NA>
## 1151 NW_019211456.1_9615304 0/0
## 1152 NW_019211456.1_10707372 0/1
## 1153 NW_019211456.1_11003117 0/0
## 1154 NW_019211456.1_11028134 0/0
## 1155 NW_019211456.1_12647375 <NA>
## 1156 NW_019211456.1_12878824 0/0
## 1157 NW_019211456.1_12880280 0/1
## 1158 NW_019211456.1_12892547 0/0
## 1159 NW_019211456.1_12893809 0/0
## 1160 NW_019211456.1_13032008 0/1
## 1161 NW_019211456.1_13066662 <NA>
## 1162 NW_019211456.1_13659623 <NA>
## 1163 NW_019211456.1_13681742 0/0
## 1164 NW_019211456.1_15021677 0/0
## 1165 NW_019211456.1_15035896 <NA>
## 1166 NW_019211456.1_15048085 <NA>
## 1167 NW_019211456.1_15048146 <NA>
## 1168 NW_019211456.1_15053283 <NA>
## 1169 NW_019211456.1_15054416 0/0
## 1170 NW_019211456.1_15055335 0/0
## 1171 NW_019211456.1_15058545 <NA>
## 1172 NW_019211456.1_15061215 <NA>
## 1173 NW_019211456.1_15093465 0/0
## 1174 NW_019211456.1_15103042 <NA>
## 1175 NW_019211456.1_15105080 <NA>
## 1176 NW_019211456.1_15106684 0/0
## 1177 NW_019211456.1_15139477 <NA>
## 1178 NW_019211456.1_15141612 <NA>
## 1179 NW_019211456.1_15165693 0/0
## 1180 NW_019211456.1_15201102 <NA>
## 1181 NW_019211456.1_15207348 <NA>
## 1182 NW_019211456.1_15224287 0/1
## 1183 NW_019211456.1_15237881 0/0
## 1184 NW_019211456.1_15243567 0/0
## 1185 NW_019211456.1_15243642 <NA>
## 1186 NW_019211456.1_15261402 0/0
## 1187 NW_019211456.1_16320717 <NA>
## 1188 NW_019211456.1_16605371 0/0
## 1189 NW_019211456.1_16625728 0/0
## 1190 NW_019211456.1_17066297 0/0
## 1191 NW_019211456.1_17069158 <NA>
## 1192 NW_019211456.1_17085080 <NA>
## 1193 NW_019211456.1_17089494 0/0
## 1194 NW_019211456.1_17093759 <NA>
## 1195 NW_019211456.1_17094459 0/0
## 1196 NW_019211456.1_17107652 <NA>
## 1197 NW_019211456.1_17196021 0/0
## 1198 NW_019211456.1_17206444 0/0
## 1199 NW_019211456.1_17206854 0/0
## 1200 NW_019211456.1_17216707 <NA>
## 1201 NW_019211456.1_17216735 <NA>
## 1202 NW_019211456.1_17230296 0/0
## 1203 NW_019211456.1_17238107 0/0
## 1204 NW_019211456.1_17271127 <NA>
## 1205 NW_019211456.1_17272472 0/0
## 1206 NW_019211456.1_17303572 0/0
## 1207 NW_019211456.1_17303779 0/0
## 1208 NW_019211456.1_18084169 0/0
## 1209 NW_019211456.1_18088405 0/0
## 1210 NW_019211456.1_18091807 0/0
## 1211 NW_019211456.1_18094347 <NA>
## 1212 NW_019211456.1_18097440 0/0
## 1213 NW_019211456.1_18106795 <NA>
## 1214 NW_019211456.1_18109673 0/0
## 1215 NW_019211456.1_18110081 <NA>
## 1216 NW_019211456.1_18111667 <NA>
## 1217 NW_019211456.1_18117542 <NA>
## 1218 NW_019211456.1_18129541 0/0
## 1219 NW_019211456.1_18141664 0/0
## 1220 NW_019211456.1_18142270 <NA>
## 1221 NW_019211456.1_18168105 0/0
## 1222 NW_019211456.1_18168755 <NA>
## 1223 NW_019211456.1_18177475 0/0
## 1224 NW_019211456.1_18192142 0/0
## 1225 NW_019211456.1_18207954 0/0
## 1226 NW_019211456.1_18242958 0/0
## 1227 NW_019211456.1_18254612 0/0
## 1228 NW_019211456.1_18257798 <NA>
## 1229 NW_019211456.1_18268466 0/0
## 1230 NW_019211456.1_18284274 <NA>
## 1231 NW_019211456.1_18288521 0/0
## 1232 NW_019211456.1_18348868 <NA>
## 1233 NW_019211456.1_18349330 <NA>
## 1234 NW_019211456.1_18352718 0/0
## 1235 NW_019211456.1_18357740 0/0
## 1236 NW_019211456.1_18360835 <NA>
## 1237 NW_019211456.1_18377185 <NA>
## 1238 NW_019211456.1_18391231 0/0
## 1239 NW_019211456.1_18405490 0/0
## 1240 NW_019211456.1_18406228 0/0
## 1241 NW_019211456.1_18406565 0/0
## 1242 NW_019211456.1_18409154 0/0
## 1243 NW_019211456.1_18435890 <NA>
## 1244 NW_019211456.1_18452545 <NA>
## 1245 NW_019211456.1_18456498 0/0
## 1246 NW_019211456.1_18456594 0/0
## 1247 NW_019211456.1_18474923 <NA>
## 1248 NW_019211456.1_18476816 <NA>
## 1249 NW_019211456.1_18481702 0/0
## 1250 NW_019211456.1_18481713 0/0
## 1251 NW_019211456.1_18485302 0/0
## 1252 NW_019211456.1_18502236 0/0
## 1253 NW_019211456.1_18508056 0/0
## 1254 NW_019211456.1_18524257 0/0
## 1255 NW_019211456.1_18526248 0/0
## 1256 NW_019211456.1_18527486 <NA>
## 1257 NW_019211456.1_18541153 <NA>
## 1258 NW_019211456.1_18545531 <NA>
## 1259 NW_019211456.1_18546450 <NA>
## 1260 NW_019211456.1_18577461 0/1
## 1261 NW_019211456.1_18577756 0/0
## 1262 NW_019211456.1_18577759 0/0
## 1263 NW_019211456.1_19439475 0/0
## 1264 NW_019211456.1_20181164 0/0
## 1265 NW_019211456.1_20221719 <NA>
## 1266 NW_019211456.1_20225265 0/0
## 1267 NW_019211456.1_20245567 0/0
## 1268 NW_019211456.1_20262180 <NA>
## 1269 NW_019211456.1_20635312 0/0
## 1270 NW_019211456.1_20723040 0/0
## 1271 NW_019211456.1_20731173 0/0
## 1272 NW_019211456.1_20738096 <NA>
## 1273 NW_019211456.1_20738622 0/0
## 1274 NW_019211456.1_20751785 0/0
## 1275 NW_019211456.1_20907700 0/0
## 1276 NW_019211456.1_20925775 0/0
## 1277 NW_019211456.1_20931498 0/0
## 1278 NW_019211456.1_21507157 0/0
## 1279 NW_019211456.1_21539353 0/0
## 1280 NW_019211456.1_21544242 0/0
## 1281 NW_019211456.1_21548005 0/0
## 1282 NW_019211456.1_21578834 <NA>
## 1283 NW_019211456.1_21582779 <NA>
## 1284 NW_019211456.1_21590017 <NA>
## 1285 NW_019211456.1_21623040 <NA>
## 1286 NW_019211456.1_21623240 <NA>
## 1287 NW_019211456.1_22371006 <NA>
## 1288 NW_019211456.1_22845656 <NA>
## 1289 NW_019211456.1_22921703 0/0
## 1290 NW_019211456.1_22928635 <NA>
## 1291 NW_019211456.1_22943032 0/0
## 1292 NW_019211456.1_22962417 0/0
## 1293 NW_019211456.1_22980371 <NA>
## 1294 NW_019211456.1_23001515 0/0
## 1295 NW_019211456.1_23034338 0/0
## 1296 NW_019211456.1_23043887 0/0
## 1297 NW_019211456.1_23062534 <NA>
## 1298 NW_019211456.1_23089635 0/1
## 1299 NW_019211456.1_23116200 <NA>
## 1300 NW_019211456.1_23126151 <NA>
## 1301 NW_019211456.1_23140667 <NA>
## 1302 NW_019211456.1_23185651 0/0
## 1303 NW_019211456.1_23231210 0/0
## 1304 NW_019211456.1_23240943 <NA>
## 1305 NW_019211456.1_23247083 0/0
## 1306 NW_019211456.1_23283676 <NA>
## 1307 NW_019211456.1_23283737 <NA>
## 1308 NW_019211456.1_23284592 0/0
## 1309 NW_019211456.1_23310357 0/0
## 1310 NW_019211456.1_23387640 0/0
## 1311 NW_019211456.1_23388875 <NA>
## 1312 NW_019211456.1_23389214 0/0
## 1313 NW_019211456.1_23392562 0/0
## 1314 NW_019211456.1_23392934 0/0
## 1315 NW_019211456.1_23392953 0/0
## 1316 NW_019211456.1_23393223 0/0
## 1317 NW_019211456.1_23393805 0/0
## 1318 NW_019211456.1_23397535 <NA>
## 1319 NW_019211456.1_23398890 <NA>
## 1320 NW_019211456.1_23413176 <NA>
## 1321 NW_019211456.1_23413427 0/0
## 1322 NW_019211456.1_23413501 0/0
## 1323 NW_019211456.1_23432825 0/0
## 1324 NW_019211456.1_23447346 0/0
## 1325 NW_019211456.1_23503305 <NA>
## 1326 NW_019211456.1_23548980 0/0
## 1327 NW_019211456.1_23556990 <NA>
## 1328 NW_019211456.1_23569363 0/0
## 1329 NW_019211456.1_24801804 0/0
## 1330 NW_019211456.1_24825156 <NA>
## 1331 NW_019211456.1_24828112 0/0
## 1332 NW_019211456.1_24853255 0/0
## 1333 NW_019211456.1_24862822 0/0
## 1334 NW_019211456.1_24870581 <NA>
## 1335 NW_019211456.1_24876402 <NA>
## 1336 NW_019211456.1_24882138 0/0
## 1337 NW_019211456.1_24883429 0/0
## 1338 NW_019211456.1_24886256 0/0
## 1339 NW_019211456.1_24886443 0/0
## 1340 NW_019211456.1_24888644 0/0
## 1341 NW_019211456.1_24889364 0/0
## 1342 NW_019211456.1_24892386 0/0
## 1343 NW_019211456.1_24896884 0/0
## 1344 NW_019211456.1_24896930 0/0
## 1345 NW_019211456.1_24900777 <NA>
## 1346 NW_019211456.1_25978208 <NA>
## 1347 NW_019211456.1_25978248 <NA>
## 1348 NW_019211456.1_27009527 0/1
## 1349 NW_019211456.1_27186402 0/0
## 1350 NW_019211456.1_27219976 <NA>
## 1351 NW_019211456.1_28608466 0/0
## 1352 NW_019211456.1_28608768 0/0
## 1353 NW_019211456.1_28608973 0/0
## 1354 NW_019211456.1_28623307 0/0
## 1355 NW_019211456.1_28638661 <NA>
## 1356 NW_019211456.1_28679024 <NA>
## 1357 NW_019211456.1_28681580 <NA>
## 1358 NW_019211456.1_28684214 <NA>
## 1359 NW_019211456.1_28698184 0/0
## 1360 NW_019211456.1_28702632 <NA>
## 1361 NW_019211456.1_28705307 0/0
## 1362 NW_019211456.1_28709211 <NA>
## 1363 NW_019211456.1_28709934 0/0
## 1364 NW_019211456.1_28713209 <NA>
## 1365 NW_019211456.1_28738435 <NA>
## 1366 NW_019211456.1_28746482 <NA>
## 1367 NW_019211456.1_28756949 <NA>
## 1368 NW_019211456.1_28760675 <NA>
## 1369 NW_019211456.1_28771652 <NA>
## 1370 NW_019211456.1_28774982 <NA>
## 1371 NW_019211456.1_28775099 <NA>
## 1372 NW_019211456.1_28775116 <NA>
## 1373 NW_019211456.1_28776037 <NA>
## 1374 NW_019211456.1_28776666 0/0
## 1375 NW_019211456.1_28783428 0/0
## 1376 NW_019211456.1_28810341 0/0
## 1377 NW_019211456.1_28831177 <NA>
## 1378 NW_019211456.1_28835793 <NA>
## 1379 NW_019211456.1_29148731 0/0
## 1380 NW_019211456.1_29185663 <NA>
## 1381 NW_019211456.1_29851529 0/0
## 1382 NW_019211456.1_32158684 <NA>
## 1383 NW_019211456.1_40655890 0/0
## 1384 NW_019211456.1_41613285 0/0
## 1385 NW_019211456.1_42335654 0/0
## 1386 NW_019211456.1_42343263 <NA>
## 1387 NW_019211456.1_42347209 <NA>
## 1388 NW_019211456.1_42347330 0/0
## 1389 NW_019211456.1_42349826 0/0
## 1390 NW_019211456.1_42349837 0/0
## 1391 NW_019211456.1_42350229 0/0
## 1392 NW_019211456.1_42367184 0/0
## 1393 NW_019211456.1_42377338 <NA>
## 1394 NW_019211456.1_42380344 0/0
## 1395 NW_019211456.1_42382732 <NA>
## 1396 NW_019211456.1_42392011 <NA>
## 1397 NW_019211456.1_42400142 0/0
## 1398 NW_019211456.1_42412002 0/0
## 1399 NW_019211456.1_42412302 <NA>
## 1400 NW_019211456.1_42412364 <NA>
## 1401 NW_019211456.1_42412490 0/0
## 1402 NW_019211456.1_42413110 <NA>
## 1403 NW_019211456.1_42414609 <NA>
## 1404 NW_019211456.1_42415812 0/0
## 1405 NW_019211456.1_42422591 0/0
## 1406 NW_019211456.1_42476727 0/0
## 1407 NW_019211456.1_42482262 0/0
## 1408 NW_019211456.1_42485133 0/0
## 1409 NW_019211456.1_42485611 0/0
## 1410 NW_019211456.1_42515426 <NA>
## 1411 NW_019211456.1_42518980 0/0
## 1412 NW_019211456.1_42520527 <NA>
## 1413 NW_019211456.1_42522460 <NA>
## 1414 NW_019211456.1_42526371 <NA>
## 1415 NW_019211456.1_42526580 <NA>
## 1416 NW_019211456.1_42528120 0/0
## 1417 NW_019211456.1_42537823 <NA>
## 1418 NW_019211456.1_42541979 <NA>
## 1419 NW_019211456.1_42993189 <NA>
## 1420 NW_019211456.1_42993307 <NA>
## 1421 NW_019211456.1_45175731 0/0
## 1422 NW_019211456.1_45184752 <NA>
## 1423 NW_019211456.1_45185303 <NA>
## 1424 NW_019211456.1_45185562 0/0
## 1425 NW_019211456.1_45192221 0/0
## 1426 NW_019211456.1_45202115 0/0
## 1427 NW_019211456.1_45208993 0/0
## 1428 NW_019211456.1_45209008 0/0
## 1429 NW_019211456.1_45491709 <NA>
## 1430 NW_019211456.1_47124660 0/0
## 1431 NW_019211456.1_47125448 <NA>
## 1432 NW_019211456.1_47141944 0/0
## 1433 NW_019211456.1_47155456 0/0
## 1434 NW_019211456.1_47165647 0/0
## 1435 NW_019211456.1_47186039 0/0
## 1436 NW_019211456.1_47263170 0/0
## 1437 NW_019211456.1_47282197 0/0
## 1438 NW_019211456.1_47294803 0/0
## 1439 NW_019211456.1_47301651 0/0
## 1440 NW_019211456.1_47328092 0/0
## 1441 NW_019211456.1_47340652 0/0
## 1442 NW_019211456.1_47343892 <NA>
## 1443 NW_019211456.1_47346605 <NA>
## 1444 NW_019211456.1_47349062 <NA>
## 1445 NW_019211456.1_47370081 <NA>
## 1446 NW_019211456.1_47375562 <NA>
## 1447 NW_019211456.1_47409520 0/0
## 1448 NW_019211456.1_47421842 0/0
## 1449 NW_019211456.1_47469223 0/0
## 1450 NW_019211456.1_47475278 0/0
## 1451 NW_019211456.1_47483686 <NA>
## 1452 NW_019211456.1_47483847 0/0
## 1453 NW_019211456.1_47488539 0/0
## 1454 NW_019211456.1_47508413 0/0
## 1455 NW_019211456.1_47514026 <NA>
## 1456 NW_019211456.1_47514172 0/0
## 1457 NW_019211456.1_47518355 0/0
## 1458 NW_019211456.1_47519263 0/0
## 1459 NW_019211456.1_47520838 <NA>
## 1460 NW_019211456.1_47535679 0/0
## 1461 NW_019211456.1_47545095 0/0
## 1462 NW_019211456.1_47552379 0/0
## 1463 NW_019211456.1_47574407 0/1
## 1464 NW_019211456.1_47588360 0/0
## 1465 NW_019211456.1_47591931 0/0
## 1466 NW_019211456.1_47604128 0/0
## 1467 NW_019211456.1_47627929 0/0
## 1468 NW_019211456.1_47629543 0/0
## 1469 NW_019211456.1_47638249 0/0
## 1470 NW_019211456.1_47644566 0/0
## 1471 NW_019211456.1_47651302 0/0
## 1472 NW_019211456.1_47651579 0/0
## 1473 NW_019211456.1_47651788 0/0
## 1474 NW_019211456.1_47657124 0/0
## 1475 NW_019211456.1_47659056 <NA>
## 1476 NW_019211456.1_47663936 0/0
## 1477 NW_019211456.1_47666195 <NA>
## 1478 NW_019211456.1_47692273 0/0
## 1479 NW_019211456.1_47731162 0/0
## 1480 NW_019211456.1_47754687 <NA>
## 1481 NW_019211456.1_47761819 0/0
## 1482 NW_019211456.1_47801183 0/0
## 1483 NW_019211456.1_47804982 <NA>
## 1484 NW_019211456.1_47805594 <NA>
## 1485 NW_019211456.1_47833970 <NA>
## 1486 NW_019211456.1_47833974 <NA>
## 1487 NW_019211456.1_47861589 0/0
## 1488 NW_019211456.1_47869361 0/0
## 1489 NW_019211456.1_47874421 0/0
## 1490 NW_019211456.1_47875404 <NA>
## 1491 NW_019211456.1_47877435 <NA>
## 1492 NW_019211456.1_47886430 0/0
## 1493 NW_019211456.1_47893392 <NA>
## 1494 NW_019211456.1_47895141 <NA>
## 1495 NW_019211456.1_47896588 <NA>
## 1496 NW_019211456.1_47900582 <NA>
## 1497 NW_019211456.1_47902566 0/0
## 1498 NW_019211456.1_47908621 <NA>
## 1499 NW_019211456.1_47919263 0/0
## 1500 NW_019211456.1_47920455 <NA>
## 1501 NW_019211456.1_47920705 <NA>
## 1502 NW_019211456.1_47927215 <NA>
## 1503 NW_019211456.1_48371353 0/0
## 1504 NW_019211456.1_48387305 0/0
## 1505 NW_019211456.1_48387575 <NA>
## 1506 NW_019211456.1_49269730 <NA>
## 1507 NW_019211456.1_49279645 0/0
## 1508 NW_019211456.1_49282315 <NA>
## 1509 NW_019211456.1_49286349 <NA>
## 1510 NW_019211456.1_49288167 <NA>
## 1511 NW_019211456.1_49418838 0/0
## 1512 NW_019211456.1_49448825 <NA>
## 1513 NW_019211456.1_49449249 <NA>
## 1514 NW_019211456.1_49505240 0/0
## 1515 NW_019211456.1_49701065 <NA>
## 1516 NW_019211456.1_49758609 0/0
## 1517 NW_019211456.1_49831647 0/0
## 1518 NW_019211456.1_50325913 0/0
## 1519 NW_019211456.1_51947413 <NA>
## 1520 NW_019211456.1_54333703 <NA>
## 1521 NW_019211457.1_133204 0/0
## 1522 NW_019211457.1_1075332 0/0
## 1523 NW_019211457.1_1087768 <NA>
## 1524 NW_019211457.1_1094772 0/0
## 1525 NW_019211457.1_1123705 <NA>
## 1526 NW_019211457.1_3268466 <NA>
## 1527 NW_019211457.1_5390569 <NA>
## 1528 NW_019211457.1_5424203 0/0
## 1529 NW_019211457.1_5435145 <NA>
## 1530 NW_019211457.1_5435387 0/0
## 1531 NW_019211457.1_5439526 <NA>
## 1532 NW_019211457.1_5444660 <NA>
## 1533 NW_019211457.1_5446232 0/0
## 1534 NW_019211457.1_5447333 0/0
## 1535 NW_019211457.1_5461792 <NA>
## 1536 NW_019211457.1_5487400 0/0
## 1537 NW_019211457.1_5529960 <NA>
## 1538 NW_019211457.1_5539256 0/0
## 1539 NW_019211457.1_5540649 0/0
## 1540 NW_019211457.1_5547972 0/0
## 1541 NW_019211457.1_5563980 0/0
## 1542 NW_019211457.1_5605089 <NA>
## 1543 NW_019211457.1_5631514 0/0
## 1544 NW_019211457.1_5633509 <NA>
## 1545 NW_019211457.1_5638034 0/0
## 1546 NW_019211457.1_5690429 <NA>
## 1547 NW_019211457.1_5699077 <NA>
## 1548 NW_019211457.1_5700088 0/0
## 1549 NW_019211457.1_5703761 <NA>
## 1550 NW_019211457.1_5704004 <NA>
## 1551 NW_019211457.1_5704332 <NA>
## 1552 NW_019211457.1_7355710 <NA>
## 1553 NW_019211457.1_7379365 0/0
## 1554 NW_019211457.1_7395096 0/0
## 1555 NW_019211457.1_7422731 <NA>
## 1556 NW_019211457.1_7423898 <NA>
## 1557 NW_019211457.1_7430860 <NA>
## 1558 NW_019211457.1_7439215 0/0
## 1559 NW_019211457.1_7455656 <NA>
## 1560 NW_019211457.1_7455936 0/0
## 1561 NW_019211457.1_7462242 <NA>
## 1562 NW_019211457.1_7462527 0/0
## 1563 NW_019211457.1_7471753 0/0
## 1564 NW_019211457.1_7483869 <NA>
## 1565 NW_019211457.1_7494182 <NA>
## 1566 NW_019211457.1_7494351 <NA>
## 1567 NW_019211457.1_7509547 0/0
## 1568 NW_019211457.1_7547221 <NA>
## 1569 NW_019211457.1_7548136 0/0
## 1570 NW_019211457.1_7567575 <NA>
## 1571 NW_019211457.1_7568075 0/0
## 1572 NW_019211457.1_7570857 <NA>
## 1573 NW_019211457.1_7571443 <NA>
## 1574 NW_019211457.1_7578642 0/0
## 1575 NW_019211457.1_7580661 <NA>
## 1576 NW_019211457.1_7580666 <NA>
## 1577 NW_019211457.1_9628162 0/0
## 1578 NW_019211457.1_11930151 <NA>
## 1579 NW_019211457.1_12787763 0/0
## 1580 NW_019211457.1_13880645 <NA>
## 1581 NW_019211457.1_16383889 <NA>
## 1582 NW_019211457.1_16386363 0/0
## 1583 NW_019211457.1_16387199 <NA>
## 1584 NW_019211457.1_16387259 0/0
## 1585 NW_019211457.1_16394617 0/0
## 1586 NW_019211457.1_16402138 0/0
## 1587 NW_019211457.1_16410263 0/0
## 1588 NW_019211457.1_16419533 <NA>
## 1589 NW_019211457.1_16429444 0/0
## 1590 NW_019211457.1_16432326 <NA>
## 1591 NW_019211457.1_16458017 <NA>
## 1592 NW_019211457.1_16477347 <NA>
## 1593 NW_019211457.1_16484039 0/0
## 1594 NW_019211457.1_16497602 0/0
## 1595 NW_019211457.1_16501696 <NA>
## 1596 NW_019211457.1_16514827 <NA>
## 1597 NW_019211457.1_16701481 <NA>
## 1598 NW_019211457.1_16866355 0/0
## 1599 NW_019211457.1_17203852 0/0
## 1600 NW_019211457.1_17424004 <NA>
## 1601 NW_019211457.1_17424511 <NA>
## 1602 NW_019211457.1_17429169 0/0
## 1603 NW_019211457.1_17436733 <NA>
## 1604 NW_019211457.1_18658856 0/0
## 1605 NW_019211457.1_18688856 0/0
## 1606 NW_019211457.1_18710827 <NA>
## 1607 NW_019211457.1_18732148 0/0
## 1608 NW_019211457.1_18745006 0/0
## 1609 NW_019211457.1_18759511 0/0
## 1610 NW_019211457.1_18777241 <NA>
## 1611 NW_019211457.1_18820754 0/0
## 1612 NW_019211457.1_18834672 <NA>
## 1613 NW_019211457.1_18856501 <NA>
## 1614 NW_019211457.1_18916077 0/0
## 1615 NW_019211457.1_18942418 <NA>
## 1616 NW_019211457.1_18955864 0/0
## 1617 NW_019211457.1_18967159 0/0
## 1618 NW_019211457.1_18976389 0/0
## 1619 NW_019211457.1_19394120 0/0
## 1620 NW_019211457.1_19455436 0/0
## 1621 NW_019211457.1_19481462 0/0
## 1622 NW_019211457.1_19494042 0/0
## 1623 NW_019211457.1_19494694 0/0
## 1624 NW_019211457.1_19497539 0/0
## 1625 NW_019211457.1_19503800 <NA>
## 1626 NW_019211457.1_21646715 <NA>
## 1627 NW_019211457.1_22464706 <NA>
## 1628 NW_019211457.1_23454754 0/0
## 1629 NW_019211457.1_24968168 0/0
## 1630 NW_019211457.1_25588770 0/0
## 1631 NW_019211457.1_25606378 0/0
## 1632 NW_019211457.1_25608425 0/0
## 1633 NW_019211457.1_27909898 <NA>
## 1634 NW_019211457.1_30253563 0/0
## 1635 NW_019211457.1_30257390 <NA>
## 1636 NW_019211457.1_30267423 <NA>
## 1637 NW_019211457.1_30311075 0/0
## 1638 NW_019211457.1_30380288 <NA>
## 1639 NW_019211457.1_30389811 <NA>
## 1640 NW_019211457.1_30432550 0/0
## 1641 NW_019211457.1_30435284 0/0
## 1642 NW_019211457.1_30439161 <NA>
## 1643 NW_019211457.1_30442378 <NA>
## 1644 NW_019211457.1_30444685 0/0
## 1645 NW_019211457.1_34163229 <NA>
## 1646 NW_019211457.1_34519920 0/0
## 1647 NW_019211457.1_34548550 0/0
## 1648 NW_019211457.1_34555605 0/0
## 1649 NW_019211457.1_34573863 <NA>
## 1650 NW_019211457.1_34618145 0/0
## 1651 NW_019211457.1_34683117 0/0
## 1652 NW_019211457.1_34775151 0/0
## 1653 NW_019211457.1_34815741 <NA>
## 1654 NW_019211457.1_34817873 0/0
## 1655 NW_019211457.1_34823714 0/0
## 1656 NW_019211457.1_34828332 0/0
## 1657 NW_019211457.1_34845237 0/0
## 1658 NW_019211457.1_34856510 <NA>
## 1659 NW_019211457.1_34857001 <NA>
## 1660 NW_019211457.1_34862112 0/0
## 1661 NW_019211457.1_34863195 <NA>
## 1662 NW_019211457.1_34864568 0/0
## 1663 NW_019211457.1_37055534 <NA>
## 1664 NW_019211457.1_37061654 0/0
## 1665 NW_019211457.1_37073286 <NA>
## 1666 NW_019211457.1_37089049 <NA>
## 1667 NW_019211457.1_37091343 <NA>
## 1668 NW_019211457.1_38131790 0/0
## 1669 NW_019211457.1_38148303 0/0
## 1670 NW_019211457.1_38164723 0/0
## 1671 NW_019211457.1_38171473 0/0
## 1672 NW_019211457.1_38171511 0/0
## 1673 NW_019211457.1_38365755 0/0
## 1674 NW_019211457.1_38381478 0/0
## 1675 NW_019211457.1_38382033 <NA>
## 1676 NW_019211457.1_38399758 0/0
## 1677 NW_019211457.1_38399807 0/0
## 1678 NW_019211457.1_38410188 0/0
## 1679 NW_019211457.1_38443887 0/0
## 1680 NW_019211457.1_38463182 0/0
## 1681 NW_019211457.1_38654930 <NA>
## 1682 NW_019211457.1_39381290 0/0
## 1683 NW_019211457.1_39616341 0/0
## 1684 NW_019211457.1_39636521 0/0
## 1685 NW_019211457.1_39651693 0/0
## 1686 NW_019211457.1_39696486 <NA>
## 1687 NW_019211457.1_39958419 <NA>
## 1688 NW_019211457.1_39958571 <NA>
## 1689 NW_019211457.1_39996121 0/0
## 1690 NW_019211457.1_40102951 0/0
## 1691 NW_019211457.1_40260277 <NA>
## 1692 NW_019211457.1_40363203 0/0
## 1693 NW_019211457.1_40368698 0/0
## 1694 NW_019211457.1_40430968 <NA>
## 1695 NW_019211457.1_40452705 0/1
## 1696 NW_019211457.1_40475997 0/0
## 1697 NW_019211457.1_40548622 0/0
## 1698 NW_019211457.1_40565372 <NA>
## 1699 NW_019211457.1_40595641 0/0
## 1700 NW_019211457.1_40596634 <NA>
## 1701 NW_019211457.1_40672507 0/0
## 1702 NW_019211457.1_40680237 <NA>
## 1703 NW_019211457.1_40687798 0/0
## 1704 NW_019211457.1_40690400 0/0
## 1705 NW_019211457.1_40706014 0/0
## 1706 NW_019211457.1_40771663 <NA>
## 1707 NW_019211457.1_43018820 <NA>
## 1708 NW_019211457.1_43026455 <NA>
## 1709 NW_019211457.1_43027748 0/0
## 1710 NW_019211457.1_43056390 0/0
## 1711 NW_019211457.1_43062236 0/0
## 1712 NW_019211457.1_43064173 <NA>
## 1713 NW_019211457.1_43075503 <NA>
## 1714 NW_019211457.1_43082306 <NA>
## 1715 NW_019211457.1_43230758 <NA>
## 1716 NW_019211457.1_43404179 <NA>
## 1717 NW_019211457.1_43412111 0/0
## 1718 NW_019211457.1_43419956 0/0
## 1719 NW_019211457.1_43423239 <NA>
## 1720 NW_019211457.1_43423554 0/0
## 1721 NW_019211457.1_43434782 <NA>
## 1722 NW_019211457.1_43441999 0/0
## 1723 NW_019211457.1_43448690 0/0
## 1724 NW_019211457.1_43460214 0/0
## 1725 NW_019211457.1_43477460 0/0
## 1726 NW_019211457.1_43491491 <NA>
## 1727 NW_019211457.1_43494064 0/0
## 1728 NW_019211457.1_43494184 0/0
## 1729 NW_019211457.1_43500273 <NA>
## 1730 NW_019211457.1_43506468 <NA>
## 1731 NW_019211457.1_43513947 0/0
## 1732 NW_019211457.1_43522475 0/0
## 1733 NW_019211457.1_43524965 0/0
## 1734 NW_019211457.1_43560657 0/0
## 1735 NW_019211457.1_43567633 0/0
## 1736 NW_019211457.1_43570115 <NA>
## 1737 NW_019211457.1_43579126 <NA>
## 1738 NW_019211457.1_43599283 0/0
## 1739 NW_019211457.1_43624623 <NA>
## 1740 NW_019211457.1_43628278 <NA>
## 1741 NW_019211457.1_43646659 0/0
## 1742 NW_019211457.1_43664111 <NA>
## 1743 NW_019211457.1_43668917 0/0
## 1744 NW_019211457.1_43672968 <NA>
## 1745 NW_019211457.1_43688371 0/0
## 1746 NW_019211457.1_43691038 0/0
## 1747 NW_019211457.1_43692692 <NA>
## 1748 NW_019211457.1_43699438 0/0
## 1749 NW_019211457.1_43715387 0/0
## 1750 NW_019211457.1_43716304 <NA>
## 1751 NW_019211457.1_43719012 0/0
## 1752 NW_019211457.1_43734279 <NA>
## 1753 NW_019211457.1_43762754 <NA>
## 1754 NW_019211457.1_43765573 0/0
## 1755 NW_019211457.1_43769507 0/0
## 1756 NW_019211457.1_43770686 <NA>
## 1757 NW_019211457.1_43775412 0/0
## 1758 NW_019211457.1_43781953 <NA>
## 1759 NW_019211457.1_43784201 0/0
## 1760 NW_019211457.1_43798996 0/0
## 1761 NW_019211457.1_43799878 <NA>
## 1762 NW_019211457.1_43819800 <NA>
## 1763 NW_019211457.1_43829987 0/0
## 1764 NW_019211457.1_43833419 <NA>
## 1765 NW_019211457.1_43841047 0/0
## 1766 NW_019211457.1_43865442 0/0
## 1767 NW_019211457.1_43866548 0/0
## 1768 NW_019211457.1_43890042 <NA>
## 1769 NW_019211457.1_43899806 <NA>
## 1770 NW_019211457.1_43911156 0/0
## 1771 NW_019211457.1_43932093 <NA>
## 1772 NW_019211457.1_43932155 <NA>
## 1773 NW_019211457.1_43937584 <NA>
## 1774 NW_019211457.1_43946568 <NA>
## 1775 NW_019211457.1_43955715 <NA>
## 1776 NW_019211457.1_43961122 0/0
## 1777 NW_019211457.1_43988546 0/0
## 1778 NW_019211457.1_44025576 <NA>
## 1779 NW_019211457.1_44064629 0/0
## 1780 NW_019211457.1_44095649 0/0
## 1781 NW_019211457.1_44111224 0/0
## 1782 NW_019211457.1_44130335 <NA>
## 1783 NW_019211457.1_44134277 <NA>
## 1784 NW_019211457.1_44141355 0/0
## 1785 NW_019211457.1_44142885 <NA>
## 1786 NW_019211457.1_44167928 0/0
## 1787 NW_019211457.1_44178254 0/0
## 1788 NW_019211457.1_44200136 0/0
## 1789 NW_019211457.1_44235626 0/0
## 1790 NW_019211457.1_44240251 <NA>
## 1791 NW_019211457.1_44255839 0/0
## 1792 NW_019211457.1_44275899 0/0
## 1793 NW_019211457.1_44295414 <NA>
## 1794 NW_019211457.1_44295502 <NA>
## 1795 NW_019211457.1_45023343 <NA>
## 1796 NW_019211457.1_45034118 0/0
## 1797 NW_019211457.1_45035398 <NA>
## 1798 NW_019211457.1_45042444 0/0
## 1799 NW_019211457.1_45048100 0/0
## 1800 NW_019211457.1_45063030 0/0
## 1801 NW_019211457.1_45071546 <NA>
## 1802 NW_019211457.1_45073516 <NA>
## 1803 NW_019211457.1_50318127 0/0
## 1804 NW_019211457.1_50962362 <NA>
## 1805 NW_019211457.1_51443167 0/0
## 1806 NW_019211457.1_51478363 0/0
## 1807 NW_019211457.1_51513257 0/0
## 1808 NW_019211457.1_51546372 <NA>
## 1809 NW_019211457.1_51581192 <NA>
## 1810 NW_019211457.1_51584412 <NA>
## 1811 NW_019211457.1_51588475 0/0
## 1812 NW_019211457.1_51683102 0/0
## 1813 NW_019211457.1_51705140 <NA>
## 1814 NW_019211457.1_51710311 0/0
## 1815 NW_019211457.1_51761327 0/0
## 1816 NW_019211457.1_51763547 0/0
## 1817 NW_019211457.1_51769826 <NA>
## 1818 NW_019211457.1_51775782 0/0
## 1819 NW_019211457.1_51779339 0/0
## 1820 NW_019211457.1_51873098 0/0
## 1821 NW_019211457.1_51944703 0/0
## 1822 NW_019211457.1_51961139 <NA>
## 1823 NW_019211457.1_51961441 <NA>
## 1824 NW_019211457.1_51976089 0/0
## 1825 NW_019211457.1_52045948 0/0
## 1826 NW_019211457.1_52069480 0/0
## 1827 NW_019211457.1_52104315 <NA>
## 1828 NW_019211457.1_52183824 <NA>
## 1829 NW_019211457.1_52184827 0/0
## 1830 NW_019211457.1_52196943 <NA>
## 1831 NW_019211457.1_52249527 <NA>
## 1832 NW_019211457.1_52258187 0/0
## 1833 NW_019211457.1_52302222 <NA>
## 1834 NW_019211457.1_52321663 0/0
## 1835 NW_019211457.1_52353812 0/0
## 1836 NW_019211457.1_52363162 <NA>
## 1837 NW_019211457.1_52365285 0/0
## 1838 NW_019211457.1_52365614 0/0
## 1839 NW_019211457.1_52387359 0/0
## 1840 NW_019211457.1_52402281 <NA>
## 1841 NW_019211457.1_52404236 0/0
## 1842 NW_019211457.1_52417174 <NA>
## 1843 NW_019211457.1_52417275 0/0
## 1844 NW_019211457.1_52417534 0/0
## 1845 NW_019211457.1_52443043 <NA>
## 1846 NW_019211457.1_52454244 <NA>
## 1847 NW_019211457.1_52469489 0/0
## 1848 NW_019211457.1_52478758 0/0
## 1849 NW_019211457.1_52478780 0/0
## 1850 NW_019211457.1_52479560 0/0
## 1851 NW_019211457.1_52479570 0/0
## 1852 NW_019211457.1_52482140 0/0
## 1853 NW_019211457.1_52661003 0/0
## 1854 NW_019211458.1_1298382 <NA>
## 1855 NW_019211458.1_1583182 <NA>
## 1856 NW_019211458.1_2084304 0/0
## 1857 NW_019211458.1_2563007 <NA>
## 1858 NW_019211458.1_2912793 <NA>
## 1859 NW_019211458.1_3486581 0/0
## 1860 NW_019211458.1_3501170 0/0
## 1861 NW_019211458.1_3517571 0/0
## 1862 NW_019211458.1_3538441 0/0
## 1863 NW_019211458.1_3563440 0/0
## 1864 NW_019211458.1_3594149 <NA>
## 1865 NW_019211458.1_3603498 <NA>
## 1866 NW_019211458.1_3615957 <NA>
## 1867 NW_019211458.1_3632231 0/0
## 1868 NW_019211458.1_3634791 0/0
## 1869 NW_019211458.1_3637222 <NA>
## 1870 NW_019211458.1_3638327 <NA>
## 1871 NW_019211458.1_3638631 <NA>
## 1872 NW_019211458.1_3639626 <NA>
## 1873 NW_019211458.1_3640124 <NA>
## 1874 NW_019211458.1_3660510 0/0
## 1875 NW_019211458.1_3674364 <NA>
## 1876 NW_019211458.1_3677885 <NA>
## 1877 NW_019211458.1_3682080 0/0
## 1878 NW_019211458.1_3697615 0/0
## 1879 NW_019211458.1_3705441 0/0
## 1880 NW_019211458.1_3717645 <NA>
## 1881 NW_019211458.1_3761906 <NA>
## 1882 NW_019211458.1_3904036 0/0
## 1883 NW_019211458.1_3915474 <NA>
## 1884 NW_019211458.1_3932901 0/0
## 1885 NW_019211458.1_3958477 <NA>
## 1886 NW_019211458.1_3972721 <NA>
## 1887 NW_019211458.1_4030552 <NA>
## 1888 NW_019211458.1_4073158 <NA>
## 1889 NW_019211458.1_4084498 <NA>
## 1890 NW_019211458.1_4147248 0/0
## 1891 NW_019211458.1_4161593 0/0
## 1892 NW_019211458.1_4204532 <NA>
## 1893 NW_019211458.1_4237489 <NA>
## 1894 NW_019211458.1_4259639 0/0
## 1895 NW_019211458.1_4265158 <NA>
## 1896 NW_019211458.1_4321607 0/0
## 1897 NW_019211458.1_4328941 0/0
## 1898 NW_019211458.1_4394285 <NA>
## 1899 NW_019211458.1_4503041 <NA>
## 1900 NW_019211458.1_4537007 0/0
## 1901 NW_019211458.1_4561542 <NA>
## 1902 NW_019211458.1_4562049 <NA>
## 1903 NW_019211458.1_4567805 0/0
## 1904 NW_019211458.1_4599456 <NA>
## 1905 NW_019211458.1_4601710 0/0
## 1906 NW_019211458.1_4718955 <NA>
## 1907 NW_019211458.1_4731554 <NA>
## 1908 NW_019211458.1_4736124 <NA>
## 1909 NW_019211458.1_4767724 <NA>
## 1910 NW_019211458.1_4791848 <NA>
## 1911 NW_019211458.1_4807639 <NA>
## 1912 NW_019211458.1_4853777 <NA>
## 1913 NW_019211458.1_4859954 0/0
## 1914 NW_019211458.1_4864440 0/0
## 1915 NW_019211458.1_4873032 0/0
## 1916 NW_019211458.1_4925166 0/0
## 1917 NW_019211458.1_5039431 <NA>
## 1918 NW_019211458.1_5044338 0/0
## 1919 NW_019211458.1_5063604 1/1
## 1920 NW_019211458.1_5086921 <NA>
## 1921 NW_019211458.1_5110485 0/0
## 1922 NW_019211458.1_5118361 <NA>
## 1923 NW_019211458.1_5128990 <NA>
## 1924 NW_019211458.1_5173802 0/0
## 1925 NW_019211458.1_5192566 0/0
## 1926 NW_019211458.1_5192913 <NA>
## 1927 NW_019211458.1_5199497 0/0
## 1928 NW_019211458.1_5269581 0/0
## 1929 NW_019211458.1_5299536 0/0
## 1930 NW_019211458.1_5299674 <NA>
## 1931 NW_019211458.1_5307804 <NA>
## 1932 NW_019211458.1_5374857 0/0
## 1933 NW_019211458.1_5479501 0/0
## 1934 NW_019211458.1_5494423 <NA>
## 1935 NW_019211458.1_5497967 <NA>
## 1936 NW_019211458.1_5503263 0/0
## 1937 NW_019211458.1_5504618 0/0
## 1938 NW_019211458.1_5511453 0/0
## 1939 NW_019211458.1_5523299 <NA>
## 1940 NW_019211458.1_5531127 <NA>
## 1941 NW_019211458.1_5535808 0/0
## 1942 NW_019211458.1_5546803 0/0
## 1943 NW_019211458.1_5547623 0/0
## 1944 NW_019211458.1_5557778 0/0
## 1945 NW_019211458.1_5574279 0/0
## 1946 NW_019211458.1_5577966 <NA>
## 1947 NW_019211458.1_5595097 <NA>
## 1948 NW_019211458.1_5600497 0/0
## 1949 NW_019211458.1_5600613 <NA>
## 1950 NW_019211458.1_5613882 0/0
## 1951 NW_019211458.1_5631738 <NA>
## 1952 NW_019211458.1_5680791 0/0
## 1953 NW_019211458.1_5694042 <NA>
## 1954 NW_019211458.1_5722833 0/0
## 1955 NW_019211458.1_5741658 0/0
## 1956 NW_019211458.1_5741703 <NA>
## 1957 NW_019211458.1_5804801 <NA>
## 1958 NW_019211458.1_5804897 <NA>
## 1959 NW_019211458.1_5823319 <NA>
## 1960 NW_019211458.1_5847489 <NA>
## 1961 NW_019211458.1_5903748 0/0
## 1962 NW_019211458.1_5912650 <NA>
## 1963 NW_019211458.1_5914918 0/0
## 1964 NW_019211458.1_5917146 0/0
## 1965 NW_019211458.1_5924336 0/0
## 1966 NW_019211458.1_5929745 0/0
## 1967 NW_019211458.1_5934807 <NA>
## 1968 NW_019211458.1_5935652 0/0
## 1969 NW_019211458.1_5937197 0/0
## 1970 NW_019211458.1_5938445 <NA>
## 1971 NW_019211458.1_5938601 0/0
## 1972 NW_019211458.1_5965772 0/0
## 1973 NW_019211458.1_5979869 <NA>
## 1974 NW_019211458.1_5992447 0/0
## 1975 NW_019211458.1_6009470 <NA>
## 1976 NW_019211458.1_6011168 0/0
## 1977 NW_019211458.1_6030794 0/0
## 1978 NW_019211458.1_6042570 <NA>
## 1979 NW_019211458.1_6057177 0/0
## 1980 NW_019211458.1_6058524 0/0
## 1981 NW_019211458.1_6058691 <NA>
## 1982 NW_019211458.1_6061677 0/0
## 1983 NW_019211458.1_6061843 0/0
## 1984 NW_019211458.1_6069829 0/0
## 1985 NW_019211458.1_6069908 0/0
## 1986 NW_019211458.1_6069982 0/0
## 1987 NW_019211458.1_6071819 <NA>
## 1988 NW_019211458.1_6144848 <NA>
## 1989 NW_019211458.1_6151080 <NA>
## 1990 NW_019211458.1_6191019 0/0
## 1991 NW_019211458.1_6191325 0/0
## 1992 NW_019211458.1_6202189 <NA>
## 1993 NW_019211458.1_6223063 <NA>
## 1994 NW_019211458.1_6448227 0/0
## 1995 NW_019211458.1_6565784 <NA>
## 1996 NW_019211458.1_6568869 <NA>
## 1997 NW_019211458.1_6572907 0/0
## 1998 NW_019211458.1_6576809 0/0
## 1999 NW_019211458.1_6598025 0/0
## 2000 NW_019211458.1_6600418 <NA>
## 2001 NW_019211458.1_6619620 <NA>
## 2002 NW_019211458.1_6625453 0/0
## 2003 NW_019211458.1_6640507 <NA>
## 2004 NW_019211458.1_6683613 0/0
## 2005 NW_019211458.1_6751037 0/0
## 2006 NW_019211458.1_6754083 0/0
## 2007 NW_019211458.1_6764061 0/0
## 2008 NW_019211458.1_6764939 <NA>
## 2009 NW_019211458.1_6778749 <NA>
## 2010 NW_019211458.1_6782718 0/0
## 2011 NW_019211458.1_6785629 <NA>
## 2012 NW_019211458.1_6795715 <NA>
## 2013 NW_019211458.1_6803452 <NA>
## 2014 NW_019211458.1_6811442 <NA>
## 2015 NW_019211458.1_6827923 <NA>
## 2016 NW_019211458.1_6842486 0/0
## 2017 NW_019211458.1_6885232 0/0
## 2018 NW_019211458.1_6933656 <NA>
## 2019 NW_019211458.1_6934656 0/0
## 2020 NW_019211458.1_6936668 0/0
## 2021 NW_019211458.1_6961890 <NA>
## 2022 NW_019211458.1_6961981 <NA>
## 2023 NW_019211458.1_6967444 <NA>
## 2024 NW_019211458.1_6980583 0/0
## 2025 NW_019211458.1_6981187 0/0
## 2026 NW_019211458.1_7208511 0/0
## 2027 NW_019211458.1_7235294 <NA>
## 2028 NW_019211458.1_7264863 <NA>
## 2029 NW_019211458.1_7265347 0/0
## 2030 NW_019211458.1_7631653 <NA>
## 2031 NW_019211458.1_7640452 <NA>
## 2032 NW_019211458.1_7774133 0/0
## 2033 NW_019211458.1_7786800 <NA>
## 2034 NW_019211458.1_7787879 0/0
## 2035 NW_019211458.1_7791296 <NA>
## 2036 NW_019211458.1_7799492 0/0
## 2037 NW_019211458.1_7866350 0/0
## 2038 NW_019211458.1_7888320 0/0
## 2039 NW_019211458.1_7889638 0/0
## 2040 NW_019211458.1_7954921 <NA>
## 2041 NW_019211458.1_7974357 0/0
## 2042 NW_019211458.1_7984851 0/0
## 2043 NW_019211458.1_8117696 0/0
## 2044 NW_019211458.1_8237922 <NA>
## 2045 NW_019211458.1_8239891 0/0
## 2046 NW_019211458.1_8251037 <NA>
## 2047 NW_019211458.1_8336887 <NA>
## 2048 NW_019211458.1_8369585 <NA>
## 2049 NW_019211458.1_8409414 0/0
## 2050 NW_019211458.1_8436682 0/0
## 2051 NW_019211458.1_8451711 0/0
## 2052 NW_019211458.1_8466673 0/0
## 2053 NW_019211458.1_8496190 <NA>
## 2054 NW_019211458.1_8504608 <NA>
## 2055 NW_019211458.1_8523779 0/0
## 2056 NW_019211458.1_8541409 <NA>
## 2057 NW_019211458.1_8568317 <NA>
## 2058 NW_019211458.1_8753223 0/0
## 2059 NW_019211458.1_8772588 0/0
## 2060 NW_019211458.1_8801693 0/0
## 2061 NW_019211458.1_8824576 0/0
## 2062 NW_019211458.1_8842277 0/0
## 2063 NW_019211458.1_8845053 <NA>
## 2064 NW_019211458.1_8877896 0/0
## 2065 NW_019211458.1_8887963 0/0
## 2066 NW_019211458.1_8896964 0/0
## 2067 NW_019211458.1_8898831 <NA>
## 2068 NW_019211458.1_8904369 0/0
## 2069 NW_019211458.1_8908098 <NA>
## 2070 NW_019211458.1_8908442 0/0
## 2071 NW_019211458.1_8915801 <NA>
## 2072 NW_019211458.1_8925533 <NA>
## 2073 NW_019211458.1_8929328 <NA>
## 2074 NW_019211458.1_8929373 <NA>
## 2075 NW_019211458.1_8941054 0/0
## 2076 NW_019211458.1_8954838 <NA>
## 2077 NW_019211458.1_11432208 0/0
## 2078 NW_019211458.1_11481565 <NA>
## 2079 NW_019211458.1_11481676 <NA>
## 2080 NW_019211458.1_11482583 0/0
## 2081 NW_019211458.1_11505070 0/0
## 2082 NW_019211458.1_11518442 0/0
## 2083 NW_019211458.1_11519731 0/0
## 2084 NW_019211458.1_11570514 0/0
## 2085 NW_019211458.1_11601340 0/0
## 2086 NW_019211458.1_11713555 <NA>
## 2087 NW_019211458.1_11713592 <NA>
## 2088 NW_019211458.1_11815375 0/1
## 2089 NW_019211458.1_12922078 <NA>
## 2090 NW_019211458.1_14116016 <NA>
## 2091 NW_019211458.1_15139477 0/0
## 2092 NW_019211458.1_15824535 0/0
## 2093 NW_019211458.1_15827758 0/0
## 2094 NW_019211458.1_15836665 <NA>
## 2095 NW_019211458.1_15859437 0/0
## 2096 NW_019211458.1_15878169 <NA>
## 2097 NW_019211458.1_15900717 <NA>
## 2098 NW_019211458.1_15905441 <NA>
## 2099 NW_019211458.1_15915235 0/0
## 2100 NW_019211458.1_15953469 <NA>
## 2101 NW_019211458.1_15991683 0/0
## 2102 NW_019211458.1_16001815 0/0
## 2103 NW_019211458.1_16007541 0/0
## 2104 NW_019211458.1_16023747 0/0
## 2105 NW_019211458.1_16023907 0/0
## 2106 NW_019211458.1_16024023 0/0
## 2107 NW_019211458.1_16066803 <NA>
## 2108 NW_019211458.1_16069923 0/0
## 2109 NW_019211458.1_16070177 0/0
## 2110 NW_019211458.1_16088760 0/0
## 2111 NW_019211458.1_16095165 0/0
## 2112 NW_019211458.1_16095833 <NA>
## 2113 NW_019211458.1_19013418 <NA>
## 2114 NW_019211458.1_19018068 0/0
## 2115 NW_019211458.1_19027189 <NA>
## 2116 NW_019211458.1_19028926 0/0
## 2117 NW_019211458.1_19031120 0/0
## 2118 NW_019211458.1_19050065 0/0
## 2119 NW_019211458.1_19050171 <NA>
## 2120 NW_019211458.1_19055947 <NA>
## 2121 NW_019211458.1_19066857 <NA>
## 2122 NW_019211458.1_19068543 0/0
## 2123 NW_019211458.1_19068845 0/0
## 2124 NW_019211458.1_19078895 0/0
## 2125 NW_019211458.1_19078983 0/0
## 2126 NW_019211458.1_19079095 0/0
## 2127 NW_019211458.1_19084806 0/0
## 2128 NW_019211458.1_19084873 0/0
## 2129 NW_019211458.1_19091764 <NA>
## 2130 NW_019211458.1_19092961 0/0
## 2131 NW_019211458.1_19095912 <NA>
## 2132 NW_019211458.1_19109324 <NA>
## 2133 NW_019211458.1_19113505 0/0
## 2134 NW_019211458.1_19119879 <NA>
## 2135 NW_019211458.1_19119961 <NA>
## 2136 NW_019211458.1_19124025 0/0
## 2137 NW_019211458.1_19181055 0/0
## 2138 NW_019211458.1_19181271 0/0
## 2139 NW_019211458.1_19196850 <NA>
## 2140 NW_019211458.1_19201768 <NA>
## 2141 NW_019211458.1_19215862 0/0
## 2142 NW_019211458.1_19232339 0/0
## 2143 NW_019211458.1_19253053 0/0
## 2144 NW_019211458.1_19257350 <NA>
## 2145 NW_019211458.1_19267663 0/0
## 2146 NW_019211458.1_19286737 <NA>
## 2147 NW_019211458.1_19294933 0/0
## 2148 NW_019211458.1_19338531 0/0
## 2149 NW_019211458.1_19354473 0/0
## 2150 NW_019211458.1_19354796 <NA>
## 2151 NW_019211458.1_19365119 0/0
## 2152 NW_019211458.1_19371573 0/0
## 2153 NW_019211458.1_19396211 <NA>
## 2154 NW_019211458.1_19397292 0/0
## 2155 NW_019211458.1_19400732 0/0
## 2156 NW_019211458.1_19405277 0/0
## 2157 NW_019211458.1_19405888 0/0
## 2158 NW_019211458.1_19596512 <NA>
## 2159 NW_019211458.1_19601614 0/0
## 2160 NW_019211458.1_19602033 0/0
## 2161 NW_019211458.1_19602247 0/0
## 2162 NW_019211458.1_19604009 0/0
## 2163 NW_019211458.1_19605186 <NA>
## 2164 NW_019211458.1_19616226 0/0
## 2165 NW_019211458.1_19624725 <NA>
## 2166 NW_019211458.1_19624973 <NA>
## 2167 NW_019211458.1_19626291 <NA>
## 2168 NW_019211458.1_19638190 0/0
## 2169 NW_019211458.1_19650797 0/0
## 2170 NW_019211458.1_19659026 <NA>
## 2171 NW_019211458.1_19660483 0/0
## 2172 NW_019211458.1_19660856 <NA>
## 2173 NW_019211458.1_19678845 <NA>
## 2174 NW_019211458.1_19725910 0/0
## 2175 NW_019211458.1_19730714 0/0
## 2176 NW_019211458.1_19735907 0/0
## 2177 NW_019211458.1_19737396 <NA>
## 2178 NW_019211458.1_19740580 <NA>
## 2179 NW_019211458.1_19836095 0/0
## 2180 NW_019211458.1_19837623 <NA>
## 2181 NW_019211458.1_19841141 <NA>
## 2182 NW_019211458.1_19842654 <NA>
## 2183 NW_019211458.1_19873331 0/0
## 2184 NW_019211458.1_19880995 <NA>
## 2185 NW_019211458.1_19881862 <NA>
## 2186 NW_019211458.1_19887566 <NA>
## 2187 NW_019211458.1_19901298 <NA>
## 2188 NW_019211458.1_19908020 <NA>
## 2189 NW_019211458.1_19908744 <NA>
## 2190 NW_019211458.1_19908783 <NA>
## 2191 NW_019211458.1_19920008 0/0
## 2192 NW_019211458.1_19992993 0/1
## 2193 NW_019211458.1_20838753 <NA>
## 2194 NW_019211458.1_21027260 0/0
## 2195 NW_019211458.1_27062739 <NA>
## 2196 NW_019211458.1_27189781 0/0
## 2197 NW_019211458.1_27206565 0/0
## 2198 NW_019211458.1_27211355 <NA>
## 2199 NW_019211458.1_27216651 0/0
## 2200 NW_019211458.1_27219710 0/0
## 2201 NW_019211458.1_27223068 <NA>
## 2202 NW_019211458.1_27240750 0/0
## 2203 NW_019211458.1_27258002 0/0
## 2204 NW_019211458.1_27265899 0/0
## 2205 NW_019211458.1_27270567 <NA>
## 2206 NW_019211458.1_27273649 0/0
## 2207 NW_019211458.1_27319003 0/0
## 2208 NW_019211458.1_27339853 <NA>
## 2209 NW_019211458.1_27350350 <NA>
## 2210 NW_019211458.1_27355400 0/0
## 2211 NW_019211458.1_27802083 0/0
## 2212 NW_019211458.1_27803679 <NA>
## 2213 NW_019211458.1_27913829 0/0
## 2214 NW_019211458.1_28179038 0/0
## 2215 NW_019211458.1_28183832 0/0
## 2216 NW_019211458.1_28187161 <NA>
## 2217 NW_019211458.1_28190473 0/0
## 2218 NW_019211458.1_28219707 0/0
## 2219 NW_019211458.1_28232864 0/0
## 2220 NW_019211458.1_29252933 <NA>
## 2221 NW_019211458.1_32332413 <NA>
## 2222 NW_019211458.1_32361992 <NA>
## 2223 NW_019211458.1_32398801 0/0
## 2224 NW_019211458.1_32400417 <NA>
## 2225 NW_019211458.1_32434604 0/0
## 2226 NW_019211458.1_32723983 0/0
## 2227 NW_019211458.1_32761400 0/0
## 2228 NW_019211458.1_32801232 <NA>
## 2229 NW_019211458.1_32819244 0/0
## 2230 NW_019211458.1_32824333 <NA>
## 2231 NW_019211458.1_32841501 0/0
## 2232 NW_019211458.1_32858811 0/0
## 2233 NW_019211458.1_32860345 0/0
## 2234 NW_019211458.1_32875388 <NA>
## 2235 NW_019211458.1_32879957 <NA>
## 2236 NW_019211458.1_32882940 0/0
## 2237 NW_019211458.1_33022961 0/0
## 2238 NW_019211458.1_33031453 <NA>
## 2239 NW_019211458.1_33063753 0/0
## 2240 NW_019211458.1_33064558 <NA>
## 2241 NW_019211458.1_33107893 <NA>
## 2242 NW_019211458.1_33109485 <NA>
## 2243 NW_019211458.1_33233555 0/1
## 2244 NW_019211458.1_33522104 <NA>
## 2245 NW_019211458.1_33804852 0/0
## 2246 NW_019211458.1_34302052 <NA>
## 2247 NW_019211458.1_34308469 <NA>
## 2248 NW_019211458.1_34352594 <NA>
## 2249 NW_019211458.1_34359656 <NA>
## 2250 NW_019211458.1_34365462 <NA>
## 2251 NW_019211458.1_34463012 0/0
## 2252 NW_019211458.1_35512875 0/0
## 2253 NW_019211458.1_35528359 0/0
## 2254 NW_019211458.1_35533449 <NA>
## 2255 NW_019211458.1_35540304 <NA>
## 2256 NW_019211458.1_35582307 0/0
## 2257 NW_019211458.1_35596997 <NA>
## 2258 NW_019211458.1_35598175 0/0
## 2259 NW_019211458.1_35610287 0/0
## 2260 NW_019211458.1_35761965 0/0
## 2261 NW_019211458.1_35828090 <NA>
## 2262 NW_019211458.1_36241182 <NA>
## 2263 NW_019211458.1_36242813 <NA>
## 2264 NW_019211458.1_36257868 0/0
## 2265 NW_019211458.1_36273185 0/0
## 2266 NW_019211458.1_36274829 <NA>
## 2267 NW_019211458.1_36299576 0/0
## 2268 NW_019211458.1_36323521 <NA>
## 2269 NW_019211458.1_36339994 <NA>
## 2270 NW_019211458.1_36345371 0/0
## 2271 NW_019211458.1_36346267 0/0
## 2272 NW_019211458.1_36353126 0/0
## 2273 NW_019211458.1_36354283 <NA>
## 2274 NW_019211458.1_36865715 0/0
## 2275 NW_019211458.1_36876584 0/0
## 2276 NW_019211458.1_36876619 <NA>
## 2277 NW_019211458.1_36878418 0/0
## 2278 NW_019211458.1_36913417 <NA>
## 2279 NW_019211458.1_36914408 0/0
## 2280 NW_019211458.1_36914923 <NA>
## 2281 NW_019211458.1_36915002 0/0
## 2282 NW_019211458.1_36959919 <NA>
## 2283 NW_019211458.1_36999868 <NA>
## 2284 NW_019211458.1_37022712 0/0
## 2285 NW_019211458.1_37041696 <NA>
## 2286 NW_019211458.1_37052376 <NA>
## 2287 NW_019211458.1_37073289 0/0
## 2288 NW_019211458.1_37120474 0/0
## 2289 NW_019211458.1_37217687 <NA>
## 2290 NW_019211458.1_37227564 0/0
## 2291 NW_019211458.1_37229032 0/0
## 2292 NW_019211458.1_37230357 <NA>
## 2293 NW_019211458.1_37243315 0/0
## 2294 NW_019211458.1_37283307 0/0
## 2295 NW_019211458.1_37283314 <NA>
## 2296 NW_019211458.1_37320842 0/0
## 2297 NW_019211458.1_37351519 0/0
## 2298 NW_019211458.1_37399338 0/0
## 2299 NW_019211458.1_37408537 0/0
## 2300 NW_019211458.1_37416973 <NA>
## 2301 NW_019211458.1_37423183 0/0
## 2302 NW_019211458.1_37451577 <NA>
## 2303 NW_019211458.1_37492743 <NA>
## 2304 NW_019211458.1_37557391 0/0
## 2305 NW_019211458.1_37621585 <NA>
## 2306 NW_019211458.1_37621641 <NA>
## 2307 NW_019211458.1_37679835 <NA>
## 2308 NW_019211458.1_37680170 <NA>
## 2309 NW_019211458.1_37682511 0/0
## 2310 NW_019211458.1_37684867 0/0
## 2311 NW_019211458.1_37714114 <NA>
## 2312 NW_019211458.1_37782487 0/0
## 2313 NW_019211458.1_37855379 0/0
## 2314 NW_019211458.1_37869307 <NA>
## 2315 NW_019211458.1_37903630 <NA>
## 2316 NW_019211458.1_37908386 <NA>
## 2317 NW_019211458.1_37950222 0/0
## 2318 NW_019211458.1_37953562 0/0
## 2319 NW_019211458.1_37960489 0/0
## 2320 NW_019211458.1_37978174 <NA>
## 2321 NW_019211458.1_38021076 0/0
## 2322 NW_019211458.1_38023027 0/0
## 2323 NW_019211458.1_38023475 0/0
## 2324 NW_019211458.1_38024944 <NA>
## 2325 NW_019211458.1_38025836 0/0
## 2326 NW_019211458.1_38027113 <NA>
## 2327 NW_019211458.1_38520712 <NA>
## 2328 NW_019211458.1_38538361 0/0
## 2329 NW_019211458.1_38548438 <NA>
## 2330 NW_019211458.1_38559480 0/0
## 2331 NW_019211458.1_38563467 0/0
## 2332 NW_019211458.1_38570601 0/0
## 2333 NW_019211458.1_38572714 0/0
## 2334 NW_019211458.1_38582213 <NA>
## 2335 NW_019211458.1_38590371 0/0
## 2336 NW_019211458.1_38636482 <NA>
## 2337 NW_019211458.1_38661992 <NA>
## 2338 NW_019211458.1_38699410 0/0
## 2339 NW_019211458.1_38719109 0/0
## 2340 NW_019211458.1_38932488 0/0
## 2341 NW_019211458.1_38932501 0/0
## 2342 NW_019211458.1_38932509 0/0
## 2343 NW_019211458.1_40638176 0/0
## 2344 NW_019211458.1_41788723 <NA>
## 2345 NW_019211458.1_41788885 0/0
## 2346 NW_019211458.1_41791574 0/0
## 2347 NW_019211458.1_41796965 0/0
## 2348 NW_019211458.1_41804698 0/0
## 2349 NW_019211458.1_41807060 0/0
## 2350 NW_019211458.1_41811488 0/0
## 2351 NW_019211458.1_41821532 0/0
## 2352 NW_019211458.1_41826848 0/0
## 2353 NW_019211458.1_41830093 0/0
## 2354 NW_019211458.1_41839896 <NA>
## 2355 NW_019211458.1_41887673 <NA>
## 2356 NW_019211459.1_510323 0/0
## 2357 NW_019211459.1_510350 0/0
## 2358 NW_019211459.1_532411 <NA>
## 2359 NW_019211459.1_624725 <NA>
## 2360 NW_019211459.1_632570 0/0
## 2361 NW_019211459.1_1347466 <NA>
## 2362 NW_019211459.1_1367453 <NA>
## 2363 NW_019211459.1_1700114 <NA>
## 2364 NW_019211459.1_1939843 <NA>
## 2365 NW_019211459.1_1941683 <NA>
## 2366 NW_019211459.1_2519543 0/0
## 2367 NW_019211459.1_2524265 <NA>
## 2368 NW_019211459.1_2542932 0/0
## 2369 NW_019211459.1_2574745 0/0
## 2370 NW_019211459.1_2598829 0/0
## 2371 NW_019211459.1_2610146 <NA>
## 2372 NW_019211459.1_2653117 0/0
## 2373 NW_019211459.1_2686546 <NA>
## 2374 NW_019211459.1_2693033 0/0
## 2375 NW_019211459.1_2696228 0/0
## 2376 NW_019211459.1_2748114 0/0
## 2377 NW_019211459.1_2795335 <NA>
## 2378 NW_019211459.1_2795797 0/0
## 2379 NW_019211459.1_2797543 0/0
## 2380 NW_019211459.1_2801120 <NA>
## 2381 NW_019211459.1_2804298 <NA>
## 2382 NW_019211459.1_2832639 0/0
## 2383 NW_019211459.1_2837319 0/0
## 2384 NW_019211459.1_2851154 0/0
## 2385 NW_019211459.1_2852395 <NA>
## 2386 NW_019211459.1_2859310 <NA>
## 2387 NW_019211459.1_2871438 0/0
## 2388 NW_019211459.1_2878083 0/0
## 2389 NW_019211459.1_2878105 0/0
## 2390 NW_019211459.1_2884416 0/0
## 2391 NW_019211459.1_2888854 0/0
## 2392 NW_019211459.1_2924150 0/0
## 2393 NW_019211459.1_2935854 <NA>
## 2394 NW_019211459.1_2947695 0/0
## 2395 NW_019211459.1_2978438 0/0
## 2396 NW_019211459.1_2979143 <NA>
## 2397 NW_019211459.1_2979163 <NA>
## 2398 NW_019211459.1_2981924 0/0
## 2399 NW_019211459.1_2998454 0/0
## 2400 NW_019211459.1_3006130 <NA>
## 2401 NW_019211459.1_3019585 <NA>
## 2402 NW_019211459.1_3032266 0/0
## 2403 NW_019211459.1_3048648 <NA>
## 2404 NW_019211459.1_3060683 0/0
## 2405 NW_019211459.1_3078000 0/0
## 2406 NW_019211459.1_3091163 <NA>
## 2407 NW_019211459.1_3105630 0/0
## 2408 NW_019211459.1_3140777 0/0
## 2409 NW_019211459.1_3162304 0/0
## 2410 NW_019211459.1_3175769 <NA>
## 2411 NW_019211459.1_3197656 0/0
## 2412 NW_019211459.1_3207550 <NA>
## 2413 NW_019211459.1_3218094 <NA>
## 2414 NW_019211459.1_3231117 0/0
## 2415 NW_019211459.1_3232596 <NA>
## 2416 NW_019211459.1_3247637 0/0
## 2417 NW_019211459.1_3251357 0/0
## 2418 NW_019211459.1_3251413 0/0
## 2419 NW_019211459.1_3263437 <NA>
## 2420 NW_019211459.1_3270237 0/0
## 2421 NW_019211459.1_3274126 <NA>
## 2422 NW_019211459.1_3277218 0/0
## 2423 NW_019211459.1_3283014 0/0
## 2424 NW_019211459.1_3285341 0/0
## 2425 NW_019211459.1_3289139 <NA>
## 2426 NW_019211459.1_3289174 0/0
## 2427 NW_019211459.1_3291738 <NA>
## 2428 NW_019211459.1_4021907 0/0
## 2429 NW_019211459.1_4228622 0/0
## 2430 NW_019211459.1_4259114 0/0
## 2431 NW_019211459.1_4262291 <NA>
## 2432 NW_019211459.1_4297197 0/1
## 2433 NW_019211459.1_4310364 <NA>
## 2434 NW_019211459.1_4320315 0/0
## 2435 NW_019211459.1_4397084 <NA>
## 2436 NW_019211459.1_4402585 <NA>
## 2437 NW_019211459.1_4422127 <NA>
## 2438 NW_019211459.1_4434550 0/0
## 2439 NW_019211459.1_4452017 <NA>
## 2440 NW_019211459.1_4457498 <NA>
## 2441 NW_019211459.1_4460795 0/0
## 2442 NW_019211459.1_4464236 <NA>
## 2443 NW_019211459.1_4473615 <NA>
## 2444 NW_019211459.1_4511023 0/0
## 2445 NW_019211459.1_4533525 0/0
## 2446 NW_019211459.1_4543765 0/0
## 2447 NW_019211459.1_4544044 0/0
## 2448 NW_019211459.1_4552610 <NA>
## 2449 NW_019211459.1_4563263 0/0
## 2450 NW_019211459.1_4568845 0/0
## 2451 NW_019211459.1_4599001 <NA>
## 2452 NW_019211459.1_4604999 0/0
## 2453 NW_019211459.1_4623715 0/0
## 2454 NW_019211459.1_4628790 <NA>
## 2455 NW_019211459.1_4629483 0/0
## 2456 NW_019211459.1_4629910 0/0
## 2457 NW_019211459.1_4640567 0/0
## 2458 NW_019211459.1_4648926 <NA>
## 2459 NW_019211459.1_4666027 <NA>
## 2460 NW_019211459.1_4673162 0/0
## 2461 NW_019211459.1_4679058 0/1
## 2462 NW_019211459.1_4730980 <NA>
## 2463 NW_019211459.1_4739012 <NA>
## 2464 NW_019211459.1_4758633 0/0
## 2465 NW_019211459.1_4875501 0/0
## 2466 NW_019211459.1_4955607 0/0
## 2467 NW_019211459.1_4962353 <NA>
## 2468 NW_019211459.1_5000908 0/0
## 2469 NW_019211459.1_5022473 <NA>
## 2470 NW_019211459.1_5033023 <NA>
## 2471 NW_019211459.1_5033099 <NA>
## 2472 NW_019211459.1_5036602 0/0
## 2473 NW_019211459.1_5083040 0/0
## 2474 NW_019211459.1_5128239 0/1
## 2475 NW_019211459.1_5128514 0/0
## 2476 NW_019211459.1_5167203 <NA>
## 2477 NW_019211459.1_5214712 <NA>
## 2478 NW_019211459.1_5275437 <NA>
## 2479 NW_019211459.1_6992579 0/0
## 2480 NW_019211459.1_7023185 0/0
## 2481 NW_019211459.1_8300019 <NA>
## 2482 NW_019211459.1_8899022 <NA>
## 2483 NW_019211459.1_8933702 <NA>
## 2484 NW_019211459.1_8935178 0/0
## 2485 NW_019211459.1_8935884 0/0
## 2486 NW_019211459.1_8937580 0/0
## 2487 NW_019211459.1_8946785 <NA>
## 2488 NW_019211459.1_8958209 0/0
## 2489 NW_019211459.1_8973609 0/0
## 2490 NW_019211459.1_8982284 <NA>
## 2491 NW_019211459.1_9008014 0/0
## 2492 NW_019211459.1_9059347 0/0
## 2493 NW_019211459.1_9241092 0/0
## 2494 NW_019211459.1_10087026 <NA>
## 2495 NW_019211459.1_10107747 0/1
## 2496 NW_019211459.1_10107794 0/1
## 2497 NW_019211459.1_10108541 0/0
## 2498 NW_019211459.1_10112298 <NA>
## 2499 NW_019211459.1_10121231 <NA>
## 2500 NW_019211459.1_10135711 <NA>
## 2501 NW_019211459.1_10146712 <NA>
## 2502 NW_019211459.1_10153761 <NA>
## 2503 NW_019211459.1_10154127 <NA>
## 2504 NW_019211459.1_10160385 0/0
## 2505 NW_019211459.1_10160600 0/0
## 2506 NW_019211459.1_10160707 0/0
## 2507 NW_019211459.1_10169116 <NA>
## 2508 NW_019211459.1_10182778 0/0
## 2509 NW_019211459.1_10198062 <NA>
## 2510 NW_019211459.1_10208896 0/0
## 2511 NW_019211459.1_10241638 0/0
## 2512 NW_019211459.1_10261575 0/0
## 2513 NW_019211459.1_10262865 <NA>
## 2514 NW_019211459.1_10265595 <NA>
## 2515 NW_019211459.1_10267105 <NA>
## 2516 NW_019211459.1_10281952 <NA>
## 2517 NW_019211459.1_10330557 0/0
## 2518 NW_019211459.1_10344060 <NA>
## 2519 NW_019211459.1_10350581 0/0
## 2520 NW_019211459.1_10364949 0/0
## 2521 NW_019211459.1_10366758 <NA>
## 2522 NW_019211459.1_10379197 <NA>
## 2523 NW_019211459.1_10380616 0/0
## 2524 NW_019211459.1_10380950 <NA>
## 2525 NW_019211459.1_10382268 <NA>
## 2526 NW_019211459.1_10382577 0/0
## 2527 NW_019211459.1_10391111 0/0
## 2528 NW_019211459.1_10399252 0/0
## 2529 NW_019211459.1_10433176 <NA>
## 2530 NW_019211459.1_10500630 <NA>
## 2531 NW_019211459.1_10550951 <NA>
## 2532 NW_019211459.1_10589130 <NA>
## 2533 NW_019211459.1_10592546 0/0
## 2534 NW_019211459.1_10606334 <NA>
## 2535 NW_019211459.1_10634012 0/0
## 2536 NW_019211459.1_10640340 0/0
## 2537 NW_019211459.1_10642636 <NA>
## 2538 NW_019211459.1_10646579 <NA>
## 2539 NW_019211459.1_10650724 0/0
## 2540 NW_019211459.1_10655757 <NA>
## 2541 NW_019211459.1_10672227 0/0
## 2542 NW_019211459.1_10674343 <NA>
## 2543 NW_019211459.1_10694592 0/0
## 2544 NW_019211459.1_10723436 0/0
## 2545 NW_019211459.1_10726444 <NA>
## 2546 NW_019211459.1_10727591 <NA>
## 2547 NW_019211459.1_10728287 <NA>
## 2548 NW_019211459.1_10736383 <NA>
## 2549 NW_019211459.1_10749078 0/0
## 2550 NW_019211459.1_10753960 0/0
## 2551 NW_019211459.1_10757714 0/0
## 2552 NW_019211459.1_10797461 0/0
## 2553 NW_019211459.1_10809098 0/0
## 2554 NW_019211459.1_10816645 0/0
## 2555 NW_019211459.1_10828995 <NA>
## 2556 NW_019211459.1_10836412 0/0
## 2557 NW_019211459.1_10839862 0/0
## 2558 NW_019211459.1_10840342 <NA>
## 2559 NW_019211459.1_10844485 <NA>
## 2560 NW_019211459.1_10852954 <NA>
## 2561 NW_019211459.1_10870232 <NA>
## 2562 NW_019211459.1_10873062 <NA>
## 2563 NW_019211459.1_10893246 <NA>
## 2564 NW_019211459.1_10903958 <NA>
## 2565 NW_019211459.1_10913313 <NA>
## 2566 NW_019211459.1_10919810 0/0
## 2567 NW_019211459.1_10920443 0/0
## 2568 NW_019211459.1_10928329 0/0
## 2569 NW_019211459.1_10938695 <NA>
## 2570 NW_019211459.1_10981677 0/0
## 2571 NW_019211459.1_11003838 <NA>
## 2572 NW_019211459.1_11039657 <NA>
## 2573 NW_019211459.1_11048749 0/0
## 2574 NW_019211459.1_11052886 <NA>
## 2575 NW_019211459.1_11061956 0/0
## 2576 NW_019211459.1_11069914 <NA>
## 2577 NW_019211459.1_11071223 0/0
## 2578 NW_019211459.1_11073068 <NA>
## 2579 NW_019211459.1_11086706 0/0
## 2580 NW_019211459.1_11099083 <NA>
## 2581 NW_019211459.1_11126677 0/0
## 2582 NW_019211459.1_11154036 <NA>
## 2583 NW_019211459.1_11157686 0/0
## 2584 NW_019211459.1_11223902 0/0
## 2585 NW_019211459.1_11225067 0/0
## 2586 NW_019211459.1_11316868 <NA>
## 2587 NW_019211459.1_11353270 <NA>
## 2588 NW_019211459.1_11354712 <NA>
## 2589 NW_019211459.1_11357258 <NA>
## 2590 NW_019211459.1_11361710 0/0
## 2591 NW_019211459.1_11374565 <NA>
## 2592 NW_019211459.1_11412125 0/0
## 2593 NW_019211459.1_11414731 0/0
## 2594 NW_019211459.1_11415983 0/0
## 2595 NW_019211459.1_11416038 0/0
## 2596 NW_019211459.1_11420923 0/0
## 2597 NW_019211459.1_11432749 0/0
## 2598 NW_019211459.1_11432839 0/0
## 2599 NW_019211459.1_11452783 <NA>
## 2600 NW_019211459.1_11458297 <NA>
## 2601 NW_019211459.1_11458978 <NA>
## 2602 NW_019211459.1_11472264 <NA>
## 2603 NW_019211459.1_11475340 <NA>
## 2604 NW_019211459.1_11475505 0/0
## 2605 NW_019211459.1_11504195 <NA>
## 2606 NW_019211459.1_12521869 <NA>
## 2607 NW_019211459.1_12563651 0/0
## 2608 NW_019211459.1_12567074 0/0
## 2609 NW_019211459.1_12578233 <NA>
## 2610 NW_019211459.1_14820804 <NA>
## 2611 NW_019211459.1_14831426 0/0
## 2612 NW_019211459.1_14872279 0/0
## 2613 NW_019211459.1_14947379 <NA>
## 2614 NW_019211459.1_14950937 <NA>
## 2615 NW_019211459.1_14951266 <NA>
## 2616 NW_019211459.1_14952178 <NA>
## 2617 NW_019211459.1_14956802 <NA>
## 2618 NW_019211459.1_14962368 <NA>
## 2619 NW_019211459.1_14980326 <NA>
## 2620 NW_019211459.1_15000266 0/0
## 2621 NW_019211459.1_15000613 <NA>
## 2622 NW_019211459.1_15003096 0/0
## 2623 NW_019211459.1_15003526 <NA>
## 2624 NW_019211459.1_15003533 <NA>
## 2625 NW_019211459.1_15008193 <NA>
## 2626 NW_019211459.1_15011154 <NA>
## 2627 NW_019211459.1_15024701 0/0
## 2628 NW_019211459.1_15036510 0/0
## 2629 NW_019211459.1_15077907 0/0
## 2630 NW_019211459.1_15124384 <NA>
## 2631 NW_019211459.1_15155897 0/0
## 2632 NW_019211459.1_15158179 <NA>
## 2633 NW_019211459.1_15159187 <NA>
## 2634 NW_019211459.1_15177454 <NA>
## 2635 NW_019211459.1_15193479 <NA>
## 2636 NW_019211459.1_15198257 <NA>
## 2637 NW_019211459.1_15198999 0/0
## 2638 NW_019211459.1_15200811 <NA>
## 2639 NW_019211459.1_15216901 0/0
## 2640 NW_019211459.1_15216918 0/0
## 2641 NW_019211459.1_15217012 0/0
## 2642 NW_019211459.1_15220950 0/0
## 2643 NW_019211459.1_15224199 0/0
## 2644 NW_019211459.1_15240515 0/0
## 2645 NW_019211459.1_15240547 0/0
## 2646 NW_019211459.1_15249165 0/0
## 2647 NW_019211459.1_15258611 <NA>
## 2648 NW_019211459.1_15259791 0/0
## 2649 NW_019211459.1_15264885 0/0
## 2650 NW_019211459.1_15266513 0/0
## 2651 NW_019211459.1_15298609 <NA>
## 2652 NW_019211459.1_15308432 0/0
## 2653 NW_019211459.1_15322608 0/0
## 2654 NW_019211459.1_15326372 0/0
## 2655 NW_019211459.1_15349553 0/0
## 2656 NW_019211459.1_15377453 0/0
## 2657 NW_019211459.1_15377886 <NA>
## 2658 NW_019211459.1_15378537 <NA>
## 2659 NW_019211459.1_15395119 0/0
## 2660 NW_019211459.1_15417953 0/0
## 2661 NW_019211459.1_15419511 <NA>
## 2662 NW_019211459.1_15449351 <NA>
## 2663 NW_019211459.1_15451862 0/0
## 2664 NW_019211459.1_15478785 0/0
## 2665 NW_019211459.1_16507456 0/0
## 2666 NW_019211459.1_16510496 0/0
## 2667 NW_019211459.1_16510686 0/0
## 2668 NW_019211459.1_16546534 0/0
## 2669 NW_019211459.1_16574418 <NA>
## 2670 NW_019211459.1_16579165 <NA>
## 2671 NW_019211459.1_16579796 <NA>
## 2672 NW_019211459.1_16597717 <NA>
## 2673 NW_019211459.1_16649658 0/0
## 2674 NW_019211459.1_16666296 <NA>
## 2675 NW_019211459.1_16668755 <NA>
## 2676 NW_019211459.1_16705278 0/0
## 2677 NW_019211459.1_16757856 0/0
## 2678 NW_019211459.1_16821859 <NA>
## 2679 NW_019211459.1_16849069 0/0
## 2680 NW_019211459.1_16906948 <NA>
## 2681 NW_019211459.1_17455485 0/0
## 2682 NW_019211459.1_17473661 <NA>
## 2683 NW_019211459.1_17501696 0/0
## 2684 NW_019211459.1_17502670 0/0
## 2685 NW_019211459.1_17504758 0/0
## 2686 NW_019211459.1_17506325 <NA>
## 2687 NW_019211459.1_17550842 0/0
## 2688 NW_019211459.1_17550866 0/0
## 2689 NW_019211459.1_17551060 0/0
## 2690 NW_019211459.1_17704709 0/0
## 2691 NW_019211459.1_17747882 0/0
## 2692 NW_019211459.1_17754516 <NA>
## 2693 NW_019211459.1_17757452 <NA>
## 2694 NW_019211459.1_17757936 0/0
## 2695 NW_019211459.1_17762475 0/0
## 2696 NW_019211459.1_17783200 <NA>
## 2697 NW_019211459.1_17783312 <NA>
## 2698 NW_019211459.1_17853495 0/0
## 2699 NW_019211459.1_17853976 <NA>
## 2700 NW_019211459.1_17868550 0/0
## 2701 NW_019211459.1_17873236 <NA>
## 2702 NW_019211459.1_18278233 0/0
## 2703 NW_019211459.1_18322494 0/0
## 2704 NW_019211459.1_18396482 <NA>
## 2705 NW_019211459.1_18513643 <NA>
## 2706 NW_019211459.1_18518165 <NA>
## 2707 NW_019211459.1_19392838 0/0
## 2708 NW_019211459.1_19657853 <NA>
## 2709 NW_019211459.1_19685137 <NA>
## 2710 NW_019211459.1_19711917 <NA>
## 2711 NW_019211459.1_19736608 <NA>
## 2712 NW_019211459.1_19758847 <NA>
## 2713 NW_019211459.1_19778873 0/0
## 2714 NW_019211459.1_19838151 <NA>
## 2715 NW_019211459.1_19866491 <NA>
## 2716 NW_019211459.1_19866551 <NA>
## 2717 NW_019211459.1_19866562 <NA>
## 2718 NW_019211459.1_21101757 <NA>
## 2719 NW_019211459.1_21104520 0/0
## 2720 NW_019211459.1_21105358 <NA>
## 2721 NW_019211459.1_21108346 <NA>
## 2722 NW_019211459.1_21110045 <NA>
## 2723 NW_019211459.1_21127498 <NA>
## 2724 NW_019211459.1_21138929 <NA>
## 2725 NW_019211459.1_21145296 <NA>
## 2726 NW_019211459.1_21145698 <NA>
## 2727 NW_019211459.1_21147876 <NA>
## 2728 NW_019211459.1_21147886 <NA>
## 2729 NW_019211459.1_21147979 <NA>
## 2730 NW_019211459.1_21172405 <NA>
## 2731 NW_019211459.1_21177339 <NA>
## 2732 NW_019211459.1_21223719 <NA>
## 2733 NW_019211459.1_21225338 0/0
## 2734 NW_019211459.1_21229882 <NA>
## 2735 NW_019211459.1_21230817 0/0
## 2736 NW_019211459.1_21230879 0/0
## 2737 NW_019211459.1_21233465 0/0
## 2738 NW_019211459.1_21277299 <NA>
## 2739 NW_019211459.1_21280179 <NA>
## 2740 NW_019211459.1_21291378 <NA>
## 2741 NW_019211459.1_21291829 0/0
## 2742 NW_019211459.1_21363440 0/0
## 2743 NW_019211459.1_21409066 0/0
## 2744 NW_019211459.1_21560113 <NA>
## 2745 NW_019211459.1_21567705 <NA>
## 2746 NW_019211459.1_22189151 <NA>
## 2747 NW_019211459.1_23806425 <NA>
## 2748 NW_019211459.1_23806429 <NA>
## 2749 NW_019211459.1_24412142 0/0
## 2750 NW_019211459.1_24439494 <NA>
## 2751 NW_019211459.1_24450964 0/0
## 2752 NW_019211459.1_24454464 <NA>
## 2753 NW_019211459.1_24542874 <NA>
## 2754 NW_019211459.1_24590387 0/0
## 2755 NW_019211459.1_24604843 0/0
## 2756 NW_019211459.1_24817465 <NA>
## 2757 NW_019211459.1_24875829 <NA>
## 2758 NW_019211459.1_24877313 <NA>
## 2759 NW_019211459.1_24879194 <NA>
## 2760 NW_019211459.1_24887150 <NA>
## 2761 NW_019211459.1_24902074 0/0
## 2762 NW_019211459.1_24910069 <NA>
## 2763 NW_019211459.1_24919601 <NA>
## 2764 NW_019211459.1_24920560 <NA>
## 2765 NW_019211459.1_24935585 0/0
## 2766 NW_019211459.1_24937340 0/0
## 2767 NW_019211459.1_24943479 0/0
## 2768 NW_019211459.1_24947785 0/0
## 2769 NW_019211459.1_24980552 0/0
## 2770 NW_019211459.1_25026504 <NA>
## 2771 NW_019211459.1_25145986 0/0
## 2772 NW_019211459.1_25146371 0/0
## 2773 NW_019211459.1_25302867 <NA>
## 2774 NW_019211459.1_25307360 0/0
## 2775 NW_019211459.1_25308155 0/0
## 2776 NW_019211459.1_25314277 0/0
## 2777 NW_019211459.1_25327841 <NA>
## 2778 NW_019211459.1_25339976 <NA>
## 2779 NW_019211459.1_25349143 0/0
## 2780 NW_019211459.1_25368088 0/0
## 2781 NW_019211459.1_25402092 <NA>
## 2782 NW_019211459.1_25404651 0/0
## 2783 NW_019211459.1_25421345 <NA>
## 2784 NW_019211459.1_25423717 <NA>
## 2785 NW_019211459.1_25590948 0/0
## 2786 NW_019211459.1_25695764 0/0
## 2787 NW_019211459.1_25697756 <NA>
## 2788 NW_019211459.1_26994249 0/0
## 2789 NW_019211459.1_27024200 <NA>
## 2790 NW_019211459.1_27069604 <NA>
## 2791 NW_019211459.1_27078644 <NA>
## 2792 NW_019211459.1_27088479 <NA>
## 2793 NW_019211459.1_27095147 0/0
## 2794 NW_019211459.1_27101407 <NA>
## 2795 NW_019211459.1_27130898 0/0
## 2796 NW_019211459.1_27141459 0/0
## 2797 NW_019211459.1_27171958 <NA>
## 2798 NW_019211459.1_27176129 0/0
## 2799 NW_019211459.1_27191312 <NA>
## 2800 NW_019211459.1_27195799 <NA>
## 2801 NW_019211459.1_27209794 0/0
## 2802 NW_019211459.1_27313697 <NA>
## 2803 NW_019211459.1_27337349 <NA>
## 2804 NW_019211459.1_27366708 <NA>
## 2805 NW_019211459.1_27397853 0/0
## 2806 NW_019211459.1_28727429 0/0
## 2807 NW_019211459.1_28729763 0/0
## 2808 NW_019211459.1_28729772 0/0
## 2809 NW_019211459.1_28729914 <NA>
## 2810 NW_019211459.1_28731221 <NA>
## 2811 NW_019211459.1_28731819 <NA>
## 2812 NW_019211459.1_28764056 0/0
## 2813 NW_019211459.1_28787309 0/0
## 2814 NW_019211459.1_28823791 <NA>
## 2815 NW_019211459.1_28824457 0/0
## 2816 NW_019211459.1_28840633 <NA>
## 2817 NW_019211459.1_28848882 <NA>
## 2818 NW_019211459.1_28868041 0/0
## 2819 NW_019211459.1_28890187 0/0
## 2820 NW_019211459.1_28890903 0/0
## 2821 NW_019211459.1_28918063 0/0
## 2822 NW_019211459.1_28924314 <NA>
## 2823 NW_019211459.1_28924554 <NA>
## 2824 NW_019211459.1_28935370 0/0
## 2825 NW_019211459.1_28941202 0/0
## 2826 NW_019211459.1_28946902 0/0
## 2827 NW_019211459.1_28947280 <NA>
## 2828 NW_019211459.1_28954103 0/0
## 2829 NW_019211459.1_28960351 0/0
## 2830 NW_019211459.1_28985557 0/0
## 2831 NW_019211459.1_28986755 <NA>
## 2832 NW_019211459.1_28988095 <NA>
## 2833 NW_019211459.1_28989525 0/0
## 2834 NW_019211459.1_29018094 <NA>
## 2835 NW_019211459.1_29030427 0/0
## 2836 NW_019211459.1_29429319 <NA>
## 2837 NW_019211459.1_31582812 <NA>
## 2838 NW_019211460.1_12931634 0/0
## 2839 NW_019211460.1_15044225 <NA>
## 2840 NW_019211460.1_15062783 <NA>
## 2841 NW_019211460.1_15074460 <NA>
## 2842 NW_019211460.1_15080077 <NA>
## 2843 NW_019211460.1_15093868 0/0
## 2844 NW_019211460.1_15099849 0/0
## 2845 NW_019211460.1_15105019 0/0
## 2846 NW_019211460.1_15106952 0/0
## 2847 NW_019211460.1_15111841 0/0
## 2848 NW_019211460.1_15112998 <NA>
## 2849 NW_019211460.1_15137622 0/0
## 2850 NW_019211460.1_15137703 <NA>
## 2851 NW_019211460.1_15142577 <NA>
## 2852 NW_019211460.1_15143847 <NA>
## 2853 NW_019211460.1_15156152 0/0
## 2854 NW_019211460.1_15160358 <NA>
## 2855 NW_019211460.1_15166292 <NA>
## 2856 NW_019211460.1_15166629 0/0
## 2857 NW_019211460.1_15170564 0/0
## 2858 NW_019211460.1_15178081 <NA>
## 2859 NW_019211460.1_15664598 <NA>
## 2860 NW_019211460.1_16839024 0/0
## 2861 NW_019211460.1_16863056 <NA>
## 2862 NW_019211460.1_16866004 0/0
## 2863 NW_019211460.1_16866343 <NA>
## 2864 NW_019211460.1_16866368 <NA>
## 2865 NW_019211460.1_16877968 0/0
## 2866 NW_019211460.1_16885811 0/0
## 2867 NW_019211460.1_16899481 <NA>
## 2868 NW_019211460.1_16904981 <NA>
## 2869 NW_019211460.1_16908001 0/0
## 2870 NW_019211460.1_16913047 0/0
## 2871 NW_019211460.1_16913688 0/0
## 2872 NW_019211460.1_16920377 0/0
## 2873 NW_019211460.1_16922744 0/0
## 2874 NW_019211460.1_16938607 0/0
## 2875 NW_019211460.1_16941948 <NA>
## 2876 NW_019211460.1_16941957 <NA>
## 2877 NW_019211460.1_16941983 0/0
## 2878 NW_019211460.1_16951890 0/0
## 2879 NW_019211460.1_16960747 0/0
## 2880 NW_019211460.1_16962543 <NA>
## 2881 NW_019211460.1_16969978 0/0
## 2882 NW_019211460.1_16969985 0/0
## 2883 NW_019211460.1_16989117 0/0
## 2884 NW_019211460.1_16989571 <NA>
## 2885 NW_019211460.1_18071833 <NA>
## 2886 NW_019211460.1_18105555 <NA>
## 2887 NW_019211460.1_18177915 <NA>
## 2888 NW_019211460.1_18188131 <NA>
## 2889 NW_019211460.1_18206694 <NA>
## 2890 NW_019211460.1_18207348 <NA>
## 2891 NW_019211460.1_19105952 0/0
## 2892 NW_019211460.1_19106353 <NA>
## 2893 NW_019211460.1_19108135 0/0
## 2894 NW_019211460.1_19123354 <NA>
## 2895 NW_019211460.1_19124448 <NA>
## 2896 NW_019211460.1_19131285 0/0
## 2897 NW_019211460.1_19137831 <NA>
## 2898 NW_019211460.1_19152321 <NA>
## 2899 NW_019211460.1_19153419 0/0
## 2900 NW_019211460.1_19165257 0/0
## 2901 NW_019211460.1_22806463 0/0
## 2902 NW_019211460.1_23247293 0/0
## 2903 NW_019211460.1_23344582 0/0
## 2904 NW_019211460.1_23802699 <NA>
## 2905 NW_019211460.1_24002852 0/0
## 2906 NW_019211460.1_24018316 0/0
## 2907 NW_019211460.1_24027395 <NA>
## 2908 NW_019211460.1_24042122 <NA>
## 2909 NW_019211460.1_24053598 0/0
## 2910 NW_019211460.1_24724896 <NA>
## 2911 NW_019211460.1_24825526 <NA>
## 2912 NW_019211460.1_24841003 <NA>
## 2913 NW_019211460.1_24841986 0/0
## 2914 NW_019211460.1_24878441 0/0
## 2915 NW_019211460.1_24885118 <NA>
## 2916 NW_019211460.1_24887517 <NA>
## 2917 NW_019211460.1_24915912 0/0
## 2918 NW_019211460.1_25003419 0/0
## 2919 NW_019211460.1_25014341 <NA>
## 2920 NW_019211460.1_25018135 <NA>
## 2921 NW_019211460.1_25022889 0/0
## 2922 NW_019211460.1_25030682 0/0
## 2923 NW_019211460.1_25045908 0/0
## 2924 NW_019211460.1_25127481 <NA>
## 2925 NW_019211460.1_25128531 <NA>
## 2926 NW_019211460.1_25161208 <NA>
## 2927 NW_019211460.1_25178534 0/0
## 2928 NW_019211460.1_25182753 0/0
## 2929 NW_019211460.1_25195380 0/0
## 2930 NW_019211460.1_25195395 0/0
## 2931 NW_019211460.1_25198752 <NA>
## 2932 NW_019211460.1_25198828 <NA>
## 2933 NW_019211460.1_25209591 <NA>
## 2934 NW_019211460.1_25232531 0/0
## 2935 NW_019211460.1_25247339 0/0
## 2936 NW_019211460.1_25249614 <NA>
## 2937 NW_019211460.1_25256932 0/0
## 2938 NW_019211460.1_25262705 <NA>
## 2939 NW_019211460.1_25276446 <NA>
## 2940 NW_019211460.1_25276553 <NA>
## 2941 NW_019211460.1_25276669 0/0
## 2942 NW_019211460.1_25278561 0/0
## 2943 NW_019211460.1_25287719 0/0
## 2944 NW_019211460.1_25288753 <NA>
## 2945 NW_019211460.1_25295444 <NA>
## 2946 NW_019211460.1_25307870 0/0
## 2947 NW_019211460.1_25320490 0/0
## 2948 NW_019211460.1_25323653 0/0
## 2949 NW_019211460.1_25339648 0/0
## 2950 NW_019211460.1_25350152 <NA>
## 2951 NW_019211460.1_25361016 <NA>
## 2952 NW_019211460.1_25368991 0/0
## 2953 NW_019211460.1_25379747 <NA>
## 2954 NW_019211460.1_25384542 0/0
## 2955 NW_019211460.1_25396698 <NA>
## 2956 NW_019211460.1_25416104 0/0
## 2957 NW_019211460.1_25424159 0/0
## 2958 NW_019211460.1_25433021 0/0
## 2959 NW_019211460.1_25439608 <NA>
## 2960 NW_019211460.1_25439872 0/0
## 2961 NW_019211460.1_25573640 <NA>
## 2962 NW_019211460.1_25574326 <NA>
## 2963 NW_019211460.1_25576696 <NA>
## 2964 NW_019211460.1_25577312 0/0
## 2965 NW_019211460.1_25577476 <NA>
## 2966 NW_019211460.1_25579474 <NA>
## 2967 NW_019211460.1_26171391 0/0
## 2968 NW_019211460.1_26171400 0/0
## 2969 NW_019211460.1_26182516 0/0
## 2970 NW_019211460.1_26187482 <NA>
## 2971 NW_019211460.1_26196068 0/0
## 2972 NW_019211460.1_26197331 0/0
## 2973 NW_019211460.1_26198191 0/0
## 2974 NW_019211460.1_26232280 0/0
## 2975 NW_019211460.1_26232919 0/0
## 2976 NW_019211460.1_26235727 0/0
## 2977 NW_019211460.1_26243136 0/0
## 2978 NW_019211460.1_29665252 <NA>
## 2979 NW_019211460.1_29716761 <NA>
## 2980 NW_019211460.1_31235175 <NA>
## 2981 NW_019211460.1_31244709 <NA>
## 2982 NW_019211460.1_31269312 0/0
## 2983 NW_019211460.1_31272915 <NA>
## 2984 NW_019211460.1_31320858 0/0
## 2985 NW_019211460.1_31376306 0/0
## 2986 NW_019211460.1_31380075 <NA>
## 2987 NW_019211460.1_31387222 0/0
## 2988 NW_019211460.1_31387393 <NA>
## 2989 NW_019211460.1_31389335 0/0
## 2990 NW_019211460.1_31393856 0/0
## 2991 NW_019211460.1_31398725 0/0
## 2992 NW_019211460.1_31399415 0/0
## 2993 NW_019211460.1_31410941 0/0
## 2994 NW_019211460.1_31414455 <NA>
## 2995 NW_019211460.1_31419681 0/0
## 2996 NW_019211460.1_31419694 0/0
## 2997 NW_019211460.1_31425021 <NA>
## 2998 NW_019211460.1_31427558 0/0
## 2999 NW_019211460.1_31429025 0/0
## 3000 NW_019211460.1_31431478 0/0
## 3001 NW_019211460.1_31432913 0/0
## 3002 NW_019211460.1_31436979 <NA>
## 3003 NW_019211460.1_31438471 0/0
## 3004 NW_019211460.1_31442110 <NA>
## 3005 NW_019211460.1_31488643 <NA>
## 3006 NW_019211460.1_31489789 0/0
## 3007 NW_019211460.1_31491343 <NA>
## 3008 NW_019211460.1_31492247 0/1
## 3009 NW_019211460.1_31493843 <NA>
## 3010 NW_019211460.1_31495429 0/0
## 3011 NW_019211460.1_31503669 0/0
## 3012 NW_019211460.1_31503834 0/0
## 3013 NW_019211460.1_31513600 <NA>
## 3014 NW_019211460.1_31597890 <NA>
## 3015 NW_019211460.1_31599116 0/0
## 3016 NW_019211460.1_31601492 0/0
## 3017 NW_019211460.1_31610769 <NA>
## 3018 NW_019211460.1_31648212 <NA>
## 3019 NW_019211460.1_31648968 0/0
## 3020 NW_019211460.1_31652326 <NA>
## 3021 NW_019211460.1_31674644 <NA>
## 3022 NW_019211460.1_31675965 0/0
## 3023 NW_019211460.1_31676774 0/0
## 3024 NW_019211460.1_31696355 0/0
## 3025 NW_019211460.1_31703516 <NA>
## 3026 NW_019211460.1_31713732 0/0
## 3027 NW_019211460.1_31722091 <NA>
## 3028 NW_019211460.1_31722242 0/0
## 3029 NW_019211460.1_31728935 <NA>
## 3030 NW_019211460.1_31741624 0/0
## 3031 NW_019211460.1_32089125 <NA>
## 3032 NW_019211460.1_32760140 0/0
## 3033 NW_019211460.1_32783801 0/0
## 3034 NW_019211460.1_32785599 <NA>
## 3035 NW_019211460.1_32791825 <NA>
## 3036 NW_019211460.1_32791840 <NA>
## 3037 NW_019211460.1_32997494 0/0
## 3038 NW_019211460.1_33010381 0/0
## 3039 NW_019211460.1_33032400 0/0
## 3040 NW_019211460.1_33053490 <NA>
## 3041 NW_019211460.1_33057515 <NA>
## 3042 NW_019211460.1_33060141 0/0
## 3043 NW_019211460.1_33062371 0/0
## 3044 NW_019211460.1_33062504 0/0
## 3045 NW_019211460.1_33062519 0/0
## 3046 NW_019211460.1_33066045 <NA>
## 3047 NW_019211460.1_33076181 <NA>
## 3048 NW_019211460.1_33097115 0/0
## 3049 NW_019211460.1_33287487 0/0
## 3050 NW_019211460.1_33287839 <NA>
## 3051 NW_019211460.1_33304817 0/0
## 3052 NW_019211460.1_33309273 0/0
## 3053 NW_019211460.1_33310506 0/0
## 3054 NW_019211460.1_33312790 0/0
## 3055 NW_019211460.1_33313369 <NA>
## 3056 NW_019211460.1_33359694 <NA>
## 3057 NW_019211460.1_33366977 0/0
## 3058 NW_019211460.1_33379971 0/0
## 3059 NW_019211460.1_33404378 0/0
## 3060 NW_019211460.1_33415732 <NA>
## 3061 NW_019211460.1_33428353 <NA>
## 3062 NW_019211460.1_33429143 <NA>
## 3063 NW_019211460.1_33433042 <NA>
## 3064 NW_019211460.1_33495253 0/0
## 3065 NW_019211460.1_33522404 0/0
## 3066 NW_019211460.1_33522600 0/0
## 3067 NW_019211460.1_33527995 <NA>
## 3068 NW_019211460.1_33529980 <NA>
## 3069 NW_019211460.1_33530202 0/0
## 3070 NW_019211460.1_33534021 <NA>
## 3071 NW_019211460.1_33535617 <NA>
## 3072 NW_019211460.1_33917439 <NA>
## 3073 NW_019211460.1_34014687 0/0
## 3074 NW_019211460.1_34041959 <NA>
## 3075 NW_019211460.1_34043457 0/0
## 3076 NW_019211460.1_34058147 <NA>
## 3077 NW_019211460.1_34058870 0/0
## 3078 NW_019211460.1_34062816 <NA>
## 3079 NW_019211460.1_34070680 0/0
## 3080 NW_019211460.1_34126479 <NA>
## 3081 NW_019211460.1_34137993 0/0
## 3082 NW_019211460.1_34138303 0/0
## 3083 NW_019211460.1_34142318 0/0
## 3084 NW_019211460.1_34154512 <NA>
## 3085 NW_019211460.1_34165787 0/0
## 3086 NW_019211460.1_34205199 0/1
## 3087 NW_019211460.1_34217942 0/0
## 3088 NW_019211460.1_34218229 0/0
## 3089 NW_019211460.1_34232725 <NA>
## 3090 NW_019211460.1_34245101 <NA>
## 3091 NW_019211460.1_34261994 0/0
## 3092 NW_019211460.1_34262025 0/0
## 3093 NW_019211460.1_34268652 <NA>
## 3094 NW_019211460.1_34291766 <NA>
## 3095 NW_019211460.1_34292705 <NA>
## 3096 NW_019211460.1_34294018 0/0
## 3097 NW_019211460.1_34300158 0/0
## 3098 NW_019211460.1_34300503 0/0
## 3099 NW_019211460.1_34301950 <NA>
## 3100 NW_019211460.1_34308067 0/0
## 3101 NW_019211460.1_34309708 0/0
## 3102 NW_019211460.1_34318077 0/0
## 3103 NW_019211460.1_34324796 <NA>
## 3104 NW_019211460.1_34384750 0/0
## 3105 NW_019211460.1_34586356 <NA>
## 3106 NW_019211460.1_36084650 0/0
## 3107 NW_019211460.1_36089631 0/0
## 3108 NW_019211460.1_36106639 <NA>
## 3109 NW_019211460.1_36107351 <NA>
## 3110 NW_019211460.1_36116166 <NA>
## 3111 NW_019211460.1_36143041 0/0
## 3112 NW_019211460.1_36151630 0/0
## 3113 NW_019211460.1_36156137 <NA>
## 3114 NW_019211460.1_36170384 0/0
## 3115 NW_019211460.1_36183796 0/0
## 3116 NW_019211460.1_36185700 0/0
## 3117 NW_019211460.1_36186688 0/0
## 3118 NW_019211460.1_36189677 0/0
## 3119 NW_019211460.1_36415971 <NA>
## 3120 NW_019211460.1_36427509 0/0
## 3121 NW_019211460.1_36440396 0/0
## 3122 NW_019211460.1_36495785 0/0
## 3123 NW_019211460.1_36495846 0/0
## 3124 NW_019211460.1_36496401 0/0
## 3125 NW_019211460.1_36505988 <NA>
## 3126 NW_019211460.1_36506890 0/0
## 3127 NW_019211460.1_36506925 0/0
## 3128 NW_019211460.1_36506986 0/0
## 3129 NW_019211460.1_36530456 0/0
## 3130 NW_019211460.1_36537102 0/0
## 3131 NW_019211460.1_36540823 0/0
## 3132 NW_019211460.1_36574100 0/0
## 3133 NW_019211460.1_36575932 <NA>
## 3134 NW_019211460.1_36583026 <NA>
## 3135 NW_019211460.1_36611771 0/0
## 3136 NW_019211460.1_36620185 0/0
## 3137 NW_019211460.1_36659362 <NA>
## 3138 NW_019211460.1_36682242 0/0
## 3139 NW_019211460.1_36701049 0/0
## 3140 NW_019211460.1_36719644 <NA>
## 3141 NW_019211460.1_36734940 <NA>
## 3142 NW_019211460.1_36735078 0/0
## 3143 NW_019211460.1_36741319 <NA>
## 3144 NW_019211460.1_36742951 0/0
## 3145 NW_019211460.1_36772888 0/0
## 3146 NW_019211460.1_36805687 <NA>
## 3147 NW_019211460.1_36826968 0/0
## 3148 NW_019211460.1_36827037 0/0
## 3149 NW_019211460.1_36831567 <NA>
## 3150 NW_019211460.1_36840509 0/0
## 3151 NW_019211460.1_36840579 <NA>
## 3152 NW_019211460.1_36845653 0/0
## 3153 NW_019211460.1_36850173 <NA>
## 3154 NW_019211460.1_36857838 0/0
## 3155 NW_019211460.1_36865521 <NA>
## 3156 NW_019211460.1_36871446 <NA>
## 3157 NW_019211460.1_36890857 <NA>
## 3158 NW_019211460.1_36932253 0/0
## 3159 NW_019211460.1_36960392 <NA>
## 3160 NW_019211460.1_36961695 <NA>
## 3161 NW_019211460.1_36963691 <NA>
## 3162 NW_019211460.1_36983052 0/0
## 3163 NW_019211460.1_37022696 0/0
## 3164 NW_019211460.1_37025392 0/0
## 3165 NW_019211460.1_37027355 <NA>
## 3166 NW_019211460.1_37030423 <NA>
## 3167 NW_019211460.1_37033340 0/0
## 3168 NW_019211460.1_37033695 0/0
## 3169 NW_019211460.1_38389758 0/0
## 3170 NW_019211460.1_38423641 0/0
## 3171 NW_019211460.1_38425397 <NA>
## 3172 NW_019211460.1_38501125 <NA>
## 3173 NW_019211460.1_38502540 <NA>
## 3174 NW_019211460.1_38534624 0/0
## 3175 NW_019211460.1_38581841 <NA>
## 3176 NW_019211460.1_38670263 0/0
## 3177 NW_019211460.1_38670283 0/0
## 3178 NW_019211460.1_38670376 0/0
## 3179 NW_019211460.1_38706950 <NA>
## 3180 NW_019211460.1_38721194 <NA>
## 3181 NW_019211460.1_38721858 0/0
## 3182 NW_019211460.1_38730978 <NA>
## 3183 NW_019211460.1_38757613 0/0
## 3184 NW_019211460.1_38777645 <NA>
## 3185 NW_019211460.1_38787839 <NA>
## 3186 NW_019211460.1_38858664 0/0
## 3187 NW_019211460.1_38860985 <NA>
## 3188 NW_019211460.1_38862382 0/1
## 3189 NW_019211460.1_39010532 <NA>
## 3190 NW_019211460.1_39017021 <NA>
## 3191 NW_019211460.1_39017101 <NA>
## 3192 NW_019211460.1_39026115 <NA>
## 3193 NW_019211460.1_39044960 0/0
## 3194 NW_019211460.1_39045225 <NA>
## 3195 NW_019211460.1_39076063 0/0
## 3196 NW_019211460.1_39150362 <NA>
## 3197 NW_019211460.1_39155077 0/0
## 3198 NW_019211460.1_39185495 0/0
## 3199 NW_019211460.1_39299544 <NA>
## 3200 NW_019211460.1_39301652 0/0
## 3201 NW_019211460.1_39380485 <NA>
## 3202 NW_019211460.1_39395023 <NA>
# join all grandkids of 26 families together
grnson0.0_0.0 <- reduce(list_grnson, full_join, by = "site") %>% view()
grndat0.0_0.0 <- reduce(list_grndat, full_join, by = "site") %>% view()
# it looks like the 'false sites' in one family, are not necessarily the false sites in other families. should i take them out?